#!/bin/ksh
######################################################################
#  UNIX Script Documentation Block
#                      .
# Script name:         exaqm_prephyb_AK.sh.sms
# Script description:  Run nam product generator for CMAQ
#
# Abstract: This script runs 0-48 h NAM PRDGEN for CMAQ
#
# Script history log:
# 04-06-2010 Jianping Huang modified for AK
# 01-11-2011 Jianping Huang modified for NMMB-CMAQ coupling
######################################################################
set -xa
msg="JOB $job HAS BEGUN"
postmsg "$jlogfile" "$msg"

cd $DATA

echo "Start of NMMB_prdgen Job"

typeset -Z2 fhr
export fhr=00
if [ ${cycle} = 't00z' -o ${cycle} = 't18z' ]
then
   export endfhr=48
elif [  ${cycle} = 't06z' ]
then
   export endfhr=48
else
   export endfhr=48
fi

cp $PARMaqm/aqm_255_w60_140_nmmb.109tables.ctl master3.ctl

for fhr in $modak4

do
  ic=1
  while [ $ic -lt 1000 ]
  do
    if [ -s $COMNMM/nam.${cycle}.bgrd3d$fhr.tm00 ]
    then
        ln -s -f $COMNMM/nam.${cycle}.bgrd3d$fhr.tm00 nam.${cycle}.bgrd3d$fhr.tm00
   break
    else
          let "ic=ic+1"
          sleep 10
      fi

      if [ $ic -ge 180 ]
      then
      err_exit "COULD NOT LOCATE:$COMNMM/nam.${cycle}.bgrd3d$fhr.tm00"
      fi
    done
##===============================================
## using previous -24 hr fsct bgrd3d files as a proxy
##------------------------------------------------
case $cyc in
00) export flnm3=$COMNMMm1/nam.t00z.bgrd3d24.tm00;;
06) export flnm3=$COMNMMm1/nam.t06z.bgrd3d24.tm00;;
12) export flnm3=$COMNMMm1/nam.t12z.bgrd3d24.tm00;;
18) export flnm3=$COMNMMm1/nam.t18z.bgrd3d24.tm00;;
esac
 fsz1=0
 fsz2=0
 fc=0
 if [ $fhr = 00 ] ; then
 rm -rf nam.t${cyc}z.bgrd3d00.tm00
 cp -rp $COMNMM/nam.t${cyc}z.bgrd3d00.tm00 .
 cp -rp $COMNMM/nam.t${cyc}z.bgrd3d01.tm00 .
 export flnm1=nam.t${cyc}z.bgrd3d00.tm00
 export flnm2=nam.t${cyc}z.bgrd3d01.tm00
 fsz1=`ls -l $flnm1 | awk '{print $5}'`
 fsz2=`ls -l $flnm2 | awk '{print $5}'`
 ((fc=${fsz2}*7/10))
fi
if [ $fsz1 -lt $fc ] ; then
 cp -rp $flnm1 nam.t${cyc}z.bgrd3d00.tm00_proxy
 ln -sf $flnm3 fort.11
 ln -sf nam.t${cyc}z.bgrd3d00.tm00_proxy fort.51
cat > itag  << EOF
$PDY$cyc
-24
EOF
/nwprod/util/exec/overdate.dgexgrib < itag > errfile
cp -rp nam.t${cyc}z.bgrd3d00.tm00_proxy nam.t${cyc}z.bgrd3d00.tm00
fi
##=====================================================

  export pgm=aqm_nmm_prep
  . prep_step
  export XLFRTEOPTS="unit_vars=yes"
  export XLFUNIT_10="master3.ctl"
  export XLFUNIT_21="$FIXaqm/aqm_wgt_bgrd12_140"

  cat <<EOF5 >input${fhr}.prd
nam.${cycle}.bgrd3d${fhr}.tm00
EOF5

  echo 'about to cat input'
  cat input${fhr}.prd

  startmsg
 $EXECaqm/aqm_prep_nmmb < input${fhr}.prd >> $pgmout 2>errfile
  export err=$?;err_chk

  cp  $DATA/meso.AQFNMM${fhr} $DATA/aqm.${cycle}.nmm$fhr.tm00

  if [ "$SENDCOM" = 'YES' ]
  then
  cp $DATA/meso.AQFNMM${fhr} $COMOUT/aqm.${cycle}.nmm$fhr.tm00
  fi

  let "fhr=fhr+1"
  typeset -Z2 fhr

done

echo EXITING $0

########################################################

msg='ENDED NORMALLY.'
postmsg "$jlogfile" "$msg"

################## END OF SCRIPT #######################
