Fixed mass scale issue that was preventing potential energy to be com…

…puted properly
MintonGroup · May 19, 2021 · 1b7e1fb · 1b7e1fb
1 parent 3323841
commit 1b7e1fb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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)