pvtend.helmholtz_decomposition_3d

pvtend.helmholtz_decomposition_3d(u_3d: ndarray, v_3d: ndarray, lat: ndarray, lon: ndarray, R_earth: float = 6371000.0, method: str = 'spectral') → dict[str, ndarray]

Helmholtz decomposition for each vertical level.

Applies helmholtz_decomposition() independently to every level along axis 0.

Parameters:

• u_3d – Zonal wind [m s⁻¹], shape (nlevels, nlat, nlon).

• v_3d – Meridional wind [m s⁻¹], shape (nlevels, nlat, nlon).

• lat – Latitude in degrees (ascending), shape (nlat,).

• lon – Longitude in degrees, shape (nlon,).

• R_earth – Earth radius [m].

• method – "spectral" (default) or "fd".

Returns:

Dictionary with the same keys as helmholtz_decomposition(), each (nlevels, nlat, nlon).