'reinit'
'open /cpc/home/wd52vk/samer/satellite/irwin-daily-clim365-boxes-y1980-2009.ctl'
'open /cpc/home/wd52vk/samer/satellite/irwin-SA-daily-anoms-var-boxes-jan1980-2009.ctl'
'set vpage 0 8.5 0 11.0'
'set parea 0.5 4.0 7.0 10.0'
'set dfile 1'
'set grads off'
'set mpdset hires'
'set map 1 1 3'
'set annot 1 5'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.15 0.15'
'draw string 4.25 10.5 January 1980-2009 Brightness Temperature'
'set strsiz 0.15 0.15'
'draw string 2.25 10.2 Mean'
'set lat -54 14'
'set lon -89 -21'
'set time 12z16jan1980'
'define irave=ave(irwin,t-15,t+15)'
'set z 1'
'run /cpc/home/wd52vk/grads/util/rgbset.gs'
'set clevs 260 263 266 269 272 275 278 281 284 287 290'
*'set ccols 59 49 39 37 35 33 31 71 73 75 77 79' 
'set ccols 59 49 47 45 43 41 21 22 23 24 25 26' 
'set gxout shaded'
'set ylint 10'
'set xlint 20'
'set ylab on'
'd irave'
'run /cpc/home/wd52vk/grads/util/cbarhovn.gs'
'set vpage 0 8.5 0 11.0'
'set parea 4.5 8.0 7.0 10.0'
'set dfile 1'
'set grads off'
'set mpdset hires'
'set map 1 1 3'
'set annot 1 5'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.15 0.15'
'draw string 6.25 10.2 Variance'
'set dfile 2'
'set lat -54 14'
'set lon -89 -21'
'set t 1'
'set clevs 30 60 90 120 150 180 210 240 270 300'
'set ccols 0 41 43 45 47 49 21 23 25 27 29' 
'set gxout shaded'
'set ylint 10'
'set xlint 20'
'set ylab on'
'd irvar'
'run /cpc/home/wd52vk/grads/util/cbarhovn.gs'
'gxprint irwin-mean-variance-jan1980-2009.png white'