thermal_phase_curve

jaxon.thermal.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)[source]

Compute the phase curve evaluated at phases xi.

Warning

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

Parameters

  • xi (array-like) – Orbital phase angle, must be sorted

  • hotspot_offset (float) – Angle of hotspot offset [radians]

  • omega_drag (float) – Dimensionless drag frequency

  • alpha (float) – Dimensionless fluid number

  • C_11 (float) – Spherical harmonic power in the \(m=1\,\ell=1\) mode

  • T_s (float [K]) – Stellar effective temperature

  • a_rs (float) – Semimajor axis in units of stellar radii

  • rp_a (float) – Planet radius normalized by the semimajor axis

  • A_B (float) – Bond albedo

  • theta2d (array-like) – Grid of latitude values evaluated over the surface of the sphere

  • phi2d (array-like) – Grid of longitude values evaluated over the surface of the sphere

  • filt_wavelength (array-like) – Filter transmittance curve wavelengths [m]

  • filt_transmittance (array-like) – Filter transmittance

  • f (float) – Greenhouse parameter (typically ~1/sqrt(2)).

Returns

  • fluxes (tensor-like) – System fluxes as a function of phase angle \(\xi\).

  • T (tensor-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)