pvtend.isentropic_interpolation_pressure

pvtend.isentropic_interpolation_pressure(theta_levels: ndarray, pressure_hpa: ndarray, temperature_k: ndarray, *, vertical_axis: int = 0, max_iters: int = 50, eps: float = 1e-06, bottom_up_search: bool = True) → ndarray

Compute pressure on isentropic surfaces (MetPy algorithm).

Parameters

theta_levels(n_theta,)

Target potential-temperature levels [K].

pressure_hpa(nlev,) or same shape as temperature_k

Pressure [hPa].

temperature_k : (nlev, …) — temperature [K] on isobaric levels. vertical_axis : int

Axis of temperature_k that is the vertical.

max_iters, eps : Newton-Raphson controls. bottom_up_search : bool

Search direction for bounding indices.

Returns

isen_prs(n_theta, …) — pressure [hPa] on each θ surface.

NaN where the θ level is outside the data range.