繁体   English   中英

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

How can I get the effective storm helicity and effective shear values for this metpy calculation?

提示:本站收集StackOverFlow近2千万问答,支持中英文搜索,鼠标放在语句上弹窗显示对应的参考中文或英文, 本站还提供   中文繁体   英文版本   中英对照 版本,有任何建议请联系yoyou2525@163.com。

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

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

谢谢!

2 个回复

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

来自Thompson 等人的释义。 2007 年

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

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

这是计算有效风暴相对螺旋度和有效体积剪切的有效层的快速尝试。 这使用 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
1 如何获得有效的文件权限?

我正在尝试获取文件有效权限 。 最好的方法是什么? 我正在尝试使用win32security但GetEffectiveRightsFromAcl(trustee)函数需要PyTRUSTEE参数。 而且我不知道如何正确设置它。 结果,我需要获得与在PowerShell中调用Get- ...

2 如何在pintool中获取数组的内存有效地址?

我尝试打印数组的内存位置(有效地址),发现以下内容: 1-第一种方式:使用(以读取操作为例) IARG_MEMORYREAD_EA 作为分析函数的参数并在此函数内打印此值,(内存读取的有效地址,仅在INS_IsMemoryRead为true且在IPOINT_BEFORE时有效)。 ...

3 ASLR如何有效?

我听说过这个理论。 地址空间位置随机化需要库并将它们加载到虚拟地址空间中的随机位置,以便在黑客在程序中找到漏洞的情况下,他没有预先知道的地址来执行返回到libc的攻击例如,反对。 但在考虑了几秒之后,它作为一种防御措施没有任何意义。 假设我们假设的TargetLib(libc或黑客正在 ...

4 如何使该停止标准有效?

我不知道为什么此停止条件会不断失败,但我认为这与逻辑或(||)有关。 我调试了,它似乎并没有停在我想要的地方,这是在cp(当前位置)在数组的周长上的时候打开的。 这是一种迷宫求解算法。 基本上,迷宫是由X组成的2D阵列墙。 每个带有X的单元格都是一堵墙,“当前位置”无法访问它。 ...

5 无效的有效地址计算

请看下面的程序,错误是无效的有效地址计算,我已经提到那行,请告诉我为什么这里无效的有效地址计算是程序 ...

2010-05-11 08:39:01 2 5716   assembly
6 这种标记解决方案的效果如何?

我正在图像共享网站上工作,并希望为图像实现标记。 我读过问题#20856和#2504150 我对上述问题的方法几乎没有关注。 首先,将图像链接到标签看起来很容易。 然而,通过标签关系获取图像并不容易。 这很简单,因为您必须从一个表中获取图像到标记的关系,然后使用一堆OR语句(每 ...

9 如何使用Javascript或jQuery使此代码更有效?

我想通过&lt;a&gt;标签创建10个体内链接。 在演示中,您可以看到该代码对于3个链接来说太长了,我敢肯定有一种方法可以改进此代码。 因此,我如何知道选择了哪个链接以及如何处理此信息? 演示: http : //jsbin.com/panarupo/1/edit ...

10 如何使此动画代码更有效?

此代码通过ObjectAnimator和AnimationSet在Button单击时切换两个视图(mView1和mView2)的y-position 。 在平移动画时,两个视图的alpha值都将减小并再次增大。 这只是一个设置,可以玩一下。 alpha动画是用XML定义的,而翻译动画是用代 ...

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2022 STACKOOM.COM