From 125946f35741fb0655ec112f56a37c46b91058af Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 15 Nov 2021 15:26:58 -0500 Subject: [PATCH] Added more diagnostic info to the mass conservation failure --- src/io/io.f90 | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/io/io.f90 b/src/io/io.f90 index ebdd68531..c657ec058 100644 --- 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)