- 积分
- 667
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-5-21
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
相求逐日区域平均的散度场,但是grads一直提示mdiv为无效变量,想问下大神错在哪里,在这个gs之前已经将多个nc逐日资料合并为一个dat,以下为其中u风场的ctl和gs,比较急,确实找了很久也没有找出原因。
dset D:\studing\dipin\1979-2015\jifeng\2\u1.grd
undef -9.99E+33
title daily.OLR.b19790101.e20151231
xdef 1 linear 82.5 2.5
ydef 1 linear 27.5 2.5
zdef 1 levels 600
tdef 13514 linear 01jan1979 1dy
vars 1
u 0 99 afd
endvars
'reinit'
'open D:\studing\dipin\1979-2015\jifeng\2\u1.ctl'
'open D:\studing\dipin\1979-2015\jifeng\2\v1.ctl'
'set gxout fwrite'
'set fwrite D:\studing\dipin\1979-2015\jifeng\2\mdiv.grd'
'set lev 600'
'set t 1 13514'
'define div=hdivg(u.1,v.2)'
'define mdiv=aave(div,lon=82.5,lon=100,lat=27.5,lat=37.5)'
id=1
while(id<=13514)
'set t 'id''
'd mdiv'
id=id+1
endwhile
'disable fwrite'
'print'
'reinit'
;
|
|