#!/bin/csh 

#cd /spare11/ralazear/wrf_output


gdplot2 << endgdp

GDFILE   = /spare11/ktyle/wrf/$3_wrf03.gem
GDATTIM  = F$1
GLEVEL   = 0!!!!!10
GVCORD   = none!!pres!!none!hagl
PANEL    = 0
SKIP     = !!4!!5!/3
SCALE    = 0!!-1!!0
TYPE     = f!!c!!!a
CONTUR   = 3/1
CINT     = !!2/546!2//540!1!1//1
LINE     = !!15/5/3!26/5/3!1//4/2
FINT     = 2;4;6;8;12;16;20;28;40
FLINE    = 0;23-19;17-14
HILO     = !!!!4;2/H#;L#/1016-1100;800-1008/8/2;4
HLSYM    =  
CLRBAR   = 1
WIND     = AM7//2
REFVEC   = 15;.03;.05;.8
TITLE    = 0!1/-2/^ 3-hr Precip., 1000-500-hPa thick., SLP, 10m wind (m/s)!0
TEXT     = .8/2/hw!!!!!2/22/hw
CLEAR    = YES
GAREA    = grid
IJSKIP   =  
PROJ     = lcc/30;-73;60
MAP      = 8//2 + 8//11
MSCALE   = 0
LATLON   =  
DEVICE   = psc
STNPLT   =  
SATFIL   =  
RADFIL   =  
IMCBAR   =  
LUTFIL   =  
STREAM   =  
POSN     = 0
COLORS   = 1
MARKER   = 0
GRDLBL   = 0
FILTER   = no
gdpfun   = sub(ptot^F$1,ptot^F$2)!!sub(hght@500,hght@1000)!!pmsl!kntv(wind)
\$MAPFIL  = hicnus.nws + histus.nws

run


exit

endgdp

gpend

convert -density 100x100 -rotate 90 psc.plt wrf03_$3F$1_pmsl_pre.gif
rm -rf psc.plt




