#!/bin/sh
set -e
# get surface temperature on land and ses-ice
# for 3hr: tslsi=(FARE(1)*GTT(1)+(FARE(2)*GTT(2)+FARE(6)*GTT(6))/(FARE(1)+FARE(2)+FARE(6))
# for daily: bind 8

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

pfx=$1
inp1=$2
inp2=$3
out=$4

release inp1 inp2
access  inp1 ${pfx}_${inp1}
access  inp2 ${pfx}_${inp2}

echo "C*SELLEV      1    1" > ic.sellev1
sellev inp1  gtt1   input=ic.sellev1
sellev inp2  fare1  input=ic.sellev1
mlt gtt1 fare1 gts1

echo "C*SELLEV      1    2" > ic.sellev2
sellev inp1  gtt2   input=ic.sellev2
sellev inp2  fare2  input=ic.sellev2
mlt gtt2 fare2 gts2

echo "C*SELLEV      1    6" > ic.sellev6
sellev inp1  gtt6   input=ic.sellev6
sellev inp2  fare6  input=ic.sellev6
mlt gtt6 fare6 gts6

add gts1 gts2  gts12
add gts6 gts12 gts126

add fare1 fare2  fare12
add fare6 fare12 fare126

div gts126 fare126 gts3hr

#bind 3hr data into daily data
echo "C* BINSML     8"|ccc binsml gts3hr gtsday

delt=900; yearoff=0001
deltmin=`echo "$delt"         | $AWK '{printf "%5d", $0/60.}'`
echo "C*TSTEP   $deltmin     ${yearoff}010100          MODEL       YYYYMMDD" > ic.tstep_model
tstep gtsday ${pfx}_${out} input=ic.tstep_model

release inp1 inp2 gtt* fare* gts*