问题标签 [metpy]

我的目标是使用 mpcalc.frontogenesis 方法计算给定压力水平下的锋面发生。

我已阅读 NOMADS 的相关数据：

计算相关坐标/网格间距值：

用单位标记每个数组：

然后我去计算FGEN：

这会产生以下错误：

温度、u、v数组的形状分别为：

因此，虽然 dx 和 dy 的长度可以理解为 n-1，但我不确定为什么这不起作用。

我同时使用 dry_lapse 和 wet_lapse 来尝试查找数据集的 showalter 索引。我能够让 dry_lapse 成功运行，并希望在 wet_lapse 函数中使用 dry 的温度输出。目前，我的代码如下所示：

其中 p 是数据中的压力水平，tac 是地块剖面温度，以摄氏度为单位。

代码的目的是计算包裹温度，如果它从 p[0] (1008mb) 绝热提升到 p[2] (950mb) 处的 lcl 水平，然后计算包裹温度，如果它随后被提升潮湿从 p[2] 到 [p11]（500mb 级）绝热。我已通过将它们替换为 dry_lapse 中使用的值来检查我在 mot_lapse 中使用的值，但错误无法解决。我不确定错误是否来自函数源代码中的错误，或者我是否忽略了输入中的错误，无论是干失效还是湿失效。

在笔记本上阅读并绘制卫星图像后，我发现很难重新采样和重新投影（投影 = PlateCarre）到区域范围（范围 = -1、19、-1、17）。如果海岸线绘制在图像上，我将不胜感激。如果也可以绘制名为“DNKN”（long = 8.53，lat = 12.05）的电台，我将不胜感激。读取和绘制卫星图像的代码如下。请根据需要进行编辑。谢谢。

卫星图像（'2014_06_16_16Z.jpg'）的链接是 https://drive.google.com/file/d/1fE4V9aHIklGdv7NizQommdWzAvR6DpUN/view?usp=sharing

我按照代码在此站点上绘制了横截面图。（https://unidata.github.io/MetPy/latest/examples/cross_section.html#sphx-glr-examples-cross-section-py）

在示例中，运行以下代码将产生此结果（交叉）。 cross = cross_section(data, start, end).set_coords(('lat', 'lon')) 在此处输入图像描述

但是，与示例不同的是，运行横截面代码后结果值的 x,y 坐标发生了变化。并且经度和纬度是固定的。 在此处输入图像描述 我不明白为什么会出现这样的结果。

我的代码

我的代码必须根据 nc 文件除以等压层来分析，因此我将每个层的值组合在一起。像这样，

所以，我打印NC。结果是

但是，当我运行横截面时，x，y 会更改为开始和结束，并且原始经度和纬度是固定的。

我希望x,y保持不变，并且纬度和经度在示例中的开始和结束范围内。我真的不明白为什么会这样。

感谢您的关注。

作为参考，我附上数据源和初始数据。数据来源：来自 KMA 的 LDAPS 数据

完整的 NC 数据（初始数据）

我正在尝试将基本单位添加到我的数组中以使用 Metpy 进行计算。使用函数调用metpy.units.meters（或.knots 或.degrees）不起作用。

这是我的代码片段：

其中 webSKNT 和 webDRCT 是 pandas 数据帧。

这是我不断收到的错误： AttributeError: module 'metpy.units' has no attribute 'knots'

如果我尝试，我会收到类似的错误：mpunits.meters

我下载了覆盖区域 [30°W-30°E; 的 ERA5 netcdf 数据；0-20°N]。下载的数据包含水平风分量 (u, v) 和比湿度 (q)。我需要通过使用 Metpy 函数advection来计算湿度的水平平流。

我收到这个错误

我正在使用的代码如下所示

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

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

谢谢！

在我的 Conda 环境中，当我运行时pip show MetPy，我的第一行和最后两行是：

当我运行时pip show pint，我的第一行和最后两行是：

这些应该是这些软件包的最新版本。然而，当我尝试做...

...我收到以下错误：

这让我觉得我的包更新遗漏了一些东西。我还缺少另一个依赖项吗？

我有兴趣从过去的探测中计算 CAPE。我做了一些测试计算，由于某种原因，我的 CAPE 值与当前探测提供的值不一致。

例如，查看来自巴西 Boa Vista 的最新探测（随机选择 b/c 温暖到足以在 1 月份获得可观的 CAPE）。数据显示 CAPE 为 189 J/kg，而 metapy 给我的值要低得多。具体来说，如果我使用测深中所有压力水平的所有数据，我得到 72 J/kg。相反，如果我只使用标准压力水平（1000、925 等——如下面的代码所示），我得到 1 J/kg。

我在下面复制了我的代码。有人可以帮我弄清楚如何正确地做到这一点吗？

谢谢，纪尧姆

我从 era-5 下载了 Era5 U 和 V 风组件，我正在使用 xarray 读取 .nc 文件并从数据中选择几个纬度点。在我需要使用metpy.calc函数计算风速和风向之后：

风速工作没有问题，但风向引发错误：

如果我尝试在 U 和 V 中仅使用 1 个元素计算 Wdir，令人惊讶的是它可以工作：