add single mass option for neutrinos
@@ -108,7 +108,7 @@ class Cosmology(object):
m_nu (:obj:`float`, optional): Total mass in eV of the massive
neutrinos present. Defaults to 0.
m_nu_type (:obj:`str`, optional): The type of massive neutrinos. Should
- be one of 'inverted', 'normal', 'equal' or 'list'. The default
+ be one of 'inverted', 'normal', 'equal', 'single', or 'list'. The default
of None is the same as 'normal'.
w0 (:obj:`float`, optional): First order term of dark energy equation
of state. Defaults to -1.
@@ -450,6 +450,11 @@ class Cosmology(object):
mnu_list[0] = m_nu[0]/3.
mnu_list[1] = m_nu[0]/3.
mnu_list[2] = m_nu[0]/3.
+ elif (m_nu_type == 'single'):
+ mnu_list = [0]*3
+ mnu_list[0] = m_nu[0]
+ mnu_list[1] = 0.
+ mnu_list[2] = 0.
# Check which of the neutrino species are non-relativistic today
N_nu_mass = 0
