我怎样才能得到这个metpy计算的有效风暴螺旋度和有效切变值？

Question

我正在尝试使用metpy.calc 函数来获取supercell 复合值，如下所示： Supercell Composite

但是，我似乎在任何地方都找不到如何计算有效风暴螺旋度和有效剪切力。 我的 grib2 数据中有风暴相对螺旋度，但我如何获得有效的？

谢谢！

Answer 1

从 1.0 版开始，MetPy 不具备计算有效 SRH (ESRH) 和有效体积剪切 (EBS) 的功能。

来自Thompson 等人的释义。 2007 年：

ESRH 是在“有效风暴流入层”内计算的 SRH，它是由地块特征 CAPE >= 100 J/kg 和 CIN >= -250 J/kg 的下边界定义的风暴的最低垂直区域，继续直到不再满足其中任何一个标准。 EBS 是风暴的整体切变，在其深度而非标准垂直范围内归一化。

考虑打开功能请求或拉取请求以将此功能添加到 MetPy！

Answer 2

这是计算有效风暴相对螺旋度和有效体积剪切的有效层的快速尝试。 这使用 MetPy 函数来计算包裹温度曲线以及 CAPE 和 CIN 值。

然后可以使用来自此 function 的 output 设置适当的层的深度和底部，以计算风暴相对螺旋度和体积剪切（然后计算诸如超级单体复合参数之类的东西）。

我已经在 1999 年 5 月 4 日 00 UTC 从 KOUN 上对其进行了测试，它给出了一个合理的答案，但尚未经过广泛测试。

# Effective Shear Algorithm
def effective_layer(p, T, Td, h, height_layer=False):
    '''This function determines the effective inflow layer for a convective sounding.
    
    Input:
      - p: sounding pressure with units
      - T: sounding temperature with units
      - Td: sounding dewpoint temperature with units
      - h: sounding heights with units
      
    Returns:
      - pbot/hbot, ptop/htop: pressure/height of the bottom level, pressure/height of the top level
    '''
    from metpy.calc import parcel_profile, cape_cin
    from metyp.units import units
    import numpy as np
    
    pbot = None
    
    for i in range(p.shape[0]):
        prof = parcel_profile(p[i:], T[i], Td[i])
        sbcape, sbcin = cape_cin(p[i:], T[i:], Td[i:], prof)
        if sbcape >= 100 * units('J/kg') and sbcin > -250 * units('J/kg'):
            pbot = p[i]
            hbot = h[i]
            bot_idx = i
            break
    if not pbot:
        return None, None
        
    for i in range(bot_idx+1, p.shape[0]):
        prof = parcel_profile(p[i:], T[i], Td[i])
        sbcape, sbcin = cape_cin(p[i:], T[i:], Td[i:], prof)
        if sbcape < 100 * units('J/kg') or sbcin < -250 * units('J/kg'):
            ptop = p[i]
            htop = h[i]
            break
            
    if height_layer:
        return hbot, htop
    else:
        return pbot, ptop