#!/bin/sh
set -e

# For the 15 reflectivity/height bins for CloudSat generated by the CloudSat simulator
# Mask out the missing points and combine them into a single variable with the level corresponding to a bin, i.e., ZE01 -> bin 1, etc..

if [ $# -ne 17 ] ; then
  echo "Error in $0: invalid number of arguments $#"
  exit 1
fi
pfx=$1
ze01=$2
ze02=$3
ze03=$4
ze04=$5
ze05=$6
ze06=$7
ze07=$8
ze08=$9
ze09=$10
ze10=$11
ze11=$12
ze12=$13
ze13=$14
ze14=$15
ze15=$16
out=$17

files=""
n=1
while [ $n -le 15 ] ; do
  n=`echo $n | awk '{printf "%02d", $1}'`

  eval v=\$sr$n

  access $v ${pfx}_$v
  echo "C* FMSKPLT        -1 NEXT   GT     -999.    1     1.E38    1" | ccc fmskplt $v m$v

  # Relabel the name and ibuf(7)
  echo "C*RELABL   GRID
C*         GRID             ZE               $n" | ccc relabl m$v f$n

  release $v m$v

  files="$files f$n"
  n=`expr $n + 1`
done

/home/esm/PRODUCTION/cmip5_scripts/rmergel ${pfx}_${out} $files

rm $files .ccc_cards