gen_pressure_stokes

Converts pressure fields from the spatial domain to spherical harmonic coefficients [2, 24]

Calling Sequence

from model_harmonics.gen_pressure_stokes import gen_pressure_stokes
from gravity_toolkit.associated_legendre import plm_holmes
PLM, dPLM = plm_holmes(LMAX, np.cos(th))
Ylms = gen_pressure_stokes(P, G, R, lon, lat, LMAX=LMAX, PLM=PLM, LOVE=(hl,kl,ll))

Source code

model_harmonics.gen_pressure_stokes(P, G, R, lon, lat, LMAX=60, MMAX=None, PLM=None, LOVE=None)[source]

Converts pressure fields from the spatial domain to spherical harmonic coefficients [2, 24]

Parameters:

P: np.ndarray: Pressure (Pa)
G: np.ndarray: Gravitational acceleration (m/s²)
R: np.ndarray: Radius at point (m)
lon: np.ndarray: longitude array
lat: np.ndarray: latitude array
LMAX: int, default 60: Upper bound of Spherical Harmonic Degrees
MMAX: int or NoneType, default None: Upper bound of Spherical Harmonic Orders
PLM: np.ndarray or NoneType, default None: Legendre polynomials
LOVE: tuple or NoneType, default None: Load Love numbers up to degree LMAX (hl, kl, ll)

Returns:

clm: np.ndarray: fully-normalized cosine spherical harmonic coefficients
slm: np.ndarray: fully-normalized sine spherical harmonic coefficients
l: np.ndarray: spherical harmonic degree to LMAX
m: np.ndarray: spherical harmonic order to MMAX