'open /cpc/r2/6hr/pgb/pgb.ctl'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 1.0 7.5 6.0 10.0'
'set annot 1 5'
'set map 1 1 1'
'set mproj sps'
'set frame on'
'set cterp off'
'set csmooth on'
'set lat -65 -20'
'set lon 270 340'
'set lev 1000'
'set grads off'
'set string 1 tc 7'
'set strsiz 0.12 0.12'
'draw string 4.25 10.2 SLP (hPa-1000)i & 1000-hPa Wind' 
'draw string 4.25 10.5 12 UTC 27 June 1994' 
'set time 12Z27jun1994'
'run /export/lnx286/wd52vk/grads/util/rgbset.gs'
'set xlint 20'
'set ylint 10'
'set gxout contour'
'set cint 4'
'set ccolor 1'
'set cthick 1'
'set cstyle 1'
'd presmsl/100.-1000.'
'set gxout vector'
'set arrscl 0.3 15'
'set arrowhead -0.3'
'set arrlab off'
'set ccolor 1'
'set cthick 3'
'd ugrdprs;vgrdprs'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 1.0 7.5 1.5 5.5'
'set annot 1 5'
'set map 1 1 1'
'set mproj sps'
'set frame on'
'set cterp off'
'set csmooth on'
'set lev 500'
'set grads off'
'set string 1 tc 7'
'set strsiz 0.12 0.12'
'draw string 4.25 5.7 500-hPa Height (dam) & Abs Vort Adv (shaded)'
'define rad = 3.14159/180'
'define g = 9.80665'
'set lat -70 -15'
'set lon 260 350'
'define cor = 2*7.292e-5*sin(lat*rad)'
'define dx = cdiff(lon,x)*rad*6.37e6*cos(lat*rad)'
'define dy = cdiff(lat,y)*rad*6.37e6'
'define q = hcurl(ugrdprs,vgrdprs)+cor'
'define advq = -ugrdprs*cdiff(q,x)/dx - vgrdprs*cdiff(q,y)/dy'
'set xlint 20'
'set ylint 10'
'set gxout shaded'
'set clevs -5 -4 -3 -2 -1 0 1 2 3 4 5'
'set ccols 59 48 47 45 43 41 21 23 25 27 29 79'
'set lat -65 -20'
'set lon 270 340'
'd advq*1e9'
'set gxout contour'
'set ccolor 1'
'set cstyle 3'
'set cthick 1'
'set cint 1'
*'d q*100000'
'set gxout contour'
'set cint 6'
'set ccolor 1'
'set cthick 1'
'set cstyle 1'
'set clab off'
'd hgtprs/10.'
'set gxout contour'
'set clevs 504 522 540 558 576'
'set ccolor 1'
'set cthick 1'
'set cstyle 1'
'set clab on'
'd hgtprs/10.'
'run /export/lnx286/wd52vk/grads/util/cbarhovn.gs'
'printim slp-uv1000-500z-advq-12Z27jun1994.gif white'