# Licensed under a 3-clause BSD style license - see LICENSE.rst
"""
Astronomical and physics constants in SI units.  See :mod:`astropy.constants`
for a complete listing of constants defined in Astropy.
"""

import numpy as np

from .constant import Constant
from .config import codata

# ASTRONOMICAL CONSTANTS


class IAU2015(Constant):
    default_reference = 'IAU 2015'
    _registry = {}
    _has_incompatible_units = set()


# DISTANCE

# Astronomical Unit (did not change from 2012)
au = IAU2015('au', "Astronomical Unit", 1.49597870700e11, 'm', 0.0,
             "IAU 2012 Resolution B2", system='si')

# Parsec

pc = IAU2015('pc', "Parsec", au.value / np.radians(1. / 3600.), 'm',
             au.uncertainty / np.radians(1. / 3600.),
             "Derived from au + IAU 2015 Resolution B 2 note [4]", system='si')

# Kiloparsec
kpc = IAU2015('kpc', "Kiloparsec",
              1000. * au.value / np.radians(1. / 3600.), 'm',
              1000. * au.uncertainty / np.radians(1. / 3600.),
              "Derived from au + IAU 2015 Resolution B 2 note [4]", system='si')

# Luminosity
L_bol0 = IAU2015('L_bol0', "Luminosity for absolute bolometric magnitude 0",
                 3.0128e28, "W", 0.0, "IAU 2015 Resolution B 2", system='si')


# SOLAR QUANTITIES

# Solar luminosity
L_sun = IAU2015('L_sun', "Nominal solar luminosity", 3.828e26,
                'W', 0.0, "IAU 2015 Resolution B 3", system='si')

# Solar mass parameter
GM_sun = IAU2015('GM_sun', 'Nominal solar mass parameter', 1.3271244e20,
                 'm3 / (s2)', 0.0, "IAU 2015 Resolution B 3", system='si')

# Solar mass (derived from mass parameter and gravitational constant)
M_sun = IAU2015('M_sun', "Solar mass", GM_sun.value / codata.G.value,
                'kg', ((codata.G.uncertainty / codata.G.value) *
                       (GM_sun.value / codata.G.value)),
                f"IAU 2015 Resolution B 3 + {codata.G.reference}",
                system='si')

# Solar radius
R_sun = IAU2015('R_sun', "Nominal solar radius", 6.957e8, 'm', 0.0,
                "IAU 2015 Resolution B 3", system='si')


# OTHER SOLAR SYSTEM QUANTITIES

# Jupiter mass parameter
GM_jup = IAU2015('GM_jup', 'Nominal Jupiter mass parameter', 1.2668653e17,
                 'm3 / (s2)', 0.0, "IAU 2015 Resolution B 3", system='si')

# Jupiter mass (derived from mass parameter and gravitational constant)
M_jup = IAU2015('M_jup', "Jupiter mass", GM_jup.value / codata.G.value,
                'kg', ((codata.G.uncertainty / codata.G.value) *
                       (GM_jup.value / codata.G.value)),
                f"IAU 2015 Resolution B 3 + {codata.G.reference}",
                system='si')

# Jupiter equatorial radius
R_jup = IAU2015('R_jup', "Nominal Jupiter equatorial radius", 7.1492e7,
                'm', 0.0, "IAU 2015 Resolution B 3", system='si')

# Earth mass parameter
GM_earth = IAU2015('GM_earth', 'Nominal Earth mass parameter', 3.986004e14,
                   'm3 / (s2)', 0.0, "IAU 2015 Resolution B 3", system='si')

# Earth mass (derived from mass parameter and gravitational constant)
M_earth = IAU2015('M_earth', "Earth mass",
                  GM_earth.value / codata.G.value,
                  'kg', ((codata.G.uncertainty / codata.G.value) *
                         (GM_earth.value / codata.G.value)),
                  f"IAU 2015 Resolution B 3 + {codata.G.reference}",
                  system='si')

# Earth equatorial radius
R_earth = IAU2015('R_earth', "Nominal Earth equatorial radius", 6.3781e6,
                  'm', 0.0, "IAU 2015 Resolution B 3", system='si')