请问湿位涡MPV2的计算公式

绿裙子

搜索了坛子里计算湿位涡的帖子，而我自己这里也有一份计算湿位涡的gs文件，仔细比较之后，发现两个gs计算湿位涡公式还是有些差异，其中计算MPV1的都没有问题，但是MPV2差异比较大，所以发帖来问问MPV2的计算公式，究竟哪一种方法可取，欢迎大家来讨论
注：公式中eqt为假相当位温，u，v表风；lev是hPa层；
公式1（坛子里找来的）
'define pi=3.14159'
'define g=9.8'
'define R=6.371e6'
'define dp=100*(lev(z-1)-lev(z+1))'
'define dv=v(z-1)-v(z+1)'
'define du=u(z-1)-u(z+1)'
'define dx=2.0*R*cos(lat*pi/180.0)*pi/180.0'
'define dy=2.0*R*pi/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define mpv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
公式2
'define pi=3.14159'
'define g=9.8'
'define R=6.371e6'
'define deqt=eqt(z-1)-eqt(z+1)'
'define dp=100*(lev(z-1)-lev(z+1))'
'define dy=cdiff(lat,y)*pi/180*R'
'define dx=cdiff(lon,x)*pi/180*cos(lat*pi/180)*R'
'define du=cdiff(u,y)'
'define dv=cdiff(v,x)'
'define mpv2=g*(dv/dx-du/dy)*deqt/dp'

qxtlyf

抢沙发！！！！！顶一顶，画一张看看差别

绿裙子

qxtlyf 发表于 2014-2-21 17:24

                               
登录/注册后可看大图

抢沙发！！！！！顶一顶，画一张看看差别

公式1画出的剖面图 D:\OpenGrADS\Contents\Cygwin\Versions\2.0.a7.oga.3\i686\2011061720-1.gif

绿裙子

这是公式2画的剖面图D:\OpenGrADS\Contents\Cygwin\Versions\2.0.a7.oga.3\i686\2011061720-2.gif

839

对照公式看看不就得了，这问题太简单，不思考容易懒惰

绿裙子

两个gs中对dx、dy的处理有不同，为什么？不懂啊，求教导

绿裙子

顶顶，求指教

north23

求指导啊      

qixiangren34

我也发现了这个问题，我一般用第二个，觉得第一个不妥。

Qoo

学习下~~呵呵