#!/bin/sh
#######################################################################
#  Job Name: exgfsmos_pac_extfcst.sh.sms
#  Purpose: To run all steps necessary to create extended range GFS MOS 
#           fcsts for the Pacific sites.
#  Remarks: 
#  HISTORY: May 24, 2005      - new job for GFS gridded MOS
#                               This job currently runs at all
#                               4 cycles, but just archives model data.
#
#######################################################################
#
echo MDLLOG: `date` - Begin job exgfsmos_pac_extfcst.sh.sms
set -x

cd $DATA

echo $PDY $cyc: Date and Cycle - echo PDY and cyc

export DAT="$PDY$cyc"

################################################################################
#  FIRST GATHER MODEL DATA AND PUT IN TDLPACK
#  LOOP THROUGH PROJECTIONS (00 THROUGH 96 ARE DONE IN GFSMOS_PAC_FCST)
################################################################################
for tau in 99 102 105 108 111\
           114 117 120 123 126 129 132 135 138 141\
           144 147 150 153 156 159 162 165 168 171\
           174 177 180 192 204 216 228 240 252 264\
           276 288 300 312 324 336 348 360 372 384
do
################################################################################
# COPY MODEL FILES TO TEMP SPACE
################################################################################

  cp $COMINGFS/gfs.$cycle.pgrbf$tau gfs.$DAT.pgrbf$tau
  cp $COMINGFS/gfs.$cycle.pgrbif$tau gfs.$DAT.pgrbif$tau

################################################################################
# RUN WGRIB TO GET INPUT LIST FOR COPYGB
# RUN COPYGB FOR NON PRECIP
# CONVERTS TO TDL GRID
################################################################################

  grid=`cat $FIXmdl/mdl_merc.finegds`

  g1=gfs.$DAT.pgrbf$tau
  x1=gfs.$DAT.pgrbif$tau
 
  $EXECutil/wgrib $g1 > wgrib.out
  grep -f $FIXmdl/mdl_pacgfs_xpcp.wgrib wgrib.out |\
  $EXECutil/copygb -a -kw -g"$grid" -i0 $g1 $x1 mdl_pacextxpcp.$cycle.pgrb 
 
################################################################################
# RUN WGRIB TO GET INPUT LIST FOR COPYGB
# RUN COPYGB FOR PRECIP
# CONVERTS TO TDL GRID
################################################################################

  grid=`cat $FIXmdl/mdl_merc.finegds`

  g1=gfs.$DAT.pgrbf$tau
  x1=gfs.$DAT.pgrbif$tau
 
  grep -f $FIXmdl/mdl_pacgfs_pcp.wgrib wgrib.out |\
  $EXECutil/copygb -a -kw -g"$grid" -i3 $g1 $x1 mdl_pacextpcp.$cycle.pgrb
 
done

################################################################################
# GET INDICES AND INVENTORY
################################################################################
$EXECutil/grbindex mdl_pacextxpcp.$cycle.pgrb mdl_pacextxpcp.$cycle.pgrbi
$EXECutil/grbindex mdl_pacextpcp.$cycle.pgrb mdl_pacextpcp.$cycle.pgrbi
$EXECutil/invindex mdl_pacextxpcp.$cycle.pgrbi > mdl_pacextxpcp.pgrbinv.$DAT
$EXECutil/invindex mdl_pacextpcp.$cycle.pgrbi > mdl_pacextpcp.pgrbinv.$DAT

################################################################################
# GRB2MDLP
# CONVERT GRIB TO TDL_PACK
################################################################################

export pgm=mdl_grb2mdlp
. prep_step

export XLFUNIT_10="ncepdate"
export XLFUNIT_11="mdl_pacextxpcp.$cycle.pgrb"
export XLFUNIT_21="mdl_pacextxpcp.$cycle.pgrbi"
export XLFUNIT_12="mdl_pacextpcp.$cycle.pgrb"
export XLFUNIT_22="mdl_pacextpcp.$cycle.pgrbi"
export XLFUNIT_30="$FIXmdl/mdl_gridlst"
export XLFUNIT_51="pac_pkgfsxraw.$DAT"
export XLFUNIT_52="grb2mdlp.lst"
echo MDLLOG: `date` - Program grb2mdlp has begun.
$EXECmdl/mdl_grb2mdlp < $PARMmdl/mdl_pacgfsx_arch.cn.$cycle >> $pgmout 2>errfile
export err=$?

#######################################################################
# COPY FILES TO COM
#######################################################################

if test $SENDCOM = 'YES'
then
  cp pac_pkgfsxraw.$DAT $COMOUT/mdl_pacgfsxpkd.$cycle
fi

echo MDLLOG: `date` - Job gfsmos_pac_extfcst has ended.
#######################################################################