Added mass loss exception commit

MintonGroup · Aug 30, 2021 · a2b8b4a · a2b8b4a
1 parent 1312cb7
commit a2b8b4a
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/io/io.f90 b/src/io/io.f90
@@ -67,6 +67,10 @@ module subroutine io_conservation_report(self, param, lterminal)
             Ecoll_error = param%Ecollisions / abs(param%Eorbit_orig)
             Etotal_error = (Eorbit_now - param%Ecollisions - param%Eorbit_orig - param%Euntracked) / abs(param%Eorbit_orig)
             Merror = (GMtot_now - param%GMtot_orig) / param%GMtot_orig
+            if (Merror < -10 * epsilon(Merror)) then
+               write(*,*) 'Mass loss! Halting!'
+               call util_exit(FAILURE)
+            end if
             write(*, EGYTERMFMT) Lerror, Ecoll_error, Etotal_error, Merror
          end if
       end associate