Commit 4a15cef5 authored by Julia Stadler's avatar Julia Stadler
Browse files

thermodynamics header and input module updated

parent 862838a1
CHANGES TO CLASS v2.7
-----------------------------------------
* this file lists all changes on order to incorporate DM-neutrino scattering
into CLASS
1. 'thermodynamics.h'
----------------------------------------
* in struct thermo
line 154: double u_nuDM
-> add the coupling strength as a new parameter, the mass is normalized to 100 GeV
line 156: short has_coupling_nuDM
-> flag to indicate weather the new coupling should be taken into account
line 178: int index_th_dmu_nuDM
-> index for the DM-neutrino scattering rate
2. 'input.c'
----------------------------------------
* in function input_read_parameters
......@@ -1293,6 +1293,11 @@ int input_read_parameters(
}
}
}
/** parameters for DM-nu scattering */
class_read_double("u_nuDM", pth->u_nuDM);
if (pth->u_nuDM>0.)
pth->has_coupling_nuDM = _TRUE_;
/** (c) define which perturbations and sources should be computed, and down to which scale */
......@@ -3018,6 +3023,9 @@ int input_default_params(
pth->compute_damping_scale = _FALSE_;
pth->u_nuDM=0.;
pth->has_coupling_nuDM=_FALSE_;
/** - perturbation structure */
ppt->has_cl_cmb_temperature = _FALSE_;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment