#!/bin/csh -f
setenv GEMVER `cat /unidata/gempak.latest`

cd /home/ralazear/scripts/trop

/bin/rm -rf /home/ralazear/scripts/trop/last.nts
/bin/rm -rf /home/ralazear/scripts/trop/gemglb.nts


set yy = `date -u +%y`
set mm = `date -u +%m`
set dd = `date -u +%d`
set hh = `date -u +%H`
set ddm = `date +%d`
 
if (${hh} < 04) then
    set mdate = ${yy}${mm}${ddm}18
    echo ${mdate}
endif
if (${hh} >= 22) then
    set mdate = ${yy}${mm}${ddm}18
    echo ${mdate}
endif
if (${hh} >= 04 && ${hh} < 10) then
    set mdate = ${yy}${mm}${dd}00
    echo ${mdate}
endif
if (${hh} >= 10 && ${hh} < 16) then
    set mdate = ${yy}${mm}${dd}06
    echo ${mdate}
endif
if (${hh} >= 16 && ${hh} < 22) then
    set mdate = ${yy}${mm}${dd}12
    echo ${mdate}
endif   



gdvint << EOF

gdoutf = /spare11/ralazear/gem/gfstrop_$1.gem
gdfile = /data7/nmc2/${mdate}_gfs_003.gem
latlon = 0
gdatt = F$1
gvco = pres/thta
glev = 260-400-5
maxgrd = 29999
garea = $2
run


exit


EOF


$GEMEXE/gpend


gdpvsf << EOG
gdfile = /spare11/ralazear/gem/gfstrop_$1.gem
gdoutf = /spare11/ralazear/gem/gfstrop_$1.gem
gfunc = mul(pvor(pres,obs),9.8e5)
gdattim = F$1
gvcord = THTA
startl = 400
stopl = 300
desire = 1.5
gdoutl = 15
gvoutc = TROPA
gpack = 
glist = uwnd;vwnd;thta;pres
pmax = 700


run


exit


EOG


$GEMEXE/gpend