如何计算了指数从ERA5压力水平数据?

Question

我采取了一个雷雨情况我需要计算指数,但我没有得到任何代码相关的稳定性指数在Python语言,代码可在毕业生,NCL对我来说很难理解。我用MetPy代码但是他们不是计算,KI和其他指标……提前谢谢下面的代码是错误的,我知道我使用了公式,但告诉我正确的方法。

nc1 =数据集(pl_120415.nc) lat1 = np.array (nc1.variables['自由'])lon1 = np.array (nc1.variables['经度'])列弗= np.array (nc1.variables['水平'])时间= nc1。('时间']变量r = nc1。变量[' r '] [:,:, 25:28 20:25] t = nc1。[' t ']变量u = nc1。变量[' u '] [:,:, 25:28 20:25] v = nc1。变量[' v '] [:] w = nc1。[' w ']变量t_500 = nc1。变量[' t '][21日::,:# # #领域空气温度在500 Pressure_level tp_500 = nc1。变量[' t '][: 21日25:28,20:25)# #在specfific位置温度(包裹温度)李= np.subtract (t_500、tp_500 dtype = np.float32)打印(“李”,李)

Answer 1

Metpy 1.0有一个lifted_index函数,它需要一个剖面温度参数,进而可以使用计算parcel_profile函数。需要有露点温度,你要么必须直接从另一个含水量或计算变量。我不确定你什么r变量,所以我不确定是否能帮你计算。而不是使用numpy数组,使用xarray处理netcdf输入变量往往是更容易和更有效。