#!/bin/sh
set -e
# compute sea ice snow thickness from daily SNO as it is done in model

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

pfx=$1
snow=$2
sicn=$3
out=$4

mtype=`uname -s|tr '[A-Z]' '[a-z]'`
case $mtype in
    aix) BINDIR=~ccc101/gcm2nc/cmip5_scripts/bin_aix ;;
  linux) BINDIR=~ccc101/gcm2nc/cmip5_scripts/bin_linux ;;
      *) echo "Unkown kernel type $mtype"
         exit 1 ;;
esac
PATH="${BINDIR}:$PATH"
export PATH

release snow sicn
access snow ${pfx}_${snow}
access sicn ${pfx}_${sicn}

snowdens snow dens
div snow dens zn

echo "C* FMSKPLT        -1 NEXT   GT      0.00    1     0.E00    1    1" | ccc fmskplt zn zn1 sicn
echo "C*MONAVG    YYYYMMDD" | ccc monavg zn1 ${pfx}_${out}

release snow sicn dens zn zn1