Added non-cb binding energy back to the total orbital energy

MintonGroup · Jan 11, 2023 · 2b33502 · 2b33502
1 parent d1f96be
commit 2b33502
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/swiftest/swiftest_io.f90 b/src/swiftest/swiftest_io.f90
@@ -141,7 +141,7 @@ module subroutine swiftest_io_conservation_report(self, param, lterminal)
          te_now = nbody_system%te
          L_orbit_now(:) = nbody_system%L_orbit(:)
          L_spin_now(:) = nbody_system%L_spin(:)
-         E_orbit_now = ke_orbit_now + pe_now
+         E_orbit_now = ke_orbit_now + pe_now + be_now
          L_total_now(:) = nbody_system%L_total(:) + nbody_system%L_escape(:)
          GMtot_now = nbody_system%GMtot + nbody_system%GMescape