swh:1:snp:87b6d113741240534c66635318ee8d2457a8c506
Tip revision: 9d45bdf6c08ad0d5cc831306f7b02640d10471de authored by Jane Lin on 21 October 2016, 09:52:03 UTC
Update mcmc_ages_2.py
Update mcmc_ages_2.py
Tip revision: 9d45bdf
asteroseismic.py
from numpy import sqrt
#reference values; Teff from IAU 2015 Resolution B3
#asteroseismic values from Huber et al. 2011
delta_nu_sol = 135.1 #sigma=0.1
nu_max_sol = 3090.0 #sigma=30.
Teff_sol = 5772.0
def delta_nu_func(M_Mo,Teff,L_Lo):
return delta_nu_sol*sqrt(M_Mo)*pow(Teff/Teff_sol,3)*pow(L_Lo,-0.75)
def nu_max_func(M_Mo,Teff,L_Lo):
return nu_max_sol*M_Mo*pow(Teff/Teff_sol,3.5)/L_Lo
