Added more diagnostic info to the mass conservation failure

MintonGroup · Nov 15, 2021 · 125946f · 125946f
1 parent c2ddf3e
commit 125946f
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/io/io.f90 b/src/io/io.f90
@@ -72,6 +72,10 @@ module subroutine io_conservation_report(self, param, lterminal)
             if (lterminal) write(*, EGYTERMFMT) Lerror, Ecoll_error, Etotal_error, Merror
             if (abs(Merror) > 100 * epsilon(Merror)) then
                write(*,*) "Severe error! Mass not conserved! Halting!"
+               write(*,*) "Merror = ", Merror
+               write(*,*) "GMtot_now : ",GMtot_now
+               write(*,*) "GMtot_orig: ",system%GMtot_orig
+               write(*,*) "Difference: ",GMtot_now - system%GMtot_orig
                call pl%xv2el(cb)
                call self%write_hdr(param%nciu, param)
                call cb%write_frame(param%nciu, param)