From 14b59e9f7da51418e9aa0f2dfa861acbc1110d7d Mon Sep 17 00:00:00 2001 From: Carlisle Wishard Date: Wed, 2 Jun 2021 11:30:28 -0400 Subject: [PATCH] symba_frag_pos scaled and unscaled after energies --- src/symba/symba_frag_pos.f90 | 13 +++++++++++++ 1 file changed, 13 insertions(+) diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index 4feb7579d..41b0b9320 100644 --- 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