Commit bb31dffd authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Merge branch 'fix_angmom_bug' into 'master'

Apply angular momentum correction the right way round in BH feedback

See merge request !1076
parents 026e9a4c 8bf0912c
......@@ -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