MetPy get_layer仅返回前3个气压值？mean_pressure_weighted存疑

get_layer的默认行为符合预期 #

get_layer函数默认会提取从最底层（最高气压值）到850 hPa的气压层。你的气压数组按1000hPa→10hPa递减排列，1000、950、900hPa均高于850hPa的阈值，因此被选中返回，这完全符合函数的设计逻辑。

如果需要提取其他范围的气压层，只需显式指定bottom和top参数，比如要获取全部气压层：

full_layer = get_layer(pre, bottom=np.min(pre), top=np.max(pre))

mean_pressure_weighted的计算范围关联 #

mean_pressure_weighted默认会调用get_layer的默认结果作为计算层，这就是为什么你的计算结果和手动取前3个值的积分结果一致——它只用到了get_layer返回的前3个气压层。

要基于全部气压层计算加权平均，需要显式传递完整的气压序列作为layer参数：

# 直接传入完整气压序列
weighted_mean = mean_pressure_weighted(pre, tem, layer=pre)

# 或通过get_layer指定完整范围后传入
full_layer = get_layer(pre, bottom=np.min(pre), top=np.max(pre))
weighted_mean = mean_pressure_weighted(pre, tem, layer=full_layer)

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