水汽通量散度剖面

酒国

我做的是水汽通量散度剖面，思路如下：先一层一层算出来并输出为grd文件，写相应的ctl，最后编写gs画图。
一层一层的水汽通量散度可以画出来，但剖面图感觉就不对了，不知道哪里出错了，希望大家帮忙看看，我把相关的程序和图附在下面。
一层一层的gs:
'reinit'
'open d:\test\20110923\fnl_20110923_00_00.ctl'
'set fwrite D:\test\20110923\q.grd'
'set gxout fwrite'
'set lon 108 120'
'set lat 12 24'
'set t 23'
'set lev 1000'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'gxprint'
'printim D:\test\20110923\1000.png white'
'disable gxprint'
'set lev 925'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'set lev 850'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'set lev 700'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'set lev 600'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'set lev 500'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'set lev 400'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'set lev 300'
'define pre=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*RHprs/100'
'define qu=(UGRDprs*q/9.8)'
'define qv=(VGRDprs*q/9.8)'
'define qdiv=hdivg(qu,qv)*1e6'
'd qdiv'
'disable fwrite'
相应的ctl文件：
DSET D:\test\20110923\q.grd
undef -9.99E33
xdef 12 linear 108 1
ydef 13 linear 12 1
zdef 8 levels 1000 925 850 700 600 500 400 300
tdef 1 linear 12z28Sep2011 6hr
vars 1
qdiv 8 99
Endvars
画图的gs:
'reinit'
'open D:\test\20110923\qdiv.ctl'
'set lat 14 20'
'set lon 114'
'set z 1 6'
'set t 1'
'd qdiv/1000000'
单层的水汽通量散度图：

水汽通量散度剖面：

四叶草

建议你用循环提取数据试试

river

请论坛搜索 剖面图 看看别人怎么画的

hm_others

建议每一层都单独设置一下xy和gxout fwrite

酒国

river 发表于 2015-11-30 00:04
请论坛搜索 剖面图 看看别人怎么画的

谢谢，看到了这篇http://bbs.06climate.com/forum.php?mod=viewthread&tid=912。。。，觉得自己把这剖面图弄复杂了

river

酒国 发表于 2015-11-30 15:43
谢谢，看到了这篇http://bbs.06climate.com/forum.php?mod=viewthread&tid=912。。。，觉得自己把这剖面 ...

搞得不是一般的复杂。提问前请先搜索，常用的到问题论坛会有很多相关的帖子的

aas314

好东西学习了

iam428jll

谢谢分享！！！

sunflower_葵

谢谢分享~~~~~~~~

taiboyu

前来学习