爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
12
返回列表 发新帖
楼主: jl2587t

[经验总结] 基于metpy库画T-LnP图(国内主流形式,并非斜温)

[复制链接]
发表于 2024-12-23 21:32:29 | 显示全部楼层
大佬,采用帖子里方法基于ERA5再分析数据尝试画了格点数据探空图,发现和micaps里差别有点大,metpy算出来的cape更是完全不对,请问是什么原因?
使用的露点温度计算公式如下:td0 = mpcalc.dewpoint_from_specific_humidity(p,t0,q0)
使用的CAPE计算公式如下:#计算假相当位温:cape, cin = metpy.calc.cape_cin(p1, t1, td1, parcel_prof,
                                which_lfc='bottom', which_el='top')
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-12-24 07:35:10 来自手机 | 显示全部楼层
球求大人的球球 发表于 2024-12-23 21:32
大佬,采用帖子里方法基于ERA5再分析数据尝试画了格点数据探空图,发现和micaps里差别有点大,metpy算出来 ...

micaps实况数据 和 ERA5格点数据,抬升起始点高度/数值不同、中间层密度不同、中间层数值不同……甚至连空间经纬度位置都不一样,探空曲线能相似就很好了,其他物理量参数一定不一样(除了0J/kg),但是总体范围、变化趋势是相似的。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-12-24 09:14:43 | 显示全部楼层
edwardli 发表于 2024-12-24 07:35
micaps实况数据 和 ERA5格点数据,抬升起始点高度/数值不同、中间层密度不同、中间层数值不同……甚至连 ...

谢谢大佬还有一个小问题。请问micaps里计算cape的算法和metpy.calc.cape_cin()是不是也不一样?因为micaps实况资料metpy算出来的cape和系统自带的还是不一致。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-12-24 09:33:09 来自手机 | 显示全部楼层
本帖最后由 edwardli 于 2024-12-25 09:29 编辑
球求大人的球球 发表于 2024-12-24 09:14
谢谢大佬还有一个小问题。请问micaps里计算cape的算法和metpy.calc.cape_cin()是不是也不一样? ...

当然不一样,你看到的所有工具,用一个探空计算得到的数据都不一样,20个工具会得到20个值。积分步长、虚温订正、相当温度订正、混合气块、LFC精度、LCL精度、EL精度…都会影响CAPE/CIN值。
甚至连MetPy自己升级前后,同一个探空数据的CAPE/cin都不一样,升级后就是虚温订正参与计算。
如上:
  • 只要趋势一致即可
  • 定性用没问题
  • 75%半定量使用(PS.  CAPE3000和3500甚至4500,在释用层面没有本质区别)。

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-12-24 10:08:19 | 显示全部楼层
非常感谢楼主分享经验
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表