我正在尝试将风场 (u,v) 分解为无旋和非发散场 (uchi, vchi, upsi, vpsi)。该字段是全局的，格式为（时间、级别、纬度、经度），并且当我使用 pyspharm python 模块中的 HelmHoltz 分解和一个没有任何缺失值的字段时有效 - 这是公开的情况可用数据（已插值到表面的再分析数据）。Pyspharm 使用球谐函数计算来计算。
但是，当我尝试使用从我自己的个人模型中获得的数据时，问题就出现了，由于山脉或高地，我无法在表面上进行插值。一旦给定级别上有一个 NaN，整个级别都会返回一个 NaN。
我最初考虑水平插值以平滑这些 NaN，但这些区域通常太大，所以我认为也许有一种不同的方法来进行亥姆霍兹分解，它只会在初始 NaN 附近返回 NaN（而不是整个级别），比如有限差分？
有谁知道我如何或是否可以为 HelmHoltz 执行有限差分，如果可以，我可以在哪里获得代码？THX