Commit 8bf0912c authored by Yannick Bahé's avatar Yannick Bahé
Browse files

Apply angular momentum correction the right way round in BH feedback

parent 026e9a4c
......@@ -479,7 +479,7 @@ __attribute__((always_inline)) INLINE static void black_holes_prepare_feedback(
r_times_v_tang * r_times_v_tang * r_times_v_tang;
const double viscous_time = 2. * M_PI * r_times_v_tang_3 /
(1e-6 * alpha_visc * G * G * BH_mass * BH_mass);
const double f_visc = max(Bondi_time / viscous_time, 1.);
const double f_visc = min(Bondi_time / viscous_time, 1.);
/* Limit the Bondi rate by the Bondi viscuous time ratio */
Bondi_rate *= f_visc;
......
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