GrADS

nc文件的维数设置、格距、缺测值的查询：q ctlinfo

2）所编写的*.gs文件； 'reinit' 'reinit'

'sdfopen F:\\1\\sst.mnmean.nc' 'set t 1 12'

'define avesst=ave(sst,t+1404+0,t+1764,12)' 'modify avesst seasonal' 'set t 1405 1764'

'define anom=sst-avesst' 'set t 1729'

'set vpage 0 5.5 2 8' 'd anom'

'set gxout shaded' 'set cmin 0' 'd anom' 'set clab on'

'set gxout contour' 'set cmax 0' 'd anom'

'cbarn 0.8 0 4 0.5' 'draw title JAN 1998' 'set t 1741'

'set vpage 5.5 11 2 8' 'd anom'

'set gxout shaded' 'set cmin 0' 'd anom' 'set clab on'

'set gxout contour' 'set cmax 0' 'd anom'

'cbarn 0.8 0 4 0.5' 'draw title JAN 1999'

'printim F:\\1\\shuang7.png white' 3）所绘图形。

实习八 描述语言的使用

1. 实习目的：

掌握循环的使用，进一步练习函数的使用。 2. 实习要求：

（1）计算夏季海温距平场；

（2）循环绘出1971年至2000年夏季海温距平图。要求海温正距平区填色，画色标，并将年份作为标题。最终将图形保存。 （3）所有命令编写于.gs文件中。 3. 实习步骤：

(1)在写字板或记事本程序中按要求编写*.gs文件； (2)启动GrADS绘图软件，运行*.gs文件 4. 实习报告： 1）所用资料

“data”文件夹下有全球海表温度资料sst.mnmean.nc。

nc文件打开：sdfopen 路径\\ sst.mnmean.nc

nc文件的维数设置、格距、缺测值的查询：q ctlinfo

2）所编写的*.gs文件 'reinit'

'sdfopen F:\\1\\sst.mnmean.nc' 'set gxout fwrite'

'set fwrite e:\\dwj\\8.6.grd' it=1410

while(it<=1758) 'set lev 0'

'set lon 90 200' 'set lat -20 60' 'set t 'it'' 'd sst' it=it+12

endwhile 'disable fwrite' ;

*.ctl文件

dset F:\\1\\8.6.grd

title NOAA Extended Reconstructsd SST V3 undef 32767

xdef 56 linear 90 2 ydef 41 linear -20 2 zdef 1 linear 1 1

tdef 30 linear 00Z01Jun1971 1yr vars 1

sst 0 -999 Monthly Means of Sea Surface Temperature ENDVARs 7月：

*.gs文件： 'reinit'

'sdfopen F:\\1\\ sst.mnmean.nc' 'set gxout fwrite'

'set fwrite e:\\dwj\\ 8.7.grd' it=1411

while(it<=1759) 'set lev 0'

'set lon 90 200' 'set lat -20 60' 'set t 'it'' 'd sst' it=it+12

endwhile 'disable fwrite' ;

*.ctl文件

dset F:\\1\\ 8.7.grd

title NOAA Extended Reconstructsd SST V3 undef 32767

xdef 56 linear 90 2 ydef 41 linear -20 2

zdef 1 linear 1 1

tdef 30 linear 00Z01Jun1971 1yr vars 1

sst 0 -999 Monthly Means of Sea Surface Temperature ENDVARs 8月：

*.gs文件： 'reinit'

'sdfopen F:\\1 \\sst.mnmean.nc' 'set gxout fwrite'

'set fwrite F:\\1\\8.8.grd' it=1412

while(it<=1760) 'set lev 0'

'set lon 90 200' 'set lat -20 60' 'set t 'it'' 'd sst' it=it+12

endwhile 'disable fwrite' ;

*.ctl文件

dset F:\\1\\8.8.grd

title NOAA Extended Reconstructsd SST V3 undef 32767

xdef 56 linear 90 2 ydef 41 linear -20 2 zdef 1 linear 1 1

tdef 30 linear 00Z01Jun1971 1yr vars 1

sst 0 -999 Monthly Means of Sea Surface Temperature ENDVARs

1971-2000年夏季（678月）： *.gs文件： 'reinit'

'open F:\\1\\8.6.ctl' 'open F:\\1\\8.7.ctl' 'open F:\\1\\8.8.ctl' 'set gxout fwrite'

'set fwrite F:\\1\\ 8.678.grd' it=1

while(it<=30) 'set lon 90 200'