From e5b28e479aceb9f6f0c8149b58c05841e4dff954 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 20 Mar 2023 11:19:23 -0400 Subject: [PATCH] Added the central body's contribution to orbital pre-collision angular momentum to the post-collision spin. --- src/symba/symba_discard.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/symba/symba_discard.f90 b/src/symba/symba_discard.f90 index aa1f52296..a46ce057d 100644 --- a/src/symba/symba_discard.f90 +++ b/src/symba/symba_discard.f90 @@ -186,7 +186,7 @@ subroutine symba_discard_conserve_energy_and_momentum(pl, nbody_system, param, i becb1 = -(3 * cb%Gmass * cb%mass) / (5 * cb%radius) ! Add planet angular momentum to central body accumulator - cb%dL(:) = Lpl(:) + cb%dL(:) + cb%dL(:) = Lpl(:) + cb%dL(:) + Lcb(:) ! Update rotation of central body to by consistent with its angular momentum if (param%lrotation) then drot0(:) = cb%L0(:)/ (cb%Ip(3) * cb%mass * cb%radius**2)