#!/bin/sh
set -e
# compute surface altitude, m 
# ryj: CMAM doesn't have MASK in data_description

if [ $# -ne 4 ] ; then
  echo "Error in $0: invalid number of arguments $#"
  exit 1
fi

pfx=$1
phis=$2
gc=$3
out=$4

release phis gc
access  phis ${pfx}_${phis}
access  gc   ${pfx}_${gc}

echo "C*RCOPY            1         1" | ccc rcopy phis phis1
echo "C*XLIN       9.80616        0.         1" | ccc xlin phis1 zs

printf "C* RCOPY           1         1" | ccc rcopy gc gc1
printf "C* FMSKPLT        -1 NEXT   NE      1.00    1        0.    1    1" | ccc fmskplt gc1 gc2
printf "C* FMSKPLT        -1 NEXT   LT     -0.50    1        0.    1    1" | ccc fmskplt gc2 gc3
printf "C* RELABL\nC* RELABL               0 MASK" | ccc relabl gc3 mask

# set to 0 over oceans, lakes are unchanged
echo "C* FMSKPLT        -1 NEXT   NE      0.00    1     0.E00    1    1" | ccc fmskplt zs ${pfx}_${out} mask
# # this set 0 over lakes as well
# echo "C* FMSKPLT        -1 NEXT   EQ      -1.0    1     0.E00    1    1" | ccc fmskplt zs ${pfx}_${out} mask

release phis gc phis1 zs mask gc1 gc2 gc3