Added the central body's contribution to orbital pre-collision angula…

…r momentum to the post-collision spin.

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)