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)