From 77a47abd7b57f7acb3ab4bf153dd84d8a97b69a3 Mon Sep 17 00:00:00 2001 From: Carlisle Wishard Date: Thu, 9 Mar 2023 16:15:41 -0500 Subject: [PATCH] deleted the angular momentum of the central body during calculation of angular momentum error during central body collisions --- 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 5f55578ed..aa1f52296 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(:) + Lcb(:) + cb%dL(:) + cb%dL(:) = Lpl(:) + cb%dL(:) ! 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)