#!/bin/sh
set -e
# compute surface altitude, m

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

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

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

echo "C*RCOPY            1         1" | ccc rcopy phis phis1
echo "C*RCOPY            1         1" | ccc rcopy mask mask1
square mask1 mask2
echo "C*XLIN       9.80616        0.         1" | ccc xlin phis1 zs
# set to 0 over oceans, lakes are unchanged
echo "C* FMSKPLT        -1 NEXT   GT      0.50    1     0.E00    1    1" | ccc fmskplt zs ${pfx}_${out} mask2
# # 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} mask1

release phis phis1 zs mask mask1 mask2 .ccc_cards