* Script to print 2 panel plot of GPCP top 5 dry (wet) July years minus climatology.
* Topics include: a) composites and b) custom color scale
* top 5 dry years are: 86,87,89,90,02
* top 5 wet years are: 96,98,00,01,84; EXCLUDING 2005

* clear GrADS settings and open data file
'reinit'
'sdfopen E:\academic\digital_FILES\NetCDF\precip.mon.mean.nc'

* define geographic coordinates 
'set lev 0'
'set lat 15 21'
'set lon 275 290'

* set page dimensions for first plot
'set vpage 0 5.5 4.5 8'
'set grads off'
'set gxout shaded'

* define custom colors (shades of brown)
'set rgb 16 153 52 4'
'set rgb 17 217 95 14'
'set rgb 18 254 153 41'
'set rgb 19 254 217 142'
'set rgb 20 255 255 212'
'set rgb 21 255 255 204'
'set rgb 22 194 230 153'
'set rgb 23 120 198 121'
'set rgb 24 49 163 84'
'set rgb 25 0 104 55'
'set rgb 26 255 255 255'

* define the contour intervals and assign colors to each interval
'set clevs -1.25 -1 -.75 -.5 -.25 0 .25 .5 .75 1 1.25'
'set rbcols 16   17     18   19     20    26    26    21    22    23    24   25'

* create a 'dry' composite and the climatology 
'set csmooth on'
'define gpcpdry =(precip(t=91)+precip(t=103)+precip(t=127)+precip(t=139)+precip(t=283))/5'
'define gpcp = ave (precip, t=7, t=343, 12)'

* display varible, title, and color bar
'd gpcpdry-gpcp'
'draw title July GPCP dry composite minus climatology'
'run cbarn'

* apply similar steps above  to a different composite with a different color scale (shades of green)
'set lev 0'
'set lat 15 21'
'set lon 275 290'
'set vpage 5.5 11 4.5 8'
'set grads off'
'define gpcpwet =(precip(t=211)+precip(t=235)+precip(t=259)+precip(t=271)+precip(t=67))/5'
'define gpcp =ave (precip, t=7, t=343, 12)'
'set clevs -1.25 -1 -.75 -.5 -.25 0 .25 .5 .75 1 1.25'
'set rbcols 16   17     18   19     20    26    26    21    22    23    24   25'
'd gpcpwet-gpcp'
'draw title July GPCP wet composite minus climatology'
'run cbarn'

'set vpage off'
'printim c:\testdata\msd_images\GPCP_composite_thesis.jpg gif x800 x600 white'
'disable print'