疑问：700hpa水汽通量散度

气象局

我根据清风http://bbs.06climate.com/forum.php?的帖子，做了修改，程序如下：
* This is a script for displaying moisture convergence
* Written by Michael Maxwell
*   
* rh    = Relative Humidity in %
* t     = Temp at *set level in degrees Kelvin
* tc    = Temp in degrees C
* td    = Dewpoint at *set level in degrees C
* e     = Vapor pressure
* mixr  = Mixing ratio
* u     = U-wind in m/s
* v     = V-wind in m/s
* mconv = moisture convergence/divergence. convergence is positive and divergence is negative.
'reinit'
'sdfopen e:\zl\2012t.nc'
'sdfopen e:\zl\2012rhum.nc'
'sdfopen e:\zl\2012u.nc'
'sdfopen e:\zl\2012v.nc'
'set lev 700'
'set lat 23 30'
'set lon 101 111'
'set t 5'
'set mpdset cnworld'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set grads off'
'set timelab off'
'set grid off'
'define tc=(t.1-273.16)'
'define td=tc-((14.55+0.114*tc)*(1-0.01*r.2)+pow((2.5+0.007*tc)*(1-0.01*r.2),3)+(15.9+0.37*tc)*pow((1-0.01*r.2),14))'
'define vapr=6.112*exp((17.67*td)/(td+243.5))'
'define e=vapr*1.001+(lev-100)/900*0.0034'
'define mixr=0.62137*(e/(lev-e))*1000/9.8'
#通过此gs文件算出的水汽通量的单位是g/(cm.s),一般为 十几到几十
'define qx=u.3*mixr'
'define qy=v.4*mixr'
'define mconv=hdivg(qx,qy)*1e6'
'enable print e:\tx\2lb-shuiqi900.gmf'
'set arrscl 1 500'
*'set lon 124'
*'set cmin 30'
'set gxout shaded'
'd mag(qx,qy)'
*'set cmin 30'
'set gxout contour'
'd mag(qx,qy)'
*'d theta'
*'d skip(u.3*2.5,4);skip(v.4*2.5,4)'
'print'
'disable print'
'printim e:\tx\2lb-shuiqi900.png white'
'reinit'
得到的图为：

修改程序得到的图片

我觉得不太对，请大家帮我看看。
我想得到下图：

这个程序能做到不？请大家帮我分析分析，我已经进入死胡同了！

sun_shine_Xia

LZ，你上面那个图画的是水汽通量散度么？我看你d的是mag(qx,qy)，这不水汽通量的大小么？

river

去看一下实用手册吧，先了解一下什么是等值线图，什么是阴影图，什么是矢量图，还有去看看mag函数是求什么的，求散度和涡度这类常用的函数又是哪几个。这些东西搞不清楚，那你就得不到你要的东西

zhuzi1009

谢谢楼主分享，好人呐

气象局

zhuzi1009 发表于 2013-11-12 13:05

                               
登录/注册后可看大图

谢谢楼主分享，好人呐

没有了，你看错了，我是提问的哦！我自己也在学习！

气象局

river 发表于 2013-11-11 18:30

                               
登录/注册后可看大图

去看一下实用手册吧，先了解一下什么是等值线图，什么是阴影图，什么是矢量图，还有去看看mag函数是求什么的 ...

怎么说了，我也看了的啊，但是只能说走一步学一步问一步吧，毕竟不是在学校！还是谢谢你啊，从你的发的帖子里我学了很多东西！

气象局

sun_shine_Xia 发表于 2013-11-11 17:37

                               
登录/注册后可看大图

LZ，你上面那个图画的是水汽通量散度么？我看你d的是mag(qx,qy)，这不水汽通量的大小么？

是啊！但是我试着改了上面的d都得的不到图，我知道mag是去平均的哦，是一个值，但是我要的图应该怎样画了？

sun_shine_Xia

气象局 发表于 2013-11-12 13:21

                               
登录/注册后可看大图

是啊！但是我试着改了上面的d都得的不到图，我知道mag是去平均的哦，是一个值，但是我要的图应该怎样画了 ...

mag求的是向量的模的大小，你要求的水汽通量散度应该是你gs文件里面的mconv变量~你d mconv试试

kaidada

river 发表于 2013-11-11 18:30
去看一下实用手册吧，先了解一下什么是等值线图，什么是阴影图，什么是矢量图，还有去看看mag函数是求什么 ...

给点实用的好吧，就知道训人

river

kaidada 发表于 2014-9-17 18:29
给点实用的好吧，就知道训人

我想问一下你给了什么实用的东西？我只是提一个建议，如果基础搞不懂，别人再怎么说也不会明白。难道要其他人直接给一个答案来？