利用metpy画t-lnp图报错

爱地理不爱物理

基于metpy画探空图时出现以下问题，其中2023050408是08时的探空资料，0514是利用14时气温和露点订正的资料，我在利用14时数据画图时没有报错，但是利用08时数据时出了这个问题，有没有大神帮忙解答下，万分感谢！图是能显示的，就是在利用My_Class包计算时出现了问题

edwardli

【metpy版本：1.4.1】

全部用你的代码，未报任何错误。都能正常出图。


你的报错，可能是：气压层不单调递增or递减。

edwardli

We've seen issues like this before with high-resolution data (alluded to in #1339). The problem is that a few levels decrease by ~0.1 hPa and break some assumptions in the code. A work-around is usually to pass the pressure levels through scipy.signal.medfilt to smooth it out ever so slightly.这是官方原文回答。

edwardli

So as we start to see 1Hz soundings in the wild, we start to see issues with our calculations and the fact that sometimes pressure, as directly reported, is not necessarily monotonically decreasing. For instance, see this sample sounding. There are a couple of spots where the pressure bounces up for a single point by 0.1 hPa. This is not significant, but it does cause things like scipy.integrate.odeint to get angry when integrating with pressure as the independent variable--thus causing cape_cin to bomb out on these soundings.  For these soundings, a med_filt is enough to solve the problem, but I'm not sure that works as a general solution. I'm not really sure what a good general solution is.