CE_atmodeller#

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

Compute chemical equilibrium profiles using the atmodeller backend.

Parameters:
p_layndarray, shape (nlay,)

Layer pressures. In the forward model this is typically in dyne cm⁻² and is converted internally to bar for the solver.

T_layndarray, shape (nlay,)

Layer temperatures in Kelvin.

paramsdict[str, ndarray]

Chemical abundance parameters containing:

  • M_to_Hfloat

    Metallicity relative to solar in dex.

  • C_to_Ofloat

    Carbon-to-oxygen ratio (dimensionless).

nlayint

Number of atmospheric layers. Used to broadcast mass constraints and initial guess to shape (nlay, ...).

Returns:
vmr_laydict[str, ndarray]

Dictionary mapping bare species names (e.g. 'H2O', 'CH4') to VMR profiles with shape (nlay,).

Notes

H and He abundances are fixed at solar values and are not scaled by metallicity. All other elements present in the species network (C, N, O, Na, K) are scaled by 10 ** M_to_H; carbon is further set by C_to_O * O.