Commit 23168db8 authored by Julia Stadler's avatar Julia Stadler
Browse files

--

parent 80c0d527
......@@ -20,7 +20,7 @@ N_ur = 3.046
---> g-cdm coupling
-------------------
--> u_nuDM(T) is given by u_nuDM_0 * a^n_nuDM
u_nuDM_0 = 1.e-4
u_nuDM_0 = 0.
n_nuDM = 0.
--> for alpha_nuDM give a list of values which apply to increasing moltipoles
......@@ -34,23 +34,21 @@ gauge = newtonian
k_pivot = 0.05
sBBN file = ./bbn/sBBN.dat
--------------------------------------------
---> interacting neutrino initial conditions
--------------------------------------------
has_nuDM_initially = 0
start_small_k_at_dmu_nuDM_over_aH = 0.001
start_large_k_at_aH_over_dmu_nuDM = 0.001
-------------------------------------------------------------------------
--> account for neutrino interactions when setting the initial conditions
-------------------------------------------------------------------------
has_nuDM_initially = 0
-----------
---> OUTPUT
-----------
output = tCl, pCl, mPk, dTk
output = tCl, mPk
extra metric transfer functions = yes
k_output_values = 0.3366, 2.0196, 3.366, 13.464, 20.196, 26.926
k_per_decade_for_pk = 100
P_k_max_h/Mpc = 50
#k_output_values = 0.3366, 2.0196, 3.366, 13.464, 20.196, 26.926, 53.806
k_per_decade_for_pk = 20
P_k_max_h/Mpc = 100
write parameters = yes
root = output/2019-05-24_initial-conditions/u1em4_old_
root = output/2019-05-24_initial-conditions/test_
---------------------------
---> some precision settings
......@@ -68,13 +66,13 @@ l_max_ur = 200
-----------------------
---> breake the silence
-----------------------
input_verbose = 1
background_verbose = 1
thermodynamics_verbose = 1
perturbations_verbose = 1
transfer_verbose = 1
primordial_verbose = 1
spectra_verbose = 1
nonlinear_verbose = 1
lensing_verbose = 1
output_verbose = 1
input_verbose = 5
background_verbose = 5
thermodynamics_verbose = 5
perturbations_verbose = 5
transfer_verbose = 5
primordial_verbose = 5
spectra_verbose = 5
nonlinear_verbose = 5
lensing_verbose = 5
output_verbose = 5
......@@ -493,18 +493,19 @@ int input_init(
if(pth->has_coupling_nuDM==_TRUE_){
printf(" -> Neutrino-Dark Matter coupling enabled with the following parameters:\n");
printf(" u_nuDM_0 = %f\n", pth->u_nuDM_0);
printf(" n_nuDM = %f\n", pth->n_nuDM);
printf(" l_max_ur = %d\n", ppr->l_max_ur);
printf(" ufa_triggr = %f\n", ppr->ur_fluid_trigger_tau_over_tau_k);
printf(" alpha_nuDM = ");
printf(" u_nuDM_0 = %f\n", pth->u_nuDM_0);
printf(" n_nuDM = %f\n", pth->n_nuDM);
printf(" l_max_ur = %d\n", ppr->l_max_ur);
printf(" ufa_triggr = %f\n", ppr->ur_fluid_trigger_tau_over_tau_k);
printf(" alpha_nuDM = ");
for (i=0; i<ppr->l_max_ur; i++)
printf("%f, ", ppt->alpha_nuDM[i]);
printf("\n");
printf("has_nuDM_initially = %d\n", ppr->has_nuDM_initially);
if (i<3 || ppt->alpha_nuDM[i] !=ppt->alpha_nuDM[i])
printf("%f, ", ppt->alpha_nuDM[i]);
printf("...\n");
printf(" has_nuDM_initially = %d\n", ppr->has_nuDM_initially);
if (ppr->has_nuDM_initially == _TRUE_){
printf("small k threshold = %f\n", ppr->start_small_k_at_dmu_nuDM_over_aH);
printf("large k threshold = %f\n", ppr->start_large_k_at_aH_over_dmu_nuDM);
printf(" small k threshold = %f\n", ppr->start_small_k_at_dmu_nuDM_over_aH);
printf(" large k threshold = %f\n", ppr->start_large_k_at_aH_over_dmu_nuDM);
}
}
return _SUCCESS_;
......
......@@ -2252,6 +2252,7 @@ int perturb_solve(
/** check if we have to impose a condition on (aH/kappa_nuDM) **/
account_for_nuDM_interactions == _FALSE_;
if(pth->has_coupling_nuDM==_TRUE_ && ppr->has_nuDM_initially==_TRUE_){
printf('in nuDM ini condition');
if(ppw->pvecthermo[pth->index_th_dmu_nuDM]/
ppw->pvecback[pba->index_bg_a]/
ppw->pvecback[pba->index_bg_H]
......
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