爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 30599|回复: 10

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

[复制链接]

新浪微博达人勋

发表于 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  上面都是对的,就这句出错

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 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)
这是我写的代码
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2020-3-14 10:15:43 | 显示全部楼层
亲你可以给个错误或者写个公式我外行
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-8-22 15:19:30 | 显示全部楼层
请问问题解决了吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-29 11:52:19 | 显示全部楼层
可以尝试改为vort_adv = mpcalc.advection(vor, u, v, dx, dy,) * 1e9  上面都是对的,就这句出错
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-9-15 22:12:13 | 显示全部楼层
请问楼主找出正确做法了吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-11-11 19:35:50 | 显示全部楼层
请问楼主有没有完整的程序吗?可以分享给学习一下吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-2-22 00:03:07 | 显示全部楼层
请问楼主解决了嘛
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-3-28 22:52:15 | 显示全部楼层
中文也太搞了吧 识别不了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-5-12 21:12:16 | 显示全部楼层
涡度平流是不是地转涡度的平流呀,你用的相对涡度的公式
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-5-12 21:43:08 | 显示全部楼层
不好意思,刚刚百度了一下,涡度平流的概念还是没查清楚,不过大部分说的是相对涡度的平流,我刚刚说错了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表