#!/bin/csh -f

set day   = `date -u +%d`
set CURFRT = `\ls -t /data/wxp/surface/front/*frt.wmo | head -1`
setenv DISPLAY :1.0
source /opt2/nawips/Gemenviron

@ THETA = 285

rm /scratch/mjames/isenUS${THETA}.gif

set DEV = 'gf|/scratch/mjames/isenustemp.gif|800;600|||NP'

   while($THETA <= 320)

   # Set map background to white
   /opt2/nawips/bin/sol/gpcolor <<EOF0
      color    = 101=black
      device = $DEV
      run

   exit
EOF0

   /opt2/nawips/bin/sol/gdplot << EOF2
   GDFILE= $HRS/cur_isen
   GDATTIM=last
   GLEVEL=$THETA
   GVCORD=thta
   PANEL=0
   SKIP=/-1;1 
   SCALE=999 
   GFUNC=pres 
   CTYPE=C  
   CONTUR= 
   CINT=50
   LINE=2/1/2 
   FINT=0 
   FLINE=2-12 
   HLSYM=
   CLRBAR=
   GVECT=WND 
   WIND=BK20/1\
   REFVEC=
   TITLE=31/-4/${THETA}K Theta Surface: ~
   TEXT=1/2//hw
   CLEAR=y
   AREA = us-
   GAREA = us
   PROJ=str/45;-100;0
   MAP=26/1
   LATLON=0
   DEVICE = $DEV
   STNPLT=
   SATFIL=
   RADFIL=
   run

   
   e  
EOF2

/opt2/nawips/bin/sol/gdcntr <<EOF1
PANEL   = 0
PROJ    =
GDATTIM = last
GVCORD  = thta
GDFILE  = $HRS/cur_isen
DEVICE  = $DEV
GLEVEL  = $THETA
GFUNC   = relh
CINT    = 20
LINE    = 3/1/2
\$MAPFIL=hipowo.cia
MAP     = 26/1
AREA = us-
GAREA = us
TITLE=31/-2/GREEN: Relative Humidity   RED: Pressure
PROJ=str/45;-100;0
CLEAR   = n
SCALE   = 0
LATLON  = 0
HLSYM   =
CLRBAR  = 32
CONTUR  = 2
SKIP    = 0
FINT    =
FLINE   =
CTYPE   = C
run
#this next blank line is VERY important
#it acts as an 'enter' to gdcntr

exit
EOF1



  /opt2/nawips/bin/sol/gpend

  cp /scratch/mjames/isenustemp.gif /scratch/mjames/isenUS${THETA}.gif
  rm /scratch/mjames/isenustemp.gif

  /usr/local/bin/combine -geometry +10+509 /scratch/mjames/isenUS${THETA}.gif /home/mjames/scripts/som.gif -compose Over /scratch/mjames/isenUS${THETA}.gif
  /usr/local/bin/convert /scratch/mjames/isenUS${THETA}.gif /scratch/mjames/isenUS${THETA}.png

  @ THETA = $THETA + 5

end

rm /home/mjames/last.nts
rm /home/mjames/gemglb.nts
exit