#!/bin/ksh
# Plot a precip analysis.
# /nwprod/gempak/.gempak was executed in the calling routine
# (so we are not running it unnecessarily when making multiple plots).
#
set -x

GBFILE=$1
GFFILE=$2
DATE=$3
AC=$4
TITLE="$5 ${AC}h Accum (mm) Ending $DATE"

# GDPFUN needs to be p01m/p03m/p24m, etc., so pad $AC with leading zero, if
# necessary.
typeset -R2 -Z AC
 
GDFILE=${GBFILE}.grd
 
# If for any reason - e.g. job hiccups - the gif file already exists,
# GEMPAK will produce something like "file.gif.2".  We don't want that,
# so check for existing gif.
 
if [ -e $GFFILE ]; then
  rm -f $GFFILE
fi

nagrib << EOFgrib
 GBFILE   = $GBFILE
 INDXFL   =
 GDOUTF   = $GDFILE
 PROJ     =
 GRDAREA  =
 KXKY     =
 MAXGRD   = 4999
 CPYFIL   = GDS
 GAREA    = dset
 OUTPUT   = T
 GBTBLS   = wmogrib.tbl
 GBDIAG   =
 PDSEXT   = NO
 
r
 
ex
 
EOFgrib
 
$GEMEXE/gdplot2 <<EOFplot
  GDFILE   = $GDFILE
  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

exit