diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index 71c38b17f..ac01e0aff 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -506,7 +506,7 @@ subroutine set_fragment_tangential_velocities() call util_crossproduct(v_r_unit(:, i), v_t_unit(:, i), L(:)) A(4:6, i) = m_frag(i) * rmag(i) * L(:) else !For the remining bodies, distribute the angular momentum equally amongs them - v_t_mag(i) = 0.25_DP * L_orb_mag / (m_frag(i) * rmag(i) * nfrag) + v_t_mag(i) = 0.9_DP * L_orb_mag / (m_frag(i) * rmag(i) * nfrag) v_frag(:, i) = v_t_mag(i) * v_t_unit(:, i) L_lin_others(:) = L_lin_others(:) + m_frag(i) * v_frag(:, i) call util_crossproduct(x_frag(:, i), v_frag(:, i), L(:))