'reinit'
*'open /cpc/analysis/cdas/month/prs/mean/prs.grib.mean.clim.y1981-2010.ctl'
'open /cpc/analysis/cdas/month/prs/mean/prs.grib.mean.y1979-cur.ctl'
'enable print zlat'
'set vpage 0 8.5 0 11'
'set parea 1.0 7.5 6.5 10.0'
'set grads off'
'set looping off' 
'set cterp off'
'set csmooth on'
'set map 1 1 8'
'set string 1 tc 7'
'set strsiz 0.15 0.15'
'draw string 4.25 10.5 FMA 2012'
'draw string 4.25 10.2 Divergent Circulation 50-60W'
'set lat -45 45'
'set lon 295 315'
'define rad = 3.14159/180.'
'define dy = cdiff(lat,y)*rad*6.37e6'
'set time mar2012'
'set z 1 12'
'define u = ave(ugrd,t-1,t+1)'
'define v = ave(vgrd,t-1,t+1)'
'define div = hdivg(u,v)'
'define w = ave(VVEL*864,t-1,t+1)'
'define chi = ave(VPOT,t-1,t+1)'
'set lon 300'
'define avediv = ave(div,lon=300,lon=310)'
'define avew = -ave(w,lon=300,lon=310)'
'define avechi = ave(chi,lon=300,lon=310)'
'define vchi = cdiff(avechi,y)/dy'
'set gxout shaded'
'set zlog on'
'run /cpc/home/wd52vk/grads/util/rgbset.gs'
'set clevs -4 -3 -2 -1 -0.5 0 0.5 1.0 2 3 4'
'set ccols 59 49 47 45 43 41 21 23 25 27 29 79'
'set lat -40 40'
'd smth9(avediv*1000000)'
'set gxout contour'
'set ccolor 1'
'set cthick 2'
'set clevs -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4'
'd smth9(avediv*1000000)'
'set gxout vector'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set cthick 5'
'set ccolor 1'
'set arrlab off'
'd vchi;avew/20'
'run /cpc/home/wd52vk/grads/util/cbarhov.gs'