'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 nps'
'set frame off'
'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 10.5 500-hPa Height (dam) & Abs Vort Adv (shaded)'
'draw string 4.25 10.2 00 UTC 14 March 1993' 
'set time 00Z14mar1993'
'define rad = 3.14159/180'
'define g = 9.80665'
'set lat 15 70'
'set lon -150 -40'
'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'
'run /export/lnx286/wd52vk/grads/util/rgbset.gs'
'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 20 65'
'set lon -140 -50'
'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 500z-advq-00Z14mar1993.gif white'