python绘制任意两点剖面

metpy中的cross_section提供了非常便捷的绘制剖面的方法

具体可见官网：Cross Section Analysis — MetPy 1.1 (unidata.github.io)
 

如果你只需要简单地画个剖面，好比第一张图，那么请看下面的代码：

# %%
import xarray as xr
import proplot as pplt
from metpy.interpolate import cross_section

# %%
era5fnm = r'F:/era5/era5_2021071900_2021072023/era5_2021071900UTC.nc'
era5f = xr.open_dataset(era5fnm)
era5f = era5f.metpy.parse_cf().squeeze()

# %%
# create cross section
start = (12.0, 96.0)
end = (27.0, 120.0)
cross = cross_section(era5f, start, end)

# %%