petitRADTRANS.physics

Module Contents

Functions

guillot_global(P, kappa_IR, gamma, grav, T_int, T_equ)

Returns a temperature array, in units of K,

guillot_global_ret(P, delta, gamma, T_int, T_equ)

guillot_modif(P, delta, gamma, T_int, T_equ, ptrans, alpha)

isothermal(P, T)

running_mean(x, N)

make_press_temp(rad_trans_params)

make_press_temp_iso(rad_trans_params)

PT_ret_model(T3, delta, alpha, tint, press, FeH, CO[, ...])

Self-luminous retrieval P-T model.

petitRADTRANS.physics.guillot_global(P, kappa_IR, gamma, grav, T_int, T_equ)

Returns a temperature array, in units of K, of the same dimensions as the pressure P (in bar). For this the temperature model of Guillot (2010) is used (his Equation 29).

Args:
P:

numpy array of floats, containing the input pressure in bars.

kappa_IR (float):

The infrared opacity in units of \(\rm cm^2/s\).

gamma (float):

The ratio between the visual and infrated opacity.

grav (float):

The planetary surface gravity in units of \(\rm cm/s^2\).

T_int (float):

The planetary internal temperature (in units of K).

T_equ (float):

The planetary equilibrium temperature (in units of K).

petitRADTRANS.physics.guillot_global_ret(P, delta, gamma, T_int, T_equ)
petitRADTRANS.physics.guillot_modif(P, delta, gamma, T_int, T_equ, ptrans, alpha)
petitRADTRANS.physics.isothermal(P, T)
petitRADTRANS.physics.running_mean(x, N)
petitRADTRANS.physics.make_press_temp(rad_trans_params)
petitRADTRANS.physics.make_press_temp_iso(rad_trans_params)
petitRADTRANS.physics.PT_ret_model(T3, delta, alpha, tint, press, FeH, CO, conv=True)

Self-luminous retrieval P-T model.

Args:
T3np.array([t1, t2, t3])

temperature points to be added on top radiative Eddington structure (above tau = 0.1). Use spline interpolation, t1 < t2 < t3 < tconnect as prior.

deltafloat

proportionality factor in tau = delta * press_cgs**alpha

alphafloat

power law index in tau = delta * press_cgs**alpha For the tau model: use proximity to kappa_rosseland photosphere as prior.

tintfloat

internal temperature of the Eddington model

pressnp.ndarray

input pressure profile in bar

convbool

enforce convective adiabat yes/no

COfloat

C/O for the nabla_ad interpolation

FeHfloat

metallicity for the nabla_ad interpolation

Returns:
Tretnp.ndarray

The temperature as a function of atmospheric pressure.