diff --git a/src/io/io_conservation_report.f90 b/src/io/io_conservation_report.f90 index d240eb90a..61816acc9 100644 --- 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