count=1
while (count<21)
if (count=1)
varshort="SoilMst1"
variable="ETA Soil Mst 0-10 cm"
correction="*100"
correction2="1"
units="(mm)"
 colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
 barunits="0 6 12 18 24 30 36 42 48 54 60 66 72"
endif

if (count=2)
varshort="SoilMst2"
variable="ETA Soil Mst 10-40 cm"
correction="*300"
correction2="1"
units="(mm)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 15 30 45 60 75 90 105 120 135 150 165 180"
endif

if (count=3)
varshort="SoilMst3"
variable="ETA Soil Mst 40-100 cm"
correction="*600"
correction2="1"
units="(mm)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 30 60 90 120 150 180 210 240 270 300 330 360"
endif

if (count=4)
varshort="SoilMst4"
variable="ETA Soil Mst 100-200 cm"
correction="*1000"
correction2="1"
units="(mm)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 60 120 180 240 300 360 420 480 540 600 660 720"
endif

if (count=5)
varshort="SoilMstTot"
variable="ETA Total Column Soil Moisture\0-200 cm Layer,"
correction="*1"
correction2="1"
units="(mm)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="100 200 300 400 500 600 700 800 900 1000 1100 1200 1300"
endif

if (count=6)
varshort="SensHeatFlux"
variable="ETA Sensible Heat Flux"
correction="*1"
correction2="1"
units="(W/m`a2`n)"
barunits="-400 -360 -320 -280 -240 -200 -160 -120 -80 -40 0 40 80"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=7)
varshort="LatentHeatFlux"
variable="ETA Latent Heat Flux"
correction="*1"
correction2="1"
units="(W/m`a2`n)"
barunits="-330 -300 -270 -240 -210 -180 -150 -120 -90 -60 -30 0 30"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=8)
varshort="SoilTmp1"
variable="ETA Soil Temperature\0-10 cm Layer,"
correction="-273.15"
correction2="1"
units="(`3.`0C)"
barunits="-25 -20 -15 -10 -5 0 5 10 15 20 25 30 35"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=9)
varshort="SoilTmp2"
variable="ETA Soil Temperature\10-40 cm Layer,"
correction="-273.15"
correction2="1"
units="(`3.`0C)"
barunits="-25 -20 -15 -10 -5 0 5 10 15 20 25 30 35"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=10)
varshort="SoilTmp3"
variable="ETA Soil Temperature\40-100 cm Layer,"
correction="-273.15"
correction2="1"
units="(`3.`0C)"
barunits="-25 -20 -15 -10 -5 0 5 10 15 20 25 30 35"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=11)
varshort="SoilTmp4"
variable="ETA Soil Temperature\100-200 cm Layer,"
correction="-273.15"
correction2="1"
units="(`3.`0C)"
barunits="-25 -20 -15 -10 -5 0 5 10 15 20 25 30 35"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=12)
varshort="SnowDepth"
variable="ETA Liquid Equiv. Snow Depth"
correction="*1"
correction2="1"
units="(mm)"
barunits="0 7 14 21 28 35 42 49 56 63 70 77 84"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=13)
varshort="SoilMst1"
variable="ETA Total Runoff"
correction="*10800"
correction2="1"
units="(mm over 3 hours)"
barunits=".5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=14)
varshort="GroundHeatFlux"
variable="ETA Ground Heat Flux"
correction="*1"
correction2="1"
units="(W/m`a2`n)"
barunits="-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

if (count=15)
varshort="SoilMst1"
variable="ETA Volumetric Soil Mst\0-10 cm Layer,"
correction="*1"
correction2="1"
units="(m`a3`n/m`a3`n)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 .05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"
endif

if (count=16)
varshort="SoilMst2"
variable="ETA Volumetric Soil Mst\10-40 cm Layer,"
correction="*1"
correction2="1"
units="(m`a3`n/m`a3`n)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 .05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"
endif

if (count=17)
varshort="SoilMst3"
variable="ETA Volumetric Soil Mst\40-100 cm Layer,"
correction="*1"
correction2="1"
units="(m`a3`n/m`a3`n)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 .05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"
endif

if (count=18)
varshort="SoilMst4"
variable="ETA Volumetric Soil Mst\100-200 cm Layer,"
correction="*1"
correction2="1"
units="(m`a3`n/m`a3`n)"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 .05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"
endif

if (count=19)
varshort="TotColRel"
variable="ETA Relative Soil Mst\0-100 cm Layer,"
correction="*1"
correction2="1"
units="(%)"
colors="9 4 11 5 13 3 10 7 12 8 2 6"
barunits="0 10 20 30 40 50 60 70 80 90 100"
endif

if (count=20)
varshort="Skintemp"
variable="ETA Skin Temperature,"
correction="-273.15"
correction2="1"
units="(`3.`0C)"
barunits="-15 -10 -5 0 5 10 15 20 25 30 35 40 45"
colors="15 9 14 4 11 5 13 3 10 7 12 8 2 6"
endif

t=1
while (t<9)
if (count<13 | count=14)
"enable print /usr/local/www/LDAS/LDASpredictions/wrkdir/"varshort"."t".eta.raw"
endif

if (count=13)
"enable print /usr/local/www/LDAS/LDASpredictions/wrkdir/TotalRunoff."t".eta.raw"
endif

if (count>14 & count<19)
"enable print /usr/local/www/LDAS/LDASpredictions/wrkdir/"varshort"Vol."t".eta.raw"
endif

if (count>18)
"enable print /usr/local/www/LDAS/LDASpredictions/wrkdir/"varshort"."t".eta.raw"
endif


"set t "t
 "set grads off"
 var2=varshort           "(t=          "t")"
"var3=("var2 correction")"
"var=("var3"*"correction2")"

 "set mpdset hires"
 "set gxout grfill"
 "set clevs "barunits
 "set ccols "colors

if (variable = 'ETA Total Runoff')
"display (StormSfcRunoff(t="t")+GrndWaterRunoff(t="t"))"
else
"display var"
endif
if (t='1')
time='15Z'
"draw title "variable " "units " at "time " on "startdate
endif
if (t='2')
time='18Z'
"draw title "variable " "units " at "time " on "startdate
endif
if (t='3')
time='21Z'
"draw title "variable " "units " at "time " on "startdate
endif
if (t='4')
time='00Z'
"draw title "variable " "units " at "time " on "startdate
endif
if (t='5')
time='03Z'
"draw title "variable " "units " at "time " on "enddate1
endif
if (t='6')
time='06Z'
"draw title "variable " "units " at "time " on "enddate1
endif
if (t='7')
time='09Z'
"draw title "variable " "units " at "time " on "enddate1
endif

if (t='8')
time='12Z'
"draw title "variable " "units " at "time " on "enddate1
endif

"run /usr/local/www/LDAS/LDASpredictions/cbarn.gs"
 var2=varshort           "(t=          "t")"
"var3=("var2 correction")"
"var=("var3"*"correction2")"
 "set clevs "barunits
"set gxout contour"
"set clab off"
"set ccols 1 1 1 1 1 1 1 1 1 1 1 1 1 1"
if (variable = 'ETA Total Runoff')
"display (StormSfcRunoff(t="t")+GrndWaterRunoff(t="t"))"
else
"display var"
endif
"set strsiz .12 .12"
"draw string 9.65 .2 ldas.gsfc.nasa.gov"
"draw string 1.41 .2 LDAS Images"
"print"
"clear"
if (count<13 | count=14)
"disable print /usr/local/www/LDAS/LDASpredictions/wrkdir/"varshort"."t".eta.raw"
endif

if (count=13)
"disable print /usr/local/www/LDAS/LDASpredictions/wrkdir/TotalRunoff."t".eta.raw"
endif

if (count>14 & count<19)
"disable print /usr/local/www/LDAS/LDASpredictions/wrkdir/"varshort"Vol."t".eta.raw"
endif

if (count>18)
"disable print /usr/local/www/LDAS/LDASpredictions/wrkdir/"varshort"."t".eta.raw"
endif

t=t+1
endwhile
count=count+1
endwhile
"quit"