#!/bin/sh

#####################################################################
#
#  script for running forecasts for breeding cycle from reananlysis data.
#
# History: Jul 2004 - original author: Jeffrey Whitaker
#          Jul 2004 - Modified by Luke Lin for the IBM SP"
#
#  environment variables exported to this script:
#    FCSTDIR
#    PERT
#
#    SIGI
#    SFCI
#    SIGO
#    SFCO
#    FHMAX
#    FHOUT
#    IGEN
#    EXECDIR
#    SIGR1
#    SIGR2
#    SFCR
#    DELTIM
#    FHRES
#    FHZER
#    CO2CON
#    MTNVAR
#    TUNE1
#####################################################################

set -x
cd $FCSTDIR

#  Input files.
co2const=$CO2CON
mtnvar=$MTNVAR
tune1=$TUNE1
SIGANL=$SIGI
SFCANL=$SFCI

#  Create namelist parameter file that is used by the model.
FHOUR=$FHMAX

rm -f fcstparms
cat <<EOF1 >fcstparms
 &NAMGSM
  FHMAX=$FHOUR,FHRES=$FHRES,FHOUT=$FHOUT,FHZER=$FHZER,
  DELTIM=$DELTIM,IGEN=$IGEN,NDIGYR=4
 &END
EOF1


export pgm=cdas_fcst6228
. prep_step

#  Symbolic links for FORTRAN logical units used in the model.
ln -s $SIGANL   fort.11
if [[ -f $SIGI2 ]]; then
  ln -s "$SIGI2"  fort.12
else
  ln -s /dev/null fort.12
fi
ln -s $SFCANL   fort.13
ln -s $co2const fort.15
ln -s $mtnvar   fort.24
ln -s $tune1    fort.43
ln -s $SIGR1    fort.51
ln -s $SIGR2    fort.52
ln -s $SFCR     fort.53

#  Run the executable

startmsg
$EXECGLOBAL/cdas_fcst6228 < fcstparms  >> $pgmout 2> errfile
err=$?;export err; err_chk

echo `date` >! $FCSTDIR/$run_completed 

if test $SENDCOM = 'YES'
then
   cd $FCSTDIR
   fhr=$FHOUT
   if test $fhr -lt 10
   then
      fhr="0$fhr"
   fi
   while test $fhr -le $FHOUR
   do
      /bin/rm -f fort.111 fort.112 fort.151
      ln -s sig.f${fhr} fort.111
      ln -s $SIGI2  fort.112
      ln -s pbin.f${fhr} fort.151
      ${EXECGLOBAL}/cdas_sig2grb
      cp $FCSTDIR/sig.f${fhr} $COMOUT/ens${PERT}.$cycle.sf${fhr}
      cp $FCSTDIR/pbin.f${fhr} $COMOUT/ens${PERT}.$cycle.pf${fhr}
      cp $FCSTDIR/sfc.f${fhr} $COMOUT/ens${PERT}.$cycle.bf${fhr}
      fhr=`expr $fhr + $FHOUT`
      if test $fhr -lt 10
      then
         fhr="0$fhr"
      fi 
   done
fi

exit