#!/bin/sh
#
# args sanl sfcanl endhour [inchour]
#


if [ $# -lt 3 -o $# -gt 4 ] ; then
  echo "runs model: $0 sanl sfcanl endhour [inchour]"
  exit 6
fi

cd $DATA
set -x
date

sanl=$1
sfcanl=$2
endhour=$3

inchour=6
[ $# -eq 4 ] && inchour=$4

#Get Date from sigma analysis file
export pgm=cdas_sigmatime
. prep_step
export XLFUNIT_11="$sanl"
startmsg
sigdate=`$EXECcdas/cdas_sigmatime`
export err=$?;err_chk


# convert to double precision
$EXECcdas/cdas_sig2sig -toR8 $sanl sanl.dbl
$EXECcdas/cdas_sig2sig -toR8 $sfcanl sfcanl.dbl

hh=0
while [ $hh -lt $endhour ]
do
  hhp=`expr $hh + $inchour`

  export pgm=cdas_mrf;  . prep_step

  # input files
  if [ $hh -eq 0 ] ; then
    export XLFUNIT_11="sanl.dbl"
    export XLFUNIT_12="/dev/null"
    export XLFUNIT_14="sfcanl.dbl"
  else
    export XLFUNIT_11="sigft.${sigdate}ft$hh.dbl"
    export XLFUNIT_12="sigdt.${sigdate}ft$hh.dbl"
    export XLFUNIT_14="bges.${sigdate}ft$hh.dbl"
  fi

  export XLFUNIT_15="$FIXcdas/cdas_co2const.dbl"
  export XLFUNIT_21="heatrate"
  export XLFUNIT_24="$FIXcdas/cdas_mtnvar.6218.dbl"
  export XLFUNIT_43="$FIXcdas/cdas_tune1.season.dbl"

  # output files

  export XLFUNIT_51="sigft.${sigdate}ft$hhp.dbl"
  export XLFUNIT_52="sigdt.${sigdate}ft$hhp.dbl"
  export XLFUNIT_53="bges.${sigdate}ft$hhp.dbl"
  export XLFUNIT_61="znl.${sigdate}ift$hh"
  export XLFUNIT_62="flx.${sigdate}ift$hh"
  export XLFUNIT_63="flx.${sigdate}ft$hhp"
  export XLFUNIT_64="znl.${sigdate}ft$hhp"
  export XLFUNIT_65="dg3.${sigdate}ift$hh"
  export XLFUNIT_66="dg3.${sigdate}ft$hhp"
  export XLFUNIT_67="ken.${sigdate}ft$hhp"
  export XLFUNIT_68="flxs.${sigdate}ft$hhp"
  export XLFUNIT_69="znls.ft${sigdate}$hhp"
  export XLFUNIT_81="diabanl"
  export XLFUNIT_82="adiages"
  export XLFUNIT_83="fullges"
  export XLFUNIT_92="diagscr"
  export XLFUNIT_98="radscr"
  export XLFUNIT_99="w3out"

  startmsg
  $EXECcdas/cdas_mrf <<EOF >> $pgmout 2> errfile
 &NAMSMF ICEN2=1,NUM(5)=0,CON(7)=${inchour}.,CON(8)=3.,CON(1)=1200., /
EOF

  export err=$?; err_chk

  echo $hh
  hh=`expr $hh + $inchour`
done

#rm sfcanl.dbl sanl.dbl

if [ $endhour -eq 6 ] ; then
	flist="sigft.${sigdate}ft*.dbl sigdt.${sigdate}ft*.dbl bges.${sigdate}ft*.dbl"
else
	flist="sigft.${sigdate}ft*.dbl"
fi
for f in $flist 
do
   ff=`echo $f | sed 's/\.dbl//'`
   $EXECcdas/cdas_sig2sig -toR4 $f $ff
#   rm $f
done