整层积分的水汽通量散度

奚子惠

求助各位大神，有没有遇到这样的情况，以下的gs是qu和qv的整层积分，然后再算一个整层水汽通量的散度。
'reinit'
'open F:\qu.ctl'
'open F:\qv.ctl'
'open F:\pres.ctl'

'set fwrite F:\flux_all_u.grd'
'set gxout fwrite'
ta=1
while(ta<=35)
'set t 'ta''
'set z 1'
'set y 1 37'
'set x 1 144'
'define qua=vint(pres.3,qu,300)*9.8'
'define qva=vint(pres.3,qv.2,300)*9.8'
'd qua'
ta=ta+1
endwhile
'disable fwrite'

'set fwrite F:\flux_all_v.grd'
'set gxout fwrite'
ta=1
while(ta<=35)
'set t 'ta''
'set z 1'
'set y 1 37'
'set x 1 144'
'define qua=vint(pres.3,qu,300)*9.8'
'define qva=vint(pres.3,qv.2,300)*9.8'
'd qva'
ta=ta+1
endwhile
'disable fwrite'

'set fwrite F:\flux_all_div.grd'
'set gxout fwrite'
ta=1
while(ta<=35)
'set t 'ta''
'set z 1'
'set y 1 37'
'set x 1 144'
'define qu=vint(pres.3,qu,300)*9.8'
'define qv=vint(pres.3,qv.2,300)*9.8'
'd hdivg(qu,qv)'
ta=ta+1
endwhile
'disable fwrite'

'reinit'
;


问题是：前两部分算qu和qv可以算出来 ，但是第三部分只能算出部分，我输出看了一下，开始值还算正常，到后面大的很离谱，再后面全部变成inf。
于是 我将第三部分算散度的单独拿出来，不做循环就算第34年的散度，是可以出图的，值大概在-0.2~0.15，但是我做了35年的循环后，在grads里set t 34 ，提示：constant value。 value=-inf。请问这是怎么回事？这个循环有问题嘛？但为啥qu和qv的循环就没问题呢

奚子惠

哎 不好意思让大家见笑了 问题找到了 很弱智 真想把这个帖子删了。。囧 当局者迷老半天看不出来

sjysjxyf66

奚子惠 发表于 2015-1-21 23:22
哎 不好意思让大家见笑了 问题找到了 很弱智 真想把这个帖子删了。。囧 当局者迷老半天看不出来

把您的问题告诉大家把，请留着供大家学习，起码我要学习一下啊，谢谢

奚子惠

sjysjxyf66 发表于 2015-1-22 09:44
把您的问题告诉大家把，请留着供大家学习，起码我要学习一下啊，谢谢

'define qu=vint(pres.3,qu,300)*9.8'
'define qv=vint(pres.3,qv.2,300)*9.8'这两句要改成：

'define qua=vint(pres.3,qu,300)*9.8'
'define qva=vint(pres.3,qv.2,300)*9.8'，等号左边要换个变量名，跟右边不能一样哦

sjysjxyf66

奚子惠 发表于 2015-1-22 10:32
'define qu=vint(pres.3,qu,300)*9.8'
'define qv=vint(pres.3,qv.2,300)*9.8'这两句要改成：

在所难免，有时就是一叶障目啊

slc20050110202

感谢，多谢楼主

slc20050110202

感谢，多谢楼主

xiaohezhou

好，学习了。

金玉姻缘

楼主你好，我要计算地面到500hpa整层水汽通量散度，都需要哪些变量啊，公式在哪里呢，我查了好多书，一直没找到水汽通量散度的公式，急用啊，谢谢

奚子惠

金玉姻缘 发表于 2015-6-2 11:04
楼主你好，我要计算地面到500hpa整层水汽通量散度，都需要哪些变量啊，公式在哪里呢，我查了好多书，一直没 ...

需要u、v、比湿shum、还有地表气压pres  家园上有很多关于水汽通量的帖子可以参照一下http://bbs.06climate.com/forum.php?mod=viewthread&tid=4920