Python值错误:操作数不能与形状(732,25,17,33) (25,732,17,33)一起广播-x数组坐标切换？

Question

我试图从metpy.calc.dewpoint_from_specific_humidity(pressure，网格数据中计算露点温度，用元函数ERA5温度、比湿度)，但是我得到了一个错误：

ds = xarray.open_dataset('/dir/file.nc')
Td = metpy.calc.dewpoint_from_specific_humidity(ds.pressure_level, ds.temperature, ds.specific_humidity)

ValueError: operands could not be broadcast together with shapes (732,25,17,33) (25,732,17,33) 

形状是x数组的坐标(时间:732，压力级:25，纬度:17，经度:33)。它似乎是顺序切换和函数无法计算，因为切换的尺寸。有办法解决这个问题吗？

如果有人能帮我，那就太好了！

最好，莉娜

Answer 1

(张贴答案以确保完整。)我不知道为什么数据会翻转维度，但是在早期版本的MetPy中，我们严格依赖NumPy广播来处理数组的匹配--这并不理解数组，它只能扩展数组来尝试匹配。对于较新版本的MetPy (我认为1.0，但肯定> 1.0)，它使用X阵列广播，它可以利用命名的维度对齐用于操作的数组。