我过去曾使用pyShtools将标量函数扩展到球谐函数并合成它们的带限表示。但是，我不太确定如何为球体表面的矢量场实现这一点，其中$u(\theta, \phi)$和$v(\theta, \phi)$分别是同纬度和纵向分量。

我过去也使用spherepack将矢量场分解为极向/环形分量；但不幸的是，我找不到我要找的东西。

我将非常感谢任何帮助我指出正确的方向。