#!/bin/sh
set -e

# For the reflectances at the 5 angles compute by the PARASOL simulator
# Mask out the missing points and combine into a single level with each level
# corresponding to a different angle, i.e., PL01 -> level 1, etc..

if [ $# -ne 7 ] ; then
  echo "Error in $0: invalid number of arguments $#"
  exit 1
fi
pfx=$1
pl01=$2
pl02=$3
pl03=$4
pl04=$5
pl05=$6
out=$7

vars="pl01 pl02 pl03 pl04 pl05"

vars_join=""
lev=1
for v in $vars ; do

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

  # Relabel the name and level
  echo "C*RELABL   GRID
C*         GRID             PL${lev}" | ccc relabl m_temp f$lev

  vars_join="$vars_join f$lev"

  release ${v}
  rm m_temp
  lev=`expr $lev + 1`
done

rmerge ${pfx}_${out} ${vars_join}

rm ${vars_join} .ccc_cards