pvtend.sh_ops.vortdiv_sh

pvtend.sh_ops.vortdiv_sh(u: ndarray, v: ndarray, lat: ndarray, lon: ndarray, R_earth: float = 6371000.0) → Tuple[ndarray, ndarray]

Vorticity ζ and divergence δ from a 2-D wind field.

NH-only inputs are mirrored with the correct parity (u even, v odd).

Parameters:

• u – Zonal wind, shape (nlat, nlon).

• v – Meridional wind, shape (nlat, nlon).

• lat – Latitudes in degrees, ascending.

• lon – Longitudes in degrees.

• R_earth – Earth radius in metres.

Returns:

Tuple (vort, div) each shaped (nlat, nlon).