Python 涡度平流、温度平流

JessicaYang · 发表于 2020-3-13 13:43:27

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

用python mpcalc.advection 函数计算平流，总是出错，请大神指点

data = xr.open_dataset('/home/ymx/文档/testpy/data/ecmwf.nc')
data = data.metpy.parse_cf().squeeze()
#print(data)

#data_crs = data['t'].metpy.cartopy_crs

lon    = data['longitude']
lat    = data['latitude']
time    = data['time']
#print(time)

u       =  data['u']
v       =  data['v']

# calculate Vorticity Advection
dx, dy = mpcalc.grid_deltas_from_dataarray(u)
#dx, dy = mpcalc.lat_lon_grid_deltas(lon, lat)
vor    = mpcalc.vorticity(u, v, dx, dy)

vort_adv = mpcalc.advection(vor, [u, v], (dx, dy), dim_order='yx') * 1e9  上面都是对的，就这句出错

坎坷 · 发表于 2022-5-12 22:05:01

u_level = units.Quantity(np.array(u_level), u_level.units)
v_level = units.Quantity(np.array(v_level), v_level.units)

lat = np.array(lat)
lon = np.array(lon)

dx, dy = mpcalc.lat_lon_grid_deltas(lon, lat)

vor_var = mpcalc.vorticity(u=u_level, v=v_level, dx=dx, dy=dy)
vor_adv_var = mpcalc.advection(vor_var, u=u_level, v=v_level, dx=dx, dy=dy)
这是我写的代码

xinshou6 · 发表于 2020-3-14 10:15:43

亲你可以给个错误或者写个公式我外行

笺疏 · 发表于 2020-8-22 15:19:30

请问问题解决了吗

陸陆 · 发表于 2021-3-29 11:52:19

可以尝试改为vort_adv = mpcalc.advection(vor, u, v, dx, dy,) * 1e9 上面都是对的，就这句出错

田甜恬 · 发表于 2021-9-15 22:12:13

请问楼主找出正确做法了吗？

LIb · 发表于 2021-11-11 19:35:50

请问楼主有没有完整的程序吗？可以分享给学习一下吗？

biubiubiuemm · 发表于 2022-2-22 00:03:07

请问楼主解决了嘛

风策策度庭竹 · 发表于 2022-3-28 22:52:15

中文也太搞了吧识别不了

坎坷 · 发表于 2022-5-12 21:12:16

涡度平流是不是地转涡度的平流呀，你用的相对涡度的公式

坎坷 · 发表于 2022-5-12 21:43:08

不好意思，刚刚百度了一下，涡度平流的概念还是没查清楚，不过大部分说的是相对涡度的平流，我刚刚说错了

		自动登录	找回密码
密码			立即注册

[求助] Python 涡度平流、温度平流

登录后查看更多精彩内容~

浏览过的版块