如何计算垂直高度剖面

Question

我想从ECMWF数据计算层厚度。

我拥有的是变量[[T（在一层内）和relhum（在一层内）和压力（在层界面和层中点处）。

我在理解pressure函数的metpy.calc.thickness_hydrostatic_from_relative_humidity参数时遇到麻烦：

[当我查看文档字符串中给出的公式时，

$$ Z_2-Z_1 =-\ frac {R_d} {g} \ int_ {p_1} ^ {p_2} T_v d \ ln p $$

在我看来，该函数似乎将返回LHS，即$ Z_2-Z_1 $。

但是，在我看来[pressure，temperature，relative_humidity自变量都必须具有相同的尺寸。

我发现这令人困惑：为了使层的厚度介于$ Z_1 $和$ Z_2 $之间，我希望同时输入压力$ p_1 $和$ p_2 $。但是，通常将温度定义为满水平，即我的温度比压力小。

例如，为了计算底层的厚度，我希望给出表面压力，层顶部压力以及层内温度（和湿度）。但是当我尝试使用这样的功能时，我得到了

ValueError: operands could not be broadcast together with shapes (361,1440,79) (361,1440,78)

请帮助我了解如何正确使用此功能。

我想根据ECMWF数据计算层厚度。我所拥有的是变量T（在一层内）和relhum（在一层内）和压力（都在层界面和层上）。

Answer 1

但是，在我看来[pressure，temperature，relative_humidity自变量都必须具有相同的尺寸。