From 1b7e1fb8eb92cf134c5d7ad6a447dc898fe0bb68 Mon Sep 17 00:00:00 2001 From: David Minton Date: Wed, 19 May 2021 17:48:27 -0400 Subject: [PATCH] Fixed mass scale issue that was preventing potential energy to be computed properly --- src/symba/symba_frag_pos.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index ac973f2d5..9fdee5b53 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -434,7 +434,7 @@ subroutine calculate_system_energy(ke_orb, ke_spin, pe, Ltot, linclude_fragments symba_plwksp%helio%swiftest%status(1:npl) = INACTIVE end where - nplm = count(symba_plwksp%helio%swiftest%mass > param%mtiny) + nplm = count(symba_plwksp%helio%swiftest%mass > param%mtiny / mscale) call util_dist_index_plpl(npl_new, nplm, symba_plwksp) call symba_energy_eucl(npl_new, symba_plwksp, param%j2rp2, param%j4rp4, ke_orb, ke_spin, pe, te, Ltot)