From 12dd019e280571e3563d15909e047cdde5c82b6c Mon Sep 17 00:00:00 2001 From: David A Minton Date: Wed, 11 Aug 2021 22:46:52 -0400 Subject: [PATCH 1/2] Reduced error tolerance value --- src/fragmentation/fragmentation.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/fragmentation/fragmentation.f90 b/src/fragmentation/fragmentation.f90 index 73cef6240..72fb82d15 100644 --- a/src/fragmentation/fragmentation.f90 +++ b/src/fragmentation/fragmentation.f90 @@ -42,7 +42,7 @@ module subroutine fragmentation_initialize(system, param, family, x, v, L_spin, integer(I4B), parameter :: NFRAG_MIN = 7 !! The minimum allowable number of fragments (set to 6 because that's how many unknowns are needed in the tangential velocity calculation) real(DP) :: r_max_start, r_max_start_old, r_max, f_spin real(DP), parameter :: Ltol = 10 * epsilon(1.0_DP) - real(DP), parameter :: Etol = 1e-9_DP + real(DP), parameter :: Etol = 1e-8_DP integer(I4B), parameter :: MAXTRY = 3000 integer(I4B), parameter :: TANTRY = 3 logical, dimension(size(IEEE_ALL)) :: fpe_halting_modes, fpe_quiet_modes From e2adfd01629c4a26ec6b377715f7efb34355c128 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Wed, 11 Aug 2021 22:52:21 -0400 Subject: [PATCH 2/2] 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