python计算副高脊线位置

齐萃妮 · 发表于 2023-8-15 12:18:08

本帖最后由齐萃妮于 2023-8-15 12:21 编辑

我想计算每天的副高脊线位置，是这样定义的：在10°~50°N、90°~180°E范围内，纬向风切变线(即u=0，且∂u/∂y>0)所在纬度位置的平均值。我写了如下代码，但是结果全为nan值，我觉得应该是我的代码那里写错了，请大家帮我指正，谢谢大家。我用的是ERA5的uwnd数据。

齐萃妮 · 发表于 2023-8-15 12:22:44

纬向风切变线那里是 u=0且偏u/偏y>0

Asofhknok · 发表于 2023-8-21 09:31:09

把关键步骤的值打印出来

雨落森林 · 发表于 2023-8-21 12:17:42

这代码看着像是chatgpt写出来的，计算偏导数那步不对，不是这样算的。不过问题不大，主要是，可能数据里面缺少u==0的数据，导致你找不到，就全是nan，建议你条件不要那么苛刻地==0,改成绝对值小于0.1这种的也试试看？

齐萃妮 · 发表于 2023-9-13 10:49:16

雨落森林发表于 2023-8-21 12:17
这代码看着像是chatgpt写出来的，计算偏导数那步不对，不是这样算的。不过问题不大，主要是，可能数据里面 ...

好的谢谢~

什夜0.0 · 发表于 2023-10-5 09:55:29

请问解决了这个问题嘛

weiyya · 发表于 2024-8-27 16:09:40

雨落森林发表于 2023-8-21 12:17
这代码看着像是chatgpt写出来的，计算偏导数那步不对，不是这样算的。不过问题不大，主要是，可能数据里面 ...

请问该怎么算偏导呢

雨落森林 · 发表于 2024-8-28 11:16:58

weiyya 发表于 2024-8-27 16:09
请问该怎么算偏导呢

方法有很多，metpy的函数metpy.calc.gradient可以算，然后是easyclimate也很好用，页面是这个https://shenyulu.github.io/easyc ... ite_difference.html

奥普拉普拉斯 · 发表于 2024-11-10 13:34:55

请问这个问题解决了吗

		自动登录	找回密码
密码			立即注册

[求助] python计算副高脊线位置