1

我同时使用 dry_lapse 和 wet_lapse 来尝试查找数据集的 showalter 索引。我能够让 dry_lapse 成功运行,并希望在 wet_lapse 函数中使用 dry 的温度输出。目前,我的代码如下所示:

# Stability
dl = mpcalc.dry_lapse(p[2], tac[0], p[0])
dl = (dl.magnitude - 273.15)*units.degC # Change units to C

ml = mpcalc.moist_lapse(p[11], dl, p[2]) 

其中 p 是数据中的压力水平,tac 是地块剖面温度,以摄氏度为单位。

代码的目的是计算包裹温度,如果它从 p[0] (1008mb) 绝热提升到 p[2] (950mb) 处的 lcl 水平,然后计算包裹温度,如果它随后被提升潮湿从 p[2] 到 [p11](500mb 级)绝热。我已通过将它们替换为 dry_lapse 中使用的值来检查我在 mot_lapse 中使用的值,但错误无法解决。我不确定错误是否来自函数源代码中的错误,或者我是否忽略了输入中的错误,无论是干失效还是湿失效。

4

1 回答 1

1

问题来自使用旧版本的 MetPy,通过在本地更新我的版本解决了这个问题。

于 2020-12-31T16:17:16.333 回答