Grads 等値線をかいてみるスクリプトファイル

*-----test1.gs-----

reinit

'open mode2.ctl'
'set lat 39.0 40.0'
'set lon 140.0  142.0'
'set lev 5'
'set map 1 1 6'

'q gxinfo'

'q file'
'q dims'

'set t 1'
'set grads off'

'set gxout shaded'

'set clevs -1 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4  0.6  0.8  1.0'

'd mode'

'cbarn'

'set gxout contour'

*-----test2.gs-----

reinit

'open s1.ctl'
'set lat 35.0 40.0'
'set lon 140.0  145.0'
'set lev 5'
'set map 1 1 6'

'q gxinfo'
'q file'
'q dims'

t=1
while (t<=745)
'c'
'set t 't''
'set grads off'

'set gxout shaded'

'set clevs 30 30.25 30.5 30.75 31.0 31.25 31.5 31.75 32.0 32.25 32.5 32.75 33.0 33.25 33.5 33.75 34.0 34.25 34.5 34.75 35.0'

'd mode'

'set strsiz 0.2 0.2'

'draw string 2 6 t='t
*pull input
t=t+2

endwhile

'cbarn'

'set gxout contour'

'd mode'