#!/bin/sh
set -e
# Ocean Model Cell Thickness - use missing values over land
#grido gz BETA DA DX DY DZ, BETA*DZ,

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

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

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

grido desc beta da dx dy dz
$RTEXTBLS/lopgm_o4xa/mlt dz beta thk
echo "C*FMSKPLT         -1 NEXT   GT      0.00    1     1.E38    1" | ccc fmskplt thk ${pfx}_${out}

release desc beta da dx dy dz thk