#!/bin/sh
#############################################################################
# Name of script:    verf_precip_plotpcp.sh
# Purpose of script: This script plots the 24h model precipitation using
#                    GEMPAK
# Arguments
#  1. 'model' : Model name
#  2. 'vdate' : Verification date/hour
#  3. 'vacc'  : Verification length
#  4: 'grids' : grids to be verified on
#  5: 'cycles': model cycles
#  6: 'bucket': bucket length
#  7: 'frange': range of model forecast
#  8: 'nest'  : special mask?  (nests)
#############################################################################
set -x

cd $DATA

model=$1
fhrbgn=$2
fhrend=$3
fhrinc=$4
grid=$5
vdate=$6
AC=$7

MODEL=`echo $model | tr '[a-z]' '[A-Z]'`
fhr=$fhrbgn

typeset -Z3 fhr fhrbgn fhrend fhr0

cp /nwprod/gempak/fix/coltbl.xwp.wbg coltbl.xwp
cp /nwprod/gempak/fix/wmogrib.tbl .

while [ $fhr -le $fhrend ]
do
   for grid in $grid
   do
      pdate=`/nwprod/util/exec/ndate -$fhr $vdate`
      fhr0=`expr $fhr - $AC`
      GRIBIN=$COMIN/${model}_${pdate}_${fhr0}_${fhr}.$grid
      GMDOUT=${model}_${pdate}_${fhr0}_${fhr}.${grid}.grd
      GFFILE=${model}.v${vdate}.${fhr}h.gif
      TITLE="$MODEL"

      rm -f $GFFILE

nagrib << EOFgrib
 GBFILE   = $GRIBIN
 INDXFL   =
 GDOUTF   = $GMDOUT
 PROJ     =
 GRDAREA  =
 KXKY     =
 MAXGRD   = 4999
 CPYFIL   = GDS
 GAREA    = dset
 OUTPUT   = T
 GBTBLS   = wmogrib.tbl
 GBDIAG   =
 PDSEXT   = NO
 
r
 
ex
 
EOFgrib
 
echo "Start to plot the precip data"

$GEMEXE/gdplot2 <<EOFplot
  GDFILE   = $GMDOUT
  GLEVEL   = 0
  GVCORD   = none
  GDPFUN   = p${AC}m
  PANEL    = 0/1/1/1
  DEVICE   = gif | $GFFILE | 480;425
  SCALE    = 0
  TYPE     = f
  CONTUR   = 0
  CINT     = 2
  LINE     = 2/1/1/1
  FINT     = 0.1;2;5;10;15;20;25;35;50;75;100;125;150;175
  FLINE    = 31;23;22;21;20;19;10;17;16;15;14;29;28;24;25
  HILO     =
  HLSYM    =
  CLRBAR   = 1/V/LL/.005;.02/.9;.013/|0.50/1
  GVECT    =
  REFVEC   =
  TITLE    = 1/0/$TITLE
  TEXT     = 0.8
  CLEAR    = y
  GAREA    = 18;-124;50;-60
  SKIP     = 0
  PROJ     = lcc/35;-95;35
  MAP      = 1/1/1
  LATLON   = 0
  STNPLT   =
 
r

EOFplot
$GEMEXE/gpend

   done
   fhr=`expr $fhr + $fhrinc`
done
exit