grads计算整层水汽通量一直出错，求助大神解决啊

tmw

第一个ctl文件：dset d:/tong/tongji/sqtl/tqzcsqtl.dat
undef -9.99e+8
xdef 72 linear 0 2.5
ydef 36 linear 0 2.5
zdef 6 levels 1000 925 850 700 600 500
tdef 1 linear dec2007 1mo
vars 3
qu 6 99 eigenvactors matrix of EOF
qv 6 99 eigenvactors matrix of EOF
quv 6 99 eigenvactors matrix of EOF
endvars
第二个ctl文件dset d:/tong/tongji/sqtl/pres.12.mean.dat
undef -9.99e+8
xdef 72 linear 0 2.5
ydef 36 linear 0 2.5
zdef 1 linear 1 1
tdef 1 linear dec2007 1mo
vars 1
pres 0 99 eigenvactors matrix of EOF
endvars
整层水汽通量积分gs文件：'reinit'
'open d:/tong/tongji/sqtl/tqzcsqtl.ctl'
'open d:/tong/tongji/sqtl/pres.12.mean.ctl'
'set fwrite d:/tong/tongji/sqtl/all_quv.dat'
'set gxout fwrite'
'set lat 0 90'
'set lon 0 180'
'set t 1'

'set z 1'
'define qua=vint(pres.2,qu,500)*9.8'
'define qva=vint(pres.2,qv,500)*9.8'
'define quva=hdivg(qua,qva)'
'd qua'
'd qva'
'd quva'

'disable fwrite'
'c'
'reinit'
grads提示：error from vint :too few levels in default file...

daodang2008

看不懂，解决了没，我也要程序

雨恬魅影

同求。。。

丫丫君

你的第一个ctl有误。

贝得

xdef 72 linear 0 2.5
ydef 36 linear 0 2.5

你这样设置只能用到87.5°N和177.5°E，后面就不能那样写了
'set lat 0 90'
'set lon 0 180'

'define qva=vint(pres.2,qv,500)*9.8'
你这里为什么还要×9.8？计算整层积分的时候要自动除的，再乘回来是因为你算分层除过么？

mengmeng12345

我也想知道，请问楼主解决了没?