gen_atmosphere_stokes

Converts 3D atmospheric geopotential height and pressure difference fields from the spatial domain to spherical harmonic coefficients [2, 24]

Calling Sequence

from model_harmonics.gen_atmosphere_stokes import gen_atmosphere_stokes
from gravity_toolkit.associated_legendre import plm_holmes
PLM, dPLM = plm_holmes(LMAX, np.cos(th))
Ylms = gen_atmosphere_stokes(GPH, pressure, lon, lat, LMAX=LMAX,
    ELLIPSOID='WGS84', GEOID=geoid, PLM=PLM, LOVE=(hl,kl,ll))

Source code

model_harmonics.gen_atmosphere_stokes(GPH, pressure, lon, lat, LMAX=60, MMAX=None, ELLIPSOID=None, GEOID=None, PLM=None, LOVE=None, METHOD='BC05')[source]

Converts 3D atmospheric geopotential height and pressure difference fields from the spatial domain to spherical harmonic coefficients

Parameters:

GPH: np.ndarray

geopotential heights at model levels

pressure: np.ndarray

pressure differences between model levels

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

ELLIPSOID: str or NoneType, default None

reference ellipsoid name

GEOID: np.ndarray or NoneType, default None

geoid height

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)

METHOD: str, default ‘BC05’

Method of integrating over pressure levels

'BC05': [2]

'SW02': [24]

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