Metpy Brunt频率特性将与NaN一起返回

Question

我使用metpy计算Brunt频率，代码运行良好，但我的数组值作为一组NaN值返回。我在以下高度测量了压力和温度:75米，1015.09m，1786.01m，2607.85m，2834.74m。有人能告诉我怎么解决这个问题吗？也许我对Brunt Vaisala的概念理解是错误的。

el=np.array([75,1015.09,1786.01,2607.85,2834.74])
pres=np.array([945.79,889.87,810.35,732.69,712.33])
temp=np.array([11.27,8.03,7.8,2.09,0.28])
PT=np.array(temp*(1000/pres)**(.286))
BV=metpy.calc.brunt_vaisala_frequency(el*units.meter,PT*units.kelvin)

Answer 1

计算Brunt频率需要取位温相对于高度的垂直导数。由于您传递的是一个完全一致的高度数组(如el)，垂直导数的值是未定义的，并返回为nan。这将导致整个计算的结果也作为nan返回。

编辑更新的数据：

您发布的新数据的潜在温度值随着高度的增加而降低。根据定义，这是一个不稳定的概要文件，在这种情况下是首当其冲的-瓦塞拉频率是未知的。 --因此对于这些值，MetPy返回nan。

如果您想看到负值，可以调用metpy.calc.brunt_vaisala_frequency_squared，这样可以避免出现问题的平方根。