求python计算沙氏指数SI、cape值、抬升指数LI、K指数

小蜗

求大神分享一下python计算沙氏指数SI、cape值、抬升指数LI、K指数等源代码

Om_MyGad

插眼 si值cape 抬升凝结 k指数

niecj09

参考metpy

feiji158

https://unidata.github.io/MetPy/latest/examples/
例如：

1. data_at_point = data.metpy.sel(
   
2.     time1='2017-09-05 12:00',
   
3.     latitude=30 * units.degrees_north,
   
4.     longitude=260 * units.degrees_east
   
5. )
   
6. dewpoint = mpcalc.dewpoint_from_relative_humidity(
   
7.     data_at_point['Temperature_isobaric'],
   
8.     data_at_point['Relative_humidity_isobaric']
   
9. )
   
10. cape, cin = mpcalc.surface_based_cape_cin(
    
11.     data_at_point['isobaric3'][::-1],
    
12.     data_at_point['Temperature_isobaric'][::-1],
    
13.     dewpoint[::-1]
    
14. )
    
15. cape

复制代码

我爱气象111

楼主 你现在问题解决了吗

我爱气象111

楼主问题解决了吗

球求大人的球球

feiji158 发表于 2021-11-9 15:58
https://unidata.github.io/MetPy/latest/examples/
例如：

计算出来的cape和micaps里自带的cape差别有点大

feiji158

球求大人的球球 发表于 2024-12-23 23:14
计算出来的cape和micaps里自带的cape差别有点大

哈哈，我也是如此

球求大人的球球

feiji158 发表于 2024-12-26 16:18
哈哈，我也是如此

和实况对比了以下，好像是mpcalc.dewpoint_from_relative_humidity()计算的各层dewpoint 和micaps差的大。

啥也没研究出来

wrf-python可以只是读cape_2d，直接输出cape指数和cin指数