#! /bin/sh

#
# GIF Script : gempak_gif2.sh
#
# Log :
# D.W.Plummer/NCEP   2/97   Add log header
# D.W.Plummer/NCEP  12/97   Added NGM v. ETA plot comparison if NGM precedes ETA
# J. Carr/HPC        8/98   Changed map to medium resolution
# J. Carr/HPC        1/99   Changed contur from 1 to 2 per OM request
# J. Carr/HPC        2/99   Changed skip to 0 per OM request
# B. Gordon/NCO      5/00   Ported to IBM-SP, Standardized for production,
#                           changed gdplot_nc -> gdplot2_nc
# B. Gordon/NCO      4/01   Modified to make individual GIF images.
#

set -xa

export PS4='GIF2:$MAPAREA$SECONDS + '

PMSL=pmsl
if [ $MODEL = "RUC2" ]
then
    PMSL=mmsl
fi

LATVAL="1/1/1/1/20;20"
if [ $MAPAREA = "namer" ]
then
   LATVAL=0
fi

cd $DATA
if [ ! -d $DATA/gif2 ] ; then
  mkdir $DATA/gif2
fi
cd $DATA/gif2
sh $utilscript/setup.sh

cp $FIXgempak/datatype.tbl datatype.tbl
#
# Need the next file in the current corking directory to create
# GIF images with a white background.
#
cp $FIXgempak/coltbl.xwp.wbg coltbl.xwp

#
# The 44 km Eta grid does not output any data at 250 mb, change 250 plot to 200.
#
if [ "$RUN" = "nam" ] ; then
  lev250=200
else
  lev250=250
fi

export pgm=gdplot2_gif;. prep_step; startmsg

  $GEMEXE/gdplot2_gif << EOF
\$MAPFIL= mepowo.gsf
GDFILE  = F-${MODEL2} | ${PDY2}/${cyc}00                 
GDATTIM	= F${fhr}
PANEL	= 0
TEXT	= .8/2/2/c/sw
CONTUR	= 2
MAP	= 8     
CLEAR	= yes
CLRBAR  = 1|.8/2/2/sw

GAREA	= $GAREA 
PROJ	= $PROJ
LATLON	= 0

restore ${USHgempak}/restore/500mb_hght_absv.2.gif.nts
DEVICE  = gif | ${RUN}_500_${fhr}m.gif | 1024 ; 768
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
LATLON  = 0
r
TITLE   = 1/2/~ ${MODEL} 500 MB HGT, GEO ABS VORTICITY!0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
LATLON  = $LATVAL
r

restore ${USHgempak}/restore/250mb_hght_wnd.2.gif.nts
GLEVEL  = $lev250
DEVICE  = gif | ${RUN}_${lev250}_${fhr}m.gif | 1024 ; 768
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes 
LATLON  = 0
r
TITLE   = 1/2/~ ${MODEL} ${lev250} MB  HEIGHTS, ISOTACHS AND WIND (KTS)!0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
LATLON  = $LATVAL
r

restore ${USHgempak}/restore/250mb_hght_wnd.2.gif.nts
DEVICE  = gif | ${RUN}_300_${fhr}m.gif | 1024 ; 768
GLEVEL  = 300
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
LATLON  = 0
r
TITLE   = 1/2/~ ${MODEL} 300 MB  HEIGHTS, ISOTACHS AND WIND (KTS)!0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
LATLON  = $LATVAL
r


restore ${USHgempak}/restore/p06m_pmsl.2.gif.nts
gdpfun  = p${precip_inc}i !sm9s($PMSL)
DEVICE  = gif | ${RUN}_p${prec_name}_${fhr}m.gif | 1024 ; 768
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
LATLON  = 0
r
TITLE   =
GDPFUN  = sm9s(tmpc) ! sm9s(tmpc) ! sm9s(tmpc)
GVCORD  = pres
type    = c          ! c          ! c
CINT    = 10/-90/1   ! 10 /10/20  ! 10 /30/90
LINE    = 4/1/2/1    ! 17/1/2/1   ! 2/1/2/1
GLEVEL  = 850
CLEAR   = no
r
TITLE   = 1/2/~ ${MODEL} ${precip_inc}-HOUR TOTAL PCPN, MSLP, 850 MB TEMPS!0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
LATLON  = $LATVAL
r

EOF

####################################
# NOT NECESSARY WHEN RUNNING JUST GDPLOT2_GF
####################################
#$GEMEXE/gpend

#####################################################
# GEMPAK DOES NOT ALWAYS HAVE A NON ZERO RETURN CODE
# WHEN IT CAN NOT PRODUCE THE DESIRED GRID.  CHECK
# FOR THIS CASE HERE.
#####################################################
ls -l *${fhr}*.gif
export err=$?
if [ $err -ne 0 ] ; then
  ps -ef
  cp /tmp/xlogfile $DATA/
fi
export pgm="GEMPAK CHECK FILE";err_chk

mv *${fhr}*.gif $DATA

#