metpy垂直速度单位问题

Question

我试着从垂直压力速度估计垂直风速。我找到了metpy.calc.vertical_velocity(ω，压力，温度，mixing=0)。

https://unidata.github.io/MetPy/latest/api/generated/metpy.calc.vertical_velocity.html

但是，没有关于它需要什么单元的文档。我以Pa s-1表示欧米加，以Pa表示的压强，以K表示的温度，都以xarray格式表示。我希望输出也是一个xarray。我得到了这样的错误：“无法从'pascal * second‘转换为’无量纲‘。有没有想过为什么它需要无量纲变量，以及如何解决这个问题？

Answer 1

该错误来自于解析Udunit格式的单元字符串"Pa 1“时出现的问题。pint (MetPy使用的单元库)中的默认单元解析器将其解析为：Pa * s - 1，这是无意义的。幸运的是，在MetPy 1.0中，这应该是可行的：

import metpy.calc
import xarray as xr

omega = xr.DataArray([5], attrs={'units': 'Pa s-1'})
pressure = xr.DataArray([50000], attrs={'units': 'Pa'})
temperature = xr.DataArray([283], attrs={'units': 'K'})
metpy.calc.vertical_velocity(omega, pressure, temperature, mixing_ratio=0)

您目前可以从pip安装metpy 1.0.0rc1：

pip install --pre metpy

或者来自conda-forge：

conda install -c conda-forge/label/metpy_rc metpy=1.0.0rc1