发布时间 : 星期六 文章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'