前言
最近在计算一下涡度的垂直项时，发现一个问题，如果使用笛卡尔坐标系下的涡度计算公式，会发现对于我们从ERA5下载的垂向速度ω，它的单位是：pa/s，而对于水平的速度u、v来说，它们的单位是：m/s，这在计算速度的旋度时就会出现问题，由于单位不一致引起计算无法进行，对于后续数据处理造成干扰。
本来我在网上搜寻解决方案时，发现大部分的解答都只提供了解决思路（由公式转换），往往比较冗杂，不够简单。所辛，我在metpy这个工具包中发现了包装好的函数，只需要提供相应的数据，即可成功转换！
这里展示由压力坐标向高度坐标转换的例子，（metpy也提供了反向转换的函数，可以尝试）
函数如下：官方链接在这
metpy.calc.vertical_velocity(omega, pressure, temperature, mixing_ratio=0)
需要传入三个参数：压力坐标系下的垂向速度ω、压力p、气温t，混合比（默认为0）
omega (pint.Quantity) – Vertical velocity in terms of pressure pressure (pint.Quantity) – Total atmospheric pressure temperature (pint.Quantity) – Air temperature mixing_ratio (pint.Quantity, optional) – Mixing ratio of air
还是比较容易理解的，下面给出了一个示例，仅供参考：
import xarray as xr import numpy as np from metpy.units import units import metpy.constants as constants import metpy.calc as mpcalc ######################################################### path1=r'D:\PV\down_hourly.nc' f = xr.open_dataset(path1).sortby('latitude') path3=r'D:/PV/tem.200406.nc' f_t=xr.open_dataset(path3).sortby('latitude').transpose('time','level','latitude','longitude') ######################################################### t=f_t.t f = f.transpose('time','level','latitude','longitude') f_u = f.u f_v = f.v f_w = f.w ######################################################### lat = f_u.latitude lon = f_u.longitude time = f_u.time lev = f_u.level pressure = lev*units.hPa ######################################################### #height=mpcalc.pressure_to_height_std( pressure) w = mpcalc.vertical_velocity(f_w, lev, t, mixing_ratio=0)
计算之后发现w的单位已经转换成功了：
有兴趣的小伙伴快去试试吧~
