The default downloaded version of COSMOMC constrains neutrino fraction f_\nu rather than neutrino mass m_\nu. Here we need to specify the limits on f_\nu in the params.ini file. I would rather constrain m_\nu rather than f_\nu. So I made the following modifications.
I just reinterpreted f_\nu as m_\nu in params.ini file. So I specified the limits of m_\nu to this parameter (parameter 6).
In SetForH(Params,CMB,H0) (params_CMB.f90 subroutine) I modified the equation for CMB%omnuh2 as follows.
CMB%omnuh2 = CMB%nufrac/93.1. (\Omega_\nu h^2 = m_\nu/93.1)
Here CMB%nufrac = Params(6) is the newly interpreted neutrino mass (previously it was f_\nu) .
Can somebody please tell me that these modifications are enough so that I will get m_\nu in chains other than f_\nu.
Thanks in advance.
Neutrino mass in COSMOMC.
-
- Posts: 5
- Joined: February 24 2010
- Affiliation: IUCAA
- Contact: