https://github.com/Azeret/galIMF
Revision a0edaeea2d69a6ed92abef579a4acba9f9abdf23 authored by juzikong on 18 September 2019, 07:17:19 UTC, committed by GitHub on 18 September 2019, 07:17:19 UTC
change mass grid parameter 1.05 to 1.01.
1 parent d9b3595
Raw File
Tip revision: a0edaeea2d69a6ed92abef579a4acba9f9abdf23 authored by juzikong on 18 September 2019, 07:17:19 UTC
Update galimf.py
Tip revision: a0edaee
stellar_luminosity.py
# The function here gives the stellar bolometric luminosity relative to the sun [L_sun],
# assuming a simplified form using only the main-sequence luminosity as a function of mass [M_sun].
# See Yan et al. 2019 for details.
# The stellar luminosity should also be a function of Y_for_helium and Z_for_metal, which shall be added later.


def stellar_luminosity_function(mass):
    if mass < 0.23 ** (1 / 1.7):
        lum = 0.23 * mass ** 2.3
    elif mass < 1.96:
        lum = mass ** 4
    elif mass < (32000 / 1.4) ** (1 / 2.5):
        lum = 1.4 * mass ** 3.5
    else:
        lum = 32000 * mass
    return lum
back to top