#! /bin/sh

#
# GIF Script : gempak_gif1.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.
# L. Sager/NCO       5/05   Modified to make Eastern Pacific graphics
#

set -xa

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

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

cd $DATA
if [ ! -d $DATA/gif1 ] ; then 
  mkdir $DATA/gif1
fi
cd $DATA/gif1
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

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

restore ${USHgempak}/restore/pmsl_thkn.2.gif.nts
DEVICE  = gif | ${RUN}_slp_${fhr}m.gif | 1024 ; 768
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
GDPFUN  = p${precip_inc}i    !ldf(hght)      !ldf(hght)      !sm9s(pmsl)
LATLON  = 0
r
TITLE   = 1/2/~ ${MODEL} MSLP ,PRECIP, 1000-500 MB THICKNESS !0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
l
LATLON  = $LATVAL
r
  
restore ${USHgempak}/restore/850mb_hght_tmpc.2.gif.nts
DEVICE  = gif | ${RUN}_850_${fhr}m.gif | 1024 ; 768
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
LATLON  = 0
r
TITLE   = 1/2/~ ${MODEL} 850 HGT, TEMP AND WIND (KTS)|~@ HGT, TMP, WIND!0!0!0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
LATLON  = $LATVAL
r

restore ${USHgempak}/restore/700mb_hght_relh_omeg.2.gif.nts
DEVICE  = gif | ${RUN}_700_${fhr}m.gif | 1024 ; 768
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
LATLON  = 0
r
TITLE   = 1/2/~ ${MODEL} 700 MB HGT, REL HUMIDITY AND OMEGA|~@ HGT, RH AND OMEGA!0
TEXT    =.8/2/2/c/sw
GDPFUN  =
CLEAR   =no
LATLON  = $LATVAL
r

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

restore ${USHgempak}/restore/850mb_vort_winds.2.gif.nts    
DEVICE  = gif | ${RUN}_85v_${fhr}m.gif |  1024 ; 768 
TITLE   = 1/0/ $TITLEIT !0
CLEAR   = yes
LATLON  = 0
l
r
TITLE   = 1/2/~ GFS 850MB HGTS(M) VORTICITY(SEC-1) AND WINDS(KNOTS)  !0 
TEXT    =.8/2/2/c/sw
GDPFUN  = 
CLEAR   = no
LATLON  = $LATVAL   

r

restore ${USHgempak}/restore/pmsl_10mw_p06m.2.gif.nts
DEVICE  = gif | ${RUN}_ten_${fhr}m.gif | 1024 ; 768
GDPFUN  = p${precip_inc}i       !sm9s(pmsl)         ! kntv(wnd)
TITLE   = 1/0/ $TITLEIT !0
LATLON  = 0
r

CLEAR   = no
TITLE   = 1/2/~ GFS   MSLP 10M WIND (KNOTS) 2M TEMP (C) !0
TEXT    = .8/2/2/c/sw
LATLON  = $LATVAL
r

CLEAR   = no
TITLE   = !
GDPFUN  = sm9s(tmpc) ! sm9s(tmpc)
TYPE    = c
GLEVEL  = 2
GVCORD  = HGHT
CINT    = 10/-30/0   ! 10/10/30
LINE    = 4/1/2      !  2/1/2
LATLON  = 0
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

#