用 CDO 计算 EOF 会给出特征值计算警告 jacobi 方案没有收敛答案

【问题标题】：Calculating EOF with CDO gives Eigenvalue computation warning jacobi scheme did not converge用 CDO 计算 EOF 会给出特征值计算警告 jacobi 方案没有收敛
【发布时间】：2021-12-17 11:22:58
【问题描述】：

我需要获取每日位势高度场（数据 = geopot.nc）的第一个 EOF 的主要系数（每日）。我只使用寒冷的季节（N，D，J，F，M，A）。我的代码如下：

cdo selmon,11,12,1,2,3,4 geopot.nc geopotw.nc
cdo --reduce_dim - copy geopotw.nc geopotw3.nc #delete 1 dimension of level=1)
cdo sub geopotw3.nc -ydaymean geopotw3.nc deseason.nc
cdo mulcoslat deseason.nc deseasoncos.nc
cdo eof,3 deseasoncos.nc eval.nc eof.nc

我收到一条警告消息：统计模块（警告）：使用单侧 jacobi 方案的特征值计算未正确收敛。 7244721 对列中的 7017371 未达到要求的 1e-12 正交性。

我做错了什么？也许我错过了什么？或者也许还有其他更好的方法来实现我正在寻找的结果？例如在 R 中。感谢您的帮助！

【问题讨论】：

这可能不是问题，但您是否检查过退季步骤是否正常工作？您的文件从年中开始 - 只是为了仔细检查，您可以尝试更“强大”的命令 cdo ydaysub geopotw3.nc -ydaymean geopotw3.nc deseason.nc 看看它是否给出相同的结果？
ps：有助于检查标记遇到的数据以定位问题 ;-) 例如 eof tag="end of file"，而不是您要查找的 EOF ;-) 而 cdo 用于 EMF，我一直在尝试重命名该标签以避免这些错误标签歧义问题，但不幸的是到目前为止无济于事......
感谢您的评论。我尝试了不同的选项，但仍然有相同的警告。我还尝试了命令 cdo sub ifile -timmean ifile anom_file 和计数 eof 警告看起来像：统计模块（警告）：使用单边 jacobi 方案的特征值计算没有正确收敛。 6239278 对列中的 6116698 未达到要求的 1e-12 cdo eof 正交性：已处理 1 个变量 [1283.48s 132MB]
嗨 - 你试过 cdo ydaysub 吗？您在 ncview 中打开了去季节化的文件并检查它看起来还可以吗？
答案有帮助吗？

标签： cdo-climate

【解决方案1】：

I'm not sure if this post will help，您可以通过以下方式使用环境变量更改最大迭代次数：

export MAX_JACOBI_ITER=100

看看这是否有助于解决您的收敛问题？

【讨论】：