MandS#

exo_skryer.vert_Tp.MandS(p_lev: Array, params: Dict[str, Array]) → Tuple[Array, Array][source]#

Generate a Madhusudhan & Seager (2009) three-region T-P profile.

This profile divides the atmosphere into three regions defined by pressure boundaries and slope parameters.

Parameters:
p_levndarray, shape (nlev,)

Pressure at atmospheric levels.

paramsdict[str, ndarray]

Parameter dictionary containing:

  • a1, a2float

    Shape/slope parameters controlling the inversion strength.

  • log_10_P1, log_10_P2, log_10_P3float

    Transition pressures in log₁₀(bar). Profile is computed in log₁₀ space throughout, matching the Madhusudhan & Seager (2009) convention.

  • T_reffloat

    Reference temperature at the top of the atmosphere in Kelvin.

Returns:
T_levndarray, shape (nlev,)

Temperature at levels in Kelvin.

T_layndarray, shape (nlev-1,)

Temperature at layer midpoints in Kelvin.