'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 zlon'
'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.2 FMA 1981-2010'
'draw string 4.25 10.5 Divergent Circulation EQ-10S'
'set dfile 1'
'set lat -15 5'
'set lon 180 365'
'define rad = 3.14159/180.'
'define dx = cdiff(lon,x)*rad*6.37e6'
'set t 1 12'
'set z 1 12'
'define uc = ugrd'
'define vc = vgrd'
'define divc = hdivg(uc,vc)'
'define wc = VVEL*864'
'define chic = VPOT'
'modify uc seasonal'
'modify vc seasonal'
'modify divc seasonal'
'modify wc seasonal'
'modify chic seasonal'
'set lat -5'
'define avedivc = ave(divc,lat=-10,lat=0)'
'define avewc = -ave(wc,lat=-10,lat=0)'
'define avechic = ave(chic,lat=-10,lat=0)'
'define uchic = cdiff(avechic,x)/dx'
'modify avedivc seasonal'
'modify avewc seasonal'
'modify avechic seasonal'
'modify uchic seasonal'
'set t 3'
'define adivc = ave(avedivc,t-1,t+1)'
'define auchic = ave(uchic,t-1,t+1)'
'define awc = ave(avewc,t-1,t+1)'
'set dfile 2'
'set time feb2012 apr2012'
'set lat -15 5'
'set lon 180 365'
'define rad = 3.14159/180.'
'define dx = cdiff(lon,x)*rad*6.37e6'
'define u = ugrd'
'define v = vgrd'
'define div = hdivg(u,v)'
'define w = VVEL*864'
'define chi = VPOT'
'set lat -5'
'define avediv = ave(div,lat=-10,lat=0)'
'define avew = -ave(w,lat=-10,lat=0)'
'define avechi = ave(chi,lat=-10,lat=0)'
'define uchi = cdiff(avechi,x)/dx'
'set time mar2012'
'define adiv = ave(avediv,t-1,t+1)'
'define auchi = ave(uchi,t-1,t+1)'
'define aw = ave(avew,t-1,t+1)'
'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 lon 200 360'
'd smth9(adivc*1000000)'
'set gxout contour'
'set ccolor 1'
'set cthick 2'
'set clevs -7 -6 -5 -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4 5 6 7'
'd smth9(adivc*1000000)'
'set gxout vector'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set cthick 5'
'set ccolor 1'
'set arrlab off'
'd auchic;awc/20'
'run /cpc/home/wd52vk/grads/util/cbarhov.gs'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 1.0 7.5 2.0 5.5'
'set grads off'
'set map 1 1 7'
'set string 1 tc 7'
'set strsiz 0.15 0.15'
'draw string 4.25 5.7 FMA 2012'
'set zlog on'
'set gxout shaded'
'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 lon 200 360'
'd smth9(adiv*1000000)'
'set gxout contour'
'set ccolor 1'
'set cthick 2'
'set clevs -7 -6 -5 -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4 5 6 7'
'd smth9(adiv*1000000)'
'set gxout vector'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set cthick 5'
'set ccolor 1'
'set arrlab off'
'd auchi;aw/20'
'run /cpc/home/wd52vk/grads/util/cbarhov.gs'
'printim zlon-div-circ-0-10s-fma2012.gif white'
'clear'
'set vpage 0.0 8.5 0.0 11.0'
'set parea 1.0 7.5 6.5 10.0'
'set grads off'
'set map 1 1 7'
'set string 1 tc 7'
'set strsiz 0.15 0.15'
'draw string 4.25 10.5 Anomalous Divergent Circulation EQ-10S'
'draw string 4.25 10.2 FMA 2012'
'set zlog on'
'set gxout shaded'
'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 lon 200 360'
'd smth9((adiv-adivc)*1000000)'
'set gxout contour'
'set ccolor 1'
'set cthick 2'
'set clevs -7 -6 -5 -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4 5 6 7'
'd smth9((adiv-adivc)*1000000)'
'set gxout vector'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set cthick 5'
'set ccolor 1'
'set arrlab off'
'd auchi-auchic;(aw-awc)/10'
'run /cpc/home/wd52vk/grads/util/cbarhov.gs'
'printim zlon-div-circ-anoms-0-10s-fma2012.gif white'