CONAN.occultquad_pya.OccultQuadPy.occultquad

CONAN.occultquad_pya.OccultQuadPy.occultquad#

CONAN.occultquad_pya.OccultQuadPy.occultquad(z0, u1, u2, p, nz=None)#

Occultquad routine

Python implementation of code provided by Mandel and Agol (2002). Please cite Mandel & Agol (2002) if you make use of this routine. The original authors are not responsible for the python re-implmentation.

Produces same result as the fortran code but ~30 times slower.

Parameters:
  • z0 (array) – Impact parameter in units of stellar radius

  • u1 (int) – Linear limb-darkening coefficient (gamma_1 in paper)

  • u2 (int) – Quadratic limb-darkening coefficient (gamma_2 in paper)

  • p (int) – Occulting star (or planet) size in units of stellar radius

  • nz (int, optional) – Number of values in z0. Determined here from z0, parameter introduced to maintain API.

Returns:

  • muo1 (array) – Fraction of flux at each z0 (impact parameter) for a limb-darkened source

  • mu0 (array) – Fraction of flux at each z0 (impact parameter) for a uniform source