added in high angular momentum failure reporting and exit

MintonGroup · Jun 25, 2021 · 02c4e20 · 02c4e20
1 parent a958a36
commit 02c4e20
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/io/io_conservation_report.f90 b/src/io/io_conservation_report.f90
@@ -74,6 +74,11 @@ module subroutine io_conservation_report(t, symba_plA, npl, j2rp2, j4rp4, param,
                write(*,*) 'dpe      : ',(pe_now - pe_last) / abs(Eorbit_orig)
                write(*,*)
             end if
+            if (Lerror > 1e-10) then
+               write(*,*) 'Something has gone wrong! Angular momentum is too high!'
+               write(*,*) 'Lerror: ', Lerror
+               call util_exit(FAILURE)
+            end if
          end if
          ke_orbit_last = ke_orbit_now
          ke_spin_last = ke_spin_now