求助：计算整层水汽通量散度和文献里差别较大

张一鸣

1.问题描述：我想用ERA5月均数据计算多年季平均的水汽通量散度，尝试了几个论坛上的不同的计算顺序（a.先求各层散度，让后在从近地面积分至300；b.先从近地面积分qv和qu，最后算一个散度）感觉效果都不是很好。注：数据的维度：X480；Y240；T12；Z20

2.疑惑部分：
  a.看到大家算出来的0值区域（白色）占了很大部分，而我的大部分却不是0。希望能指点一下，感谢！！！！

3.对比图（先积分后求散度）：（左图是我是用grads画的，右图是文献中给出的图）

4.代码部分（先积分后求散度）：
'reinit'
'sdfopen G:\Data\all_div_v3\sf1.nc'
'define p1=ave(sp,t=1,t=12)'
'close 1'

'set grads off'
'set grid off'
'set parea 1.5 10.0 1.5 8.0'
'set map 1 1 5'
'set xlint 20'
'set ylint 20'
'set xlopts 1 5 0.18'
'set ylopts 1 5 0.18'
'set clopts 1 5 0.18'
'set mpdset cnworld cnriver'

'set font 4'

'sdfopen G:\Data\all_div_v3\quv1.nc'
'set z 1 20'
'define u=ave(u,t=1,t=12)'
'define v=ave(v,t=1,t=12)'
'define q=ave(q,t=1,t=12)'

'define uq=u*q'
'define vq=v*q'

'set z 1'
'define qu=vint(p1,uq,300)/1000'
'define qv=vint(p1,vq,300)/1000'
'define uvq=mag(qu,qv)'

'define divq=hdivg(qu,qv)*1e6'

'set lon 60 180'
'set lat 0 60'

'set rgb 20 24 24 112'
'set rgb 21 16 78 139'
'set rgb 22 23 116 205'
'set rgb 23 72 118 255'
'set rgb 24 91 172 237'
'set rgb 25 173 215 230'
'set rgb 26 209 237 237'
'set rgb 27 229 239 249'
'set rgb 28 242 255 255'
'set rgb 29 255 255 255'
'set rgb 30 253 245 230'
'set rgb 31 255 228 180'
'set rgb 32 243 164 96'
'set rgb 33 237 118 0'
'set rgb 34 205 102 29'
'set rgb 35 224 49 15'
'set rgb 36 237 0 0'
'set rgb 37 205 0 0'
'set rgb 38 139 0 0'

'set csmooth on'
'set clevs   -20 -16 -12 -8 -4 -2 2 4 8 12 16 20'
'set rbcols 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 '

'd re(divq,1)'
'cbarn 1 2'
'draw title 1000-300hPa divigence of water vapor flux'
'printim sdfopen G:\Data\all_div_v3\test_.png  white x1920 y1080'

膘膘

要用日资料计算

相思豆

学习学习！！

洞里骑象

我用NCL算的时候积分前还乘了dp/g这一项，dp是各层的厚度。你的代码里好像没乘这一项。

张一鸣

洞里骑象 发表于 2023-3-19 19:05
我用NCL算的时候积分前还乘了dp/g这一项，dp是各层的厚度。你的代码里好像没乘这一项。

感谢回复，就是不太理解为什么要乘以这个dp呢。

洞里骑象

张一鸣 发表于 2023-3-20 09:09
感谢回复，就是不太理解为什么要乘以这个dp呢。

如果是计算地面积分到300hPa高度处的水汽通量散度的话，计算方法应该是图片里写的那样？红框里的就是dp/g

张一鸣

洞里骑象 发表于 2023-3-20 09:42
如果是计算地面积分到300hPa高度处的水汽通量散度的话，计算方法应该是图片里写的那样？红框里的就是dp/g

谢谢大佬，我之前参考的是天气学原理第四版里的公式。这两个有什么区别吗

洞里骑象

张一鸣 发表于 2023-3-20 20:51
谢谢大佬，我之前参考的是天气学原理第四版里的公式。这两个有什么区别吗

层高1hPa的话好像确实没区别了，抱歉。顺便，我上面发的那张图上下限有点问题，好像应该是下300上1000.

张一鸣

洞里骑象 发表于 2023-3-22 11:52
层高1hPa的话好像确实没区别了，抱歉。顺便，我上面发的那张图上下限有点问题，好像应该是下300上1000.

!!!!!谢谢大佬提醒，原来积分的时候要乘每层的层厚度啊，感谢

紫藤萝

ERA5不是提供了整层积分的水汽通量散度吗 Vertical integral of divergence of moisture flux，为什么还要重新计算