GrADS RTMA Temp plot (2011-09-13)

****************************
* RTMA Plotting Script (tmp2m only for now)
*
* 2011 09 05  Script transferred to wxjoe.com
*	      Switched data source to openDAP server
* 2011 09 12  Discovered set stat on / gxout stat
*
***************************
'reinit'
'clear'
*'!date -u +%Y%m%d%t%H > date.tmp'
**** above is being created by a shell script
**** returns YYYYMMDD	HH
tmpdate=read('date.tmp')
dapdate=subwrd(tmpdate,2)
currhr=subwrd(tmpdate,3)
'set font 0'
'sdfopen http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p5'dapdate'/rtma2p5_anl_'currhr'z'
'set lat 25 55'
'set lon -125 -70'
'set grid off'
'set mpdraw off'
* data stuff
'tempf = ( (9/5) * (tmp2m - 273) ) + 32'
'dptf = ( (9/5) * (dpt2m - 273) ) + 32'
*'define maxt=max(tempf)'
* graphics stuff
'set gxout shaded'
*'set cint 4' - I think this doesn't matter
* call color bar stuff
'run r/colorbar.gs'
'set clevs 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90'
'set ccols 15 16 17 18 19 36 38 40 41 43 44 45 46 48 50 51 53 55 57 58 60 62 63 64 65 0'
* display stuff
'q time'
hr=substr(result, 8, 3)
day=substr(result, 11, 9)
'set grads off'
'set display color white'
* THE MAP SHALL BE DISPLAYED!
'd tempf'
'run r/cbaro.gs 2'
*'set gxout contour'
*'set clevs 60 70 80 90'
*'set clab off'
*'d tempf'
****
*'set gxout contour'
*'set clevs 60 70'
*'set ccols 10 13'
*'set cthick 10'
*'set clab masked'
*'set csmooth on'
*'d smth9(smth9(dptf))'
* gis
'set line 99 1 1'
'draw shp gis/counties'
'set line 1 1 6'
'draw shp gis/states'
'draw title RTMA Temperature (F)'
*'draw string 1 2'hr
*'draw string 1 5'day
'gxyat img/rtma_temp_latest.png'
'clear'
'quit'