#!/bin/sh
set -e
# compute monthly mean for SICN>0

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

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

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

echo "C*FMASK           -1 NEXT   GT      0.15" | ccc fmask sicn mask

mlt sicn mask sicm
mlt inp  sicm inps

#echo "C*MONAVG    YYYYMMDD" | ccc monavg inps inps_avg
#echo "C*MONAVG    YYYYMMDD" | ccc monavg sicm sicn_avg

div inps sicm inp_new

echo "C* FMSKPLT        -1 NEXT   GT       0.0    1     1.E38    1" | ccc fmskplt inp_new ${pfx}_${out} sicm

release inp sicn sicm mask inps inps_avg sicn_avg inp_avg