Metpy计算Brunt-Vaisala频率返回NaN值问题求助

核心问题分析 #

你的代码存在两个关键错误，直接导致结果全为NaN：

• 位温计算未用绝对温度：位温的定义要求基于开尔文（绝对温度），但你直接使用摄氏温度计算，导致位温数值偏差，进而出现位温随高度递减的不稳定层结，此时根号内为负数，计算结果返回NaN。
• 未明确物理量单位：Metpy依赖单位系统保证计算准确性，你未给气压、温度、高度标注单位，存在潜在的单位转换风险。

修正后的代码 #

import numpy as np
import metpy.calc as mpcalc
from metpy.units import units

# 带单位的原始观测数据
el = np.array([75, 1015.09, 1786.01, 2607.85, 2834.74]) * units.meter
pres = np.array([945.79, 889.87, 810.35, 732.69, 712.33]) * units.hPa
temp = np.array([11.27, 8.03, 7.8, 2.09, 0.28]) * units.degC

# 用Metpy内置函数计算位温，自动处理单位和公式
pt = mpcalc.potential_temperature(pres, temp)

# 计算Brunt-Vaisala频率
bv = mpcalc.brunt_vaisala_frequency(el, pt)
print(bv)

关键说明 #

• 优先使用Metpy内置函数：mpcalc.potential_temperature会自动完成摄氏温度转开尔文、绝热过程计算，避免手动公式出错。
• 稳定层结要求：Brunt-Vaisala频率仅在稳定层结（位温随高度递增）时返回实数值；若仍出现NaN，说明对应层为不稳定层结（位温随高度降低），属于正常物理现象。

内容的提问来源于stack exchange，提问作者Sadie