symba_frag_pos scaled and unscaled after energies

MintonGroup · Jun 2, 2021 · 14b59e9 · 14b59e9
1 parent 1fb1c06
commit 14b59e9
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
@@ -182,6 +182,13 @@ subroutine set_scale_factors()
       ke_orbit_before = ke_orbit_before / Escale
       Ltot_before = Ltot_before / Lscale
       Lmag_before = Lmag_before / Lscale
+      Etot_after = Etot_after / Escale
+      pe_after = pe_after / Escale
+      ke_spin_after = ke_spin_after / Escale
+      ke_orbit_after = ke_orbit_after / Escale
+      Ltot_after = Ltot_after / Lscale
+      Lmag_after = Lmag_after / Lscale
+
 
       return
    end subroutine set_scale_factors
@@ -223,6 +230,12 @@ subroutine restore_scale_factors()
       ke_orbit_before = ke_orbit_before * Escale
       Ltot_before = Ltot_before * Lscale
       Lmag_before = Lmag_before * Lscale 
+      Etot_after = Etot_after * Escale
+      pe_after = pe_after * Escale
+      ke_spin_after = ke_spin_after * Escale
+      ke_orbit_after = ke_orbit_after * Escale
+      Ltot_after = Ltot_after * Lscale
+      Lmag_after = Lmag_after * Lscale 
 
       mscale = 1.0_DP
       rscale = 1.0_DP