In Calculator_Cosmology.f90, it says
Code: Select all
real(mcp), parameter :: neutrino_mass_fac= 94.07_mcp !conversion factor for thermal with Neff=3 TCMB-2.7255
!93.014 for 3.046
real(mcp), parameter :: standard_neutrino_neff = 3.046_mcp
and in CosmologyParameterizations.f90 it calculates omnuh2 as follows
Code: Select all
CMB%omnuh2=Params(6)/neutrino_mass_fac*(standard_neutrino_neff/3)**0.75_mcp
where Params(6) = mnu. If I understand correctly, the factor of (3.046/3)^(3/4) is meant to account for the fact that standard N_eff is slightly greater than 3. Multiplying by this factor results in omnuh2 = mnu/93.0025 eV, but shouldn't the denominator be 93.014 as stated in the other file?