整层水汽通量的问题

jeffggg

求整层水汽通量的时候 grads提示我说图片那样   我很奇怪 pres是地面气压 我只是求了平均而已啊 哪里出错了呢？
'reinit'
'open D:\Pictures\quv.ctl'
'open D:\scilencedata\PRESpentad\aveprespen\PRES4557.ctl'
'set lat 0 50'
'set lon 50 160'
'set xlint 20'
'set ylint 10'
'set arrscl 0.5 400'
'set grads off'
'set grid off'
'set t 45'
'define P=ave(pres.2,t=45,t=57)'
'set z 1'
'define quall=vint(p,qu,300)*9.8'
'define qvall=vint(p,qv,300)*9.8'
'd quall;qvall'
*'d qvall'
'd hdivg(quall,qvall)'
;

蔚蓝的天

'open D:\scilencedata\PRESpentad\aveprespen\PRES4557.ctl'
提示这个有问题啊

叶の幽灵

建议直接写四则运算公式，p=(pres.2(t=45)+pres.2(t=...)+...)/...
还有 z=1 要放在计算的前面

金玉姻缘

你好，请问quall=vint(p,qu,300)*9.8，这个qu的维数是怎么样的，我查看了vint这个函数，他写的是被积分量应该是三维的，（x,y,p），怎么构造这个呢，这一块不懂，请指教

cgccq

没搞清楚?太奇怪了

扫眉才子

学习了。。。一看知道我是回复骗积分的{:eb334:}{:eb334:}{:eb334:}

晨心快乐

蔚蓝的天 发表于 2013-5-2 09:50
'open D:\scilencedata\PRESpentad\aveprespen\PRES4557.ctl'
提示这个有问题啊

不知道楼主自己调试过没有啊 咱们论坛好像有过一个的，要是有一些注意点提示就好啦

phys

想问一下楼主，地面气压为啥要求平均啊？

x6xin

帮顶了!!!!!