使用ERA5数据和metpy计算抬升指数

bsb0

下载了ERA5 hourly data on pressure levels的数据，选择了整层的数据，即1000hPa到1hPa，想使用mpcalc.lifted_index计算抬升指数，公式中需要先计算parcel profile（翻译为层结？），但是会报错，ValueError: place: mask and data must be the same size，请教这是什么原因呢？
PS：网上搜了一下，有人说需要将array数据转换为 Unit类型的数据；不转就会出错。使用.data*units.degC也还是会报错，不明白原因。代码如下

1. ds = xr.open_dataset('D:/history_weather/CDS20210500h.nc')
   
2. lon=ds['longitude'][:]#读取经度
   
3. lat=ds['latitude'][:]#读取纬度
   
4. lon_range = lon[(lon>117.4) & (lon<123.6)]
   
5. lat_range = lat[(lat>24.9) & (lat<32.4)]
   
6. lons,lats=np.meshgrid(lon_range,lat_range)
   
7. p = ds['level']
   
8. t = ds['t']
   
9. q = ds['q']
   
10. h = ds['r']
    
11. 
    
12. ##计算部分
    
13. e = mpcalc.vapor_pressure(p,q)
    
14. Td = mpcalc.dewpoint(e)
    
15. t0 = ds['t'][0,36]
    
16. td0 = Td[36,0]
    
17. prof = mpcalc.parcel_profile(p,t0,td0)
    
18. li = mpcalc.lifted_index(p,t,prof)

复制代码

bsb0

自己顶一下，有哪位大神可以指导一下

扫眉才子

请问楼主解决了吗？求经验

bsb0

扫眉才子 发表于 2021-10-21 15:00
请问楼主解决了吗？求经验

莫有，尝试过把某一点的高空资料和地面资料拼接成一个数组，但是也失败了