grads练习题参考

一、
data.ctl内容:
DSET ^data.grd
UNDEF 9.999E+20
TITLE 10 year data
XDEF 37 LINEAR 60.0 2.5
YDEF 17 LINEAR 0.0 2.5
TDEF 120 LINEAR jan1980 1mo
ZDEF 5 LEVELS 850 700 500 300 200
VARS 4
u 5 99 * U Winds
v 5 99 * V Winds
z 5 99 * Geopotential Heights
sst 0 99 * Sea Surface Temperature
ENDVARS

三、
data3.gs内容:
'reinit'
'open data.ctl'
'set lon 60 150'
'set lat 0 40'
'set z 3'
'set time jan1983'
'set gxout contour'
'd z'
'draw string 2.5 8.0 jan1983 500hpa Geopotential height'
'set gxout fwrite'
'set fwrite d:\llp\height.grd'
'd z'
'disable fwrite'
;

四、
data4.gs内容:
'reinit'
'open data.ctl'
'set lon 60 150'
'set lat 0 40'
'set z 2'
'enable print d:\llp\stream.gmf'
i=1
while (i<=120)
'set t ' i
'clear'
'q time'
res=subwrd(result,3)
'set gxout stream'
'd u;v'
'draw string 2.5 8.0 ' res ' 700hpa stream line '
'print'
i=i+1
endwhile
'disable print'
;

五、
data5.gs内容:
'reinit'
'open data.ctl'
'set lon 60 150'
'set lat 0 40'
'set z 4'
'set time jul1980'
'set gxout contour'
'd z'
'draw string 2.5 8.0 jul1980 300hpa Geopotential height'
'q w2xy 85 30'
x1=subwrd(result,3)
y1=subwrd(result,6)
'set string 4 c 8 0'
'draw string 'x1' 'y1' D'
;