quench_approx#

exo_skryer.vert_chem.quench_approx(p_lay: Array, T_lay: Array, params: Dict[str, Array], nlay: int) → Dict[str, Array][source]#

Compute quenched chemical abundance profiles using the FastChem 5D grid.

Two-step process:

  1. Interpolate the FastChem 5D grid at each layer to obtain chemical equilibrium (CE) abundances for all active species.

  2. Apply the quench approximation only to the species listed in quench_approx.quench_species (YAML config). For each such species, the profile is frozen at the quench level where the chemical timescale exceeds the eddy-mixing timescale; all other species retain their CE values.

Parameters:
p_layndarray, shape (nlay,)

Layer pressures in dyne cm⁻² (converted to bar internally).

T_layndarray, shape (nlay,)

Layer temperatures in Kelvin.

paramsdict[str, ndarray]

Must contain:

  • M_to_H : metallicity relative to solar [dex]

  • C_to_O : carbon-to-oxygen ratio [dimensionless]

  • Kzz : eddy diffusion coefficient [cm² s⁻¹]

  • log_10_g : log₁₀ surface gravity [cm s⁻²]

nlayint

Number of atmospheric layers (unused; kept for API compatibility).

Returns:
vmr_laydict[str, ndarray]

Species VMR profiles with shape (nlay,), plus __mu_lay__ (mean molecular weight from the FastChem grid).