diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index 2b07d094a..4b9e3909b 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -554,7 +554,7 @@ subroutine set_fragment_tangential_velocities(lerr) ! Start the first two bodies with the same rotation as the original two impactors, then distribute the remaining angular momentum among the rest do i = 1, 2 rot_frag(:, i) = rot(:, i) - L_frag_spin(:) = L_frag_spin(:) + m_frag(i) * rad_frag(i)**2 * Ip_frag(3, i) * dot_product(rot_frag(:, i), rot_frag(:, i)) + L_frag_spin(:) = L_frag_spin(:) + m_frag(i) * rad_frag(i)**2 * Ip_frag(3, i) * rot_frag(:, i) end do L_frag_orb(:) = L_frag_tot(:) - L_frag_spin(:) L_frag_spin(:) = 0.0_DP