#!/bin/sh
set -e
# compute surface altitude in m and then put it on latitude/longitude grid

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

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

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

echo "C*RCOPY            1         1" | ccc rcopy phis phis1
echo "C*XLIN       9.80616        0.         1" | ccc xlin phis1 zs
echo "C*SELECT.  STEP         0         0    1         1    1 NAME MASK" | ccc select desc 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 zs_mask mask

# put data onto a longitude/latitude grid with a 1 degree grid-spacing

echo "C* GGALL    360  180   -1    1       0.5       0.5    1" | ccc ggall zs_mask ${pfx}_${out}

release phis desc phis1 zs mask zs_mask