From e2adfd01629c4a26ec6b377715f7efb34355c128 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Wed, 11 Aug 2021 22:52:21 -0400 Subject: [PATCH] Fixed bookkeeping term for collision energy --- 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 499c527cc..3d9be1282 100644 --- a/src/symba/symba_discard.f90 +++ b/src/symba/symba_discard.f90 @@ -320,7 +320,7 @@ module subroutine symba_discard_pl(self, system, param) if (param%lenergy) then call system%get_energy_and_momentum(param) Eorbit_after = system%te - system%Ecollisions = Eorbit_after - Eorbit_before + system%Ecollisions = system%Ecollisions + (Eorbit_after - Eorbit_before) end if end if