#!/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_avg sicn_avg inp_avg

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

release inp sicn sicm mask inps inps_avg sicn_avg inp_avg