'set display color white'
'set grads off'

'open zgrads.ctl'

'set parea 1.0 7.0 1.0 10.3592'
'set xlint 10'
'set ylint 10'
'set xlopts 20 3 0.04'
'set ylopts 20 3 0.04'
'set rgb 20 168 48 16'
'set mproj scaled'
'set mpdset mres'
'set poli on'
*'set frame off'
'set grid on 1 20 3'

'set gxout model'
'set digsiz 0.05'
'set mdlopts dig3'
'set stid on'
'set lon 118 127'
'set lat 5 21'
'set rgb 20 168 48 16'
'set map 20'
'set ccolor 4'
*'d us;vs;ts;ds;ps;pdel;cld;wx;vis'
'd us;vs;ts;ds;0.0;0.0;cld;wx'

'q time'
t = subwrd(result, 3)
'set strsiz 0.10'
'set string 1 bl 12'
'draw string 1.01 10.2 SURFACE PLOT + GSM(ft6)'
'draw string 1.01 10.05 ' t

'close 1'
'open zgrib2.ctl'
'set lon 118 127'
'set lat 5 21'
'set gxout grid'
'set digsize 0.07'
'set dignum 1'
'set ccolor 8'
'd TMP2m - 273.15'

'printim zgradsph.png x1700 y2200'
'quit'