Commit b886813a authored by Julia Stadler's avatar Julia Stadler
Browse files

state before turning off the scattering for neutrinos

parent 35ecd1c5
......@@ -45,7 +45,8 @@
* in function input_read_parameters
line 1309: class_read_double("u_nuDM_0", pth->u_nuDM_0)
-> read strength of neutrino-DM cooupling
line 1311: if the coupling strength is greater than 0 set flag has_coupling_nuDM to TRUE
line 1311: if the coupling strength is greater than 0 set flag
has_coupling_nuDM to TRUE
line 1312: if the coupling strength is greater than 0 read n_nuDM
line 2682: if has_coupling_nuDM is TRUE create array for multipole coefficients
line 2684: read in alpha_nuDM from the input file
......@@ -85,6 +86,13 @@
line 953: asign index to pth->index_th_dmu_nuDM
-> only if has_coupling_nuDM is TRUE
* in function 'thermodynamics_output_titles'
line 3723: add a column for DM-neutrino scattering rate
class_store_columntitle(titles, "dmu_nuDM [Mpc^-1]", pth->has_coupling_nuDM);
* in function 'thermodynamics_output_data'
line 3770: write DM-nu scattering rate to file
class_store_double(dataptr,pvecthermo[pth->index_th_dmu_nuDM],pth->has_coupling_nuDM,storeidx);
4. 'perturbations.c'
----------------------------------------
......@@ -94,7 +102,8 @@
* in function 'perturb_derivs'
line 6918: double S_nuDM
-> the following commands are executed only if pth->has_coupling_nuDM is TRUE
line 6983: S_nuDM = 4./3.*pvecback[pba->index_bg_rho_ur]/pvecback[pba->index_bg_rho_cdm];
line 6983: S_nuDM = 4./3.*pvecback[pba->index_bg_rho_ur]
/pvecback[pba->index_bg_rho_cdm];
line 7238: if ppw->approx[ppw->index_ap_rsa] == (int)rsa_off
dy[pv->index_pt_theta_cdm] += S_nuDM*pvecthermo[pth->index_th_dmu_nuDM]
*(y[pv->index_pt_theta_ur]- y[pv->index_pt_theta_cdm]);
......
-----------------------------------------------------
---> COSMOLOGICAL PARAMETERS
---> fix to Planck 2018 values
-----------------------------------------------------
omega_cdm = 0.1186
omega_b = 0.02223
Omega_Lambda = 0.6851
h = 0.678
n_s = 0.9658
A_s = 2.199e-9
omega_cdm = 0.1186
omega_b = 0.02226
h = 0.678
tau_reio = 0.066
n_s = 0.9677
A_s = 2.139e-9
---------------------------------------------------
---> NEUTRINOS
......@@ -22,31 +23,35 @@ N_ur = 3.046
-------------------
---> g-cdm coupling
-------------------
u_nuDM_0 = 1.e-5
u_nuDM_0 = 0.
n_nuDM = 0.
alpha_nuDM = 0., 1.0, 1.5 # the firsts to values are for detla and theta and not used!
# the final value defines all alpha_l for all higher values of l
gauge = newtonian
k_pivot = 0.05
sBBN file = ./bbn/sBBN.dat
alpha_nuDM = 0., 1., 1.5 #the firsts to values are for detla and theta and not used!
-----------
---> OUTPUT
-----------
output = tCl, pCl, mPk
root = output/test_
#root = output/2018-10-24_mixed-damping/lcdm_
format = camb
k_per_decade_for_pk = 80
P_k_max_h/Mpc = 100.
output = mPk
root = output/2018-12-07_mixed-damping-study/lcdm_extra_
#format = class
k_per_decade_for_pk = 100
P_k_max_h/Mpc = 100
write parameters = yes
write background = no
write thermodynamics = no
#k_output_values =
write thermodynamics = yes
k_output_values = 100, 80, 60, 50, 40, 30
-----------------------
---> precision settings
-----------------------
ur_fluid_trigger_tau_over_tau_k = 10.
l_max_ur = 10
ur_fluid_trigger_tau_over_tau_k = 200
l_max_ur = 200
tol_perturb_integration = 1.e-7
#radiation_streaming_trigger_tau_over_tau_k = 100.
#radiation_streaming_trigger_tau_c_over_tau = 100.
-----------------------
---> breake the silence
......
......@@ -7451,7 +7451,6 @@ int perturb_derivs(double tau,
dy[pv->index_pt_shear_ur] =
-3./tau*y[pv->index_pt_shear_ur]
+2./3.*(y[pv->index_pt_theta_ur]+metric_ufa_class);
}
}
}
......
......@@ -3720,6 +3720,7 @@ int thermodynamics_output_titles(struct background * pba,
class_store_columntitle(titles,"kappa' [Mpc^-1]",_TRUE_);
//class_store_columntitle(titles,"kappa''",_TRUE_);
//class_store_columntitle(titles,"kappa'''",_TRUE_);
class_store_columntitle(titles, "dmu_nuDM [Mpc^-1]", pth->has_coupling_nuDM);
class_store_columntitle(titles,"exp(-kappa)",_TRUE_);
class_store_columntitle(titles,"g [Mpc^-1]",_TRUE_);
//class_store_columntitle(titles,"g'",_TRUE_);
......@@ -3766,6 +3767,7 @@ int thermodynamics_output_data(struct background * pba,
class_store_double(dataptr,tau,_TRUE_,storeidx);
class_store_double(dataptr,pvecthermo[pth->index_th_xe],_TRUE_,storeidx);
class_store_double(dataptr,pvecthermo[pth->index_th_dkappa],_TRUE_,storeidx);
class_store_double(dataptr,pvecthermo[pth->index_th_dmu_nuDM],pth->has_coupling_nuDM,storeidx);
//class_store_double(dataptr,pvecthermo[pth->index_th_ddkappa],_TRUE_,storeidx);
//class_store_double(dataptr,pvecthermo[pth->index_th_dddkappa],_TRUE_,storeidx);
class_store_double(dataptr,pvecthermo[pth->index_th_exp_m_kappa],_TRUE_,storeidx);
......
Supports Markdown
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