symba_frag_pos fixed failure criteria

MintonGroup · Jun 4, 2021 · 515b00c · 515b00c
1 parent 14b59e9
commit 515b00c
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
@@ -99,7 +99,7 @@ subroutine symba_frag_pos(param, symba_plA, family, x, v, L_spin, Ip, mass, radi
          if (lmerge) write(*,*) 'Failed to find radial velocities'
          if (.not.lmerge) then
             call calculate_system_energy(linclude_fragments=.true.)
-            if ((abs(dEtot / Etot_before) < Qloss) .or. (dEtot > 0.0_DP)) then
+            if ((abs(dEtot - Qloss) / Qloss > Etol) .or. (dEtot > 0.0_DP)) then
                write(*,*) 'Failed due to high energy error: ', abs(dEtot / Etot_before)
                lmerge = .true.
             else if (abs(dLmag) > Ltol) then