水汽通量等的垂直加权积分的函数

yaya4853

对于整层水汽通量及散度的垂直积分，是用np的trapz函数，还是用metpy的mean_pressure_weighted的函数，感觉两个函数都有点不太对。请大神们指教，谢谢。

灭火器

我是直接用的DataArray的integrate方法（底层就是np.trapz），跟再分析资料比较了一下结果还行。但是没考虑去掉地表以下部分

yaya4853

灭火器 发表于 2022-3-23 19:40
我是直接用的DataArray的integrate方法（底层就是np.trapz），跟再分析资料比较了一下结果还行。但是没考虑 ...

好的，谢谢。那我就还是用trapz函数好了

灭火器

1. g = 9.8
   
2. a = 6371E3
   
3. ds['level'] = ds.level * 100  # 气压单位转为Pa.
   
4. tcwv = ds.q.integrate('level') / g
   
5. uq = ds.u * ds.q
   
6. vq = ds.v * ds.q
   
7. coslat = np.cos(np.deg2rad(ds.latitude))
   
8. duqdlon = np.rad2deg(uq.differentiate('longitude'))
   
9. dvqdlat = np.rad2deg(vq.differentiate('latitude'))
   
10. div = (duqdlon / coslat + dvqdlat) / a
    
11. mfd = div.integrate('level') / g

复制代码

补个代码，我用ERA5数据算的，不知道有没有懂这个的指出问题

一大碗年糕

求问一下直接调用函数垂直积分和转化为求和来算积分差距大吗？看有的帖子里说很多都是求和来算

yaya4853

灭火器 发表于 2022-3-23 20:17
补个代码，我用ERA5数据算的，不知道有没有懂这个的指出问题

注意era5的u和v是有正负的，它是有风向分量的

灭火器

yaya4853 发表于 2022-3-24 09:08
注意era5的u和v是有正负的，它是有风向分量的

啊？本来理论公式里风速就有正负

yaya4853

灭火器 发表于 2022-3-24 09:26
啊？本来理论公式里风速就有正负

天气学原理原理里面的水汽通量公式里面风速可是带着绝对值符号的

puck66

楼主后续如何？得空给分享一下咋做的嘛