- 积分
- 574
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-9-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 呆萌萱 于 2016-4-7 21:25 编辑
我先用月平均资料计算生成年平均资料,脚本如下:
'reinit'
'sdfopen F:\benke\xiaobo\1979_2015_3.nc'
'set fwrite F:\benke\xiaobo\100ave.grd'
'set gxout fwrite'
'set lev 100'
it=13
'set t 'it''
while(it<=444)
'define tave=ave(t,t='it',t='it+11')'
'd tave'
it=it+12
endwhile
'disable fwrite'
'reinit'
;
之后我编写了100ave.grd的描述文件,脚本如下:
dset F:\benke\xiaobo\100ave.grd
title ave of year
undef -9.99e+33
*dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 levels 100
tdef 36 linear JAN1980 1yr
vars 1
tave 36 0 Temperature
endvars
我再用gs算区域平均,脚本如下:
'reinit'
'open F:\benke\xiaobo\100ave.ctl'
'set fwrite F:\benke\xiaobo\100quyu.grd'
'set gxout fwrite'
it=1
while(it<=36)
'set t 'it''
'define a=aave(tave,lon=75,lon=105,lat=25,lat=40)'
'd tloop(a)'
it=it+1
endwhile
'disable fwrite'
'print'
'reinit'
;
但是最后出现了如图的错误,这怎么解决啊???
|
|