diff --git a/src/black_holes/EAGLE/black_holes_iact.h b/src/black_holes/EAGLE/black_holes_iact.h index 3174dea57ade8a9e3a1d9db3aa54cda428e0bc9f..0db85b3924fee6b4a8d0dcc3263c1746788219cc 100644 --- a/src/black_holes/EAGLE/black_holes_iact.h +++ b/src/black_holes/EAGLE/black_holes_iact.h @@ -134,6 +134,9 @@ runner_iact_nonsym_bh_feedback(const float r2, const float *dx, const float hi, hydro_set_physical_internal_energy(pj, xpj, cosmo, u_new); hydro_set_drifted_physical_internal_energy(pj, cosmo, u_new); + + /* Impose maximal viscosity */ + hydro_set_viscosity_alpha_max_feedback(pj); } } diff --git a/src/feedback/EAGLE/feedback_iact.h b/src/feedback/EAGLE/feedback_iact.h index f1096beec2e1beb92fcdd92ed3f47da49d9fe439..4ad145e77d750ce6c6dc8deb33319c409a838ff7 100644 --- a/src/feedback/EAGLE/feedback_iact.h +++ b/src/feedback/EAGLE/feedback_iact.h @@ -300,6 +300,9 @@ runner_iact_nonsym_feedback_apply(const float r2, const float *dx, hydro_set_physical_internal_energy(pj, xpj, cosmo, u_new); hydro_set_drifted_physical_internal_energy(pj, cosmo, u_new); + /* Impose maximal viscosity */ + hydro_set_viscosity_alpha_max_feedback(pj); + /* message( */ /* "We did some heating! id %llu star id %llu probability %.5e " */ /* "random_num %.5e du %.5e du/ini %.5e", */