- 积分
- 6139
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-11-6
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2016-5-17 22:20:53
|
显示全部楼层
不好意思,是我没有说清楚
我刚刚改了一下程序,出来的数值应该对了,您看一下应该就会明白我求的什么。
几个特定年份的距平合成,是我把modify seasonal的用法搞错了。
'reinit'
'sdfopen d:\zhangjie\hgt.mon.mean.nc'
'enable print d:\biyelunwen\otherout\EAST_SUM_500HGT.22.gmf white'
'set lon 60 180'
'set lat -10 50'
'set z 6'
'a6=ave(hgt,time=1jun1971,time=1dec2010,12)'
'a7=ave(hgt,time=1jul1971,time=1dec2010,12)'
'a8=ave(hgt,time=1aug1971,time=1dec2010,12)'
'define m=(a6+a7+a8)/3'
'set grads off'
'set grid off'
'set xlopts 1 2 0.1'
'set ylopts 1 2 0.1'
'set xlint 30'
'set ylint 5'
'set gxout contour'
*'set cint 0.5'
'define a=ave(hgt-m,time=jun1952,time=aug1952)'
'define b=ave(hgt-m,time=jun1970,time=aug1970)'
'define c=ave(hgt-m,time=jun1973,time=aug1973)'
'define d=ave(hgt-m,time=jun1977,time=aug1977)'
'define e=ave(hgt-m,time=jun1983,time=aug1983)'
'define f=ave(hgt-m,time=jun1998,time=aug1998)'
'define nasst=(a+b+c+d+e+f)/6'
'd nasst'
'draw title EAST HGT'
'print'
'disable print'
;
如果有错希望继续指正,感谢!
|
|