grads RTMA tmp plot (for gradsusr mailing list)

* this is a sample script
'reinit'
'clear'
'set font 0'
'open rtma.23z.ctl'
'set lat 35 45'
'set lon -85 -70'
'set grid off'
'set mpdraw off'
* data stuff
'tempf = ( (9/5) * (tmp2m - 273) ) + 32'
'dptf = ( (9/5) * (dpt2m - 273) ) + 32'
* 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'
*timeinfo=result
hr=substr(result, 8, 2)
day=substr(result, 11, 10)
'set grads off'
'set display color white'
* THE MAP SHALL BE DISPLAYED!
* 'display smth9(tempf)' ** this went horribly wrong
'd tempf'
'run r/cbaro.gs 2'
****
'set gxout contour'
'set clevs 50 55 60 65 70'
'set ccols 15 12 10 13 3'
'set cthick 10'
'set clab masked'
'set csmooth on'
'd dptf'
* gis
*'set line 99 1 1'
*'draw shp counties'
'set line 1 1 3'
'draw shp states'
'draw title 'hr
*'draw string 1 2'hr
*'draw string 1 5'day
'printim out/RTMA_TMP_'day'_'hr'.png'
'clear'
'quit'