我同时使用 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 中使用的值,但错误无法解决。我不确定错误是否来自函数源代码中的错误,或者我是否忽略了输入中的错误,无论是干失效还是湿失效。