'open /export-5/sgi9/data/wd51jh/stnusmap/grid.ctl'
'open w.diffjul.ctl'
'open /export-5/sgi9/data/wd51jh/stnusmap/maskus.ctl'
'enable print meta'
 
'set display color white'
'clear'
'exec /export-5/sgi9/data/wd51jh/stnusmap/rgbset'
 
'define crs=maskout(oacres(x,f.2),mask.3-1)'

* top panel
'set grads off'
'set grid off'
'set mpdset mres'
'set poli on'
'set lon -129 -65'
'set lat 24 51'
'set ylint 5'
'set xlint 20'
'set gxout shaded'
*'set clevs  100 200 300 400 500 600'
*'set ccols 0  31  32  33  34  35  36'
'set clevs  -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70'
'set ccols 27  26   25   24  23   22  21  0  31 32 33 34 35  36    37'
'set gxout shaded'
'set csmooth on'
'd crs'
'set gxout contour'
*'set clevs  100 200 300 400 500 600'
'set clevs  -70 -60 -50 -40 -30 -20 -10 10 20 30 40 50 60 70'
'd crs'
'run /export-5/sgi9/data/wd51jh/stnusmap/cbar.gs'
'draw title Calculated Soil Moisture Difference (9100-6170) July'
*
'print'
'clear'
'reinit'
*
'quit'

function gettmon(mon)
 if (mon=1); titlemon='JAN'; endif;
 if (mon=2); titlemon='FEB'; endif;
 if (mon=3); titlemon='MAR'; endif;
 if (mon=4); titlemon='APR'; endif;
 if (mon=5); titlemon='MAY'; endif;
 if (mon=6); titlemon='JUN'; endif;
 if (mon=7); titlemon='JUL'; endif;
 if (mon=8); titlemon='AUG'; endif;
 if (mon=9); titlemon='SEP'; endif;
 if (mon=10); titlemon='OCT'; endif;
 if (mon=11); titlemon='NOV'; endif;
 if (mon=12); titlemon='DEC'; endif;
return titlemon