GrADS学习资料：实习八描述语言的使用（2小时）.doc

实习八 描述语言的使用（2小时）

1．实习资料：现在“data”文件夹下有全球海表温度资料sst.mnmean.nc。

2．实习要求：
根据所提供的数据资料，计算夏季海温距平场，并循环绘出1971年至2000年夏季海温距平图。要求海温正距平区填色，画色标，并将年份作为标题。最终将图形保存。所有命令编写于.gs文件中。

3. 实习目的：
掌握循环的使用，进一步练习函数的使用。

4. 实习步骤：
4.1 在写字板或记事本程序中按要求编写*.gs文件；
'reinit'
'sdfopen d:\data\sst.mnmean.nc'
'set fwrite d:\data\sum.grd'
'set gxout fwrite'
i=1405
while(i<=1764)
'set x 1 180'
'set y 1 89'
'set t 1'
'set t 'i''
'define sstsum=ave(sst,t='i+5',t='i+7')'
'd sstsum'
i=i+12
endwhile
'disable fwrite'
;

dset d:\data\sum.grd
undef 32767.00
xdef 180 linear 0 2
ydef 89 linear -88 2
zdef 1 linear 1 1
tdef 30 linear jun1971 1yr
vars 1
sstsum 0 -999 Monthly Means of Sea Surface Temperature
endvars

'reinit'
'open d:\data\sum.ctl'
'enable print d:\data\shixi8.gmf'
j=1971
i=1
while(i<=30)
'define sumave=ave(sstsum,t=1,t=30,1)'
'define anom=sstsum-sumave'
'set gxout shaded'
'set t 'i''
'set cmin 0'
'd anom'
'cbarn 1 0 5 0.35'
'draw title summeraverage'j''
'set gxout contour'
'd anom'
'print'
'c'
i=i+1
j=j+1
endwhile
'disable print'
;

4.2启动GrADS绘图软件，运行*.gs文件