From e3d1cf31e21594db3087f2097139c46ddd9efa44 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 18 Jun 2021 14:37:07 -0400 Subject: [PATCH] Corrected bad angular momentum calculation for colliding bodies with pre-existing spin. I accidentally put in the kinetic energy formula instead of angular momentum. Oopsies. --- src/symba/symba_frag_pos.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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