thermal_phase_curve

kelp.jax.thermal_phase_curve(xi, hotspot_offset, omega_drag, alpha, C_11, T_s, a_rs, rp_a, A_B, theta2d, phi2d, filt_wavelength, filt_transmittance, f)

Compute the phase curve evaluated at phases xi.

Warning

Assumes xi is sorted, and that theta2d and phi2d are linearly spaced and increasing.

Parameters:

xiarray-like

Orbital phase angle, must be sorted

hotspot_offsetfloat

Angle of hotspot offset [radians]

omega_dragfloat

Dimensionless drag frequency

alphafloat

Dimensionless fluid number

C_11float

Spherical harmonic power in the \(m=1\,\ell=1\) mode

T_sfloat [K]

Stellar effective temperature

a_rsfloat

Semimajor axis in units of stellar radii

rp_afloat

Planet radius normalized by the semimajor axis

A_Bfloat

Bond albedo

theta2darray-like

Grid of latitude values evaluated over the surface of the sphere

phi2darray-like

Grid of longitude values evaluated over the surface of the sphere

filt_wavelengtharray-like

Filter transmittance curve wavelengths [m]

filt_transmittancearray-like

Filter transmittance

ffloat

Greenhouse parameter (typically ~1/sqrt(2)).

Returns:

fluxestensor-like

System fluxes as a function of phase angle \(\xi\).

Ttensor-like

Temperature map

Examples

Users will typically create the theta2d and phi2d grids like so:

>>> # Set resolution of grid points on sphere:
>>> n_phi = 100
>>> n_theta = 10
>>> phi = np.linspace(-2 * np.pi, 2 * np.pi, n_phi, dtype=floatX)
>>> theta = np.linspace(0, np.pi, n_theta, dtype=floatX)
>>> theta2d, phi2d = np.meshgrid(theta, phi)