- 积分
- 463
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-30
- 最后登录
- 1970-1-1
|
发表于 2012-10-5 19:15:37
|
显示全部楼层
楼主您好,我之前贴上来的图应该有些问题,数据的描述文件应该没有写对
*flying_nuist
'reinit'
'open e:\...\quv.ctl'
'open e:\...\pres_mean.ctl'
'set fwrite e:\...\quv_all.grd'
'set gxout fwrite'
'set z 1'
'define quall=vint(pres.2,qu,300)*9.8'
'define qvall=vint(pres.2,qv,300)*9.8'
'd quall'
'd qvall'
'd hdivg(quall,qvall)'
'disable fwrite'
'reinit'
;
在这个脚本执行后输出的quv_all.grd中,quall,qvall的垂直方向只有一层,对应的ctl文件中 zdef应该怎么写?
dset D:\9801\mean\quvall1.grd
title mean quvall data(1998)
undef -9.96921e+36
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 levels ???
tdef 1 linear 00Z01JAN1998 1440mn
vars 2
quall 0 -999 mean qu
qvall 0 -999 mean qv
endvars
还请楼主指点
|
|