From e3d1cfbba1bff3d10832e24142de7f19d1e8534f Mon Sep 17 00:00:00 2001
From: Carlisle Wishard <cwishard@purdue.edu>
Date: Tue, 1 Jun 2021 17:21:49 -0400
Subject: [PATCH] symba_frag_pos set Escale to be ke_family and adjust all
 scale factors

---
 src/symba/symba_frag_pos.f90 | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
index cb2a1525a..b2641f90e 100644
--- a/src/symba/symba_frag_pos.f90
+++ b/src/symba/symba_frag_pos.f90
@@ -156,10 +156,10 @@ subroutine set_scale_factors()
       ! Set scale factors
       mscale = mtot !! Because of the implied G, mass is actually G*mass with units of distance**3 / time**2
       rscale = sum(radius(:))
-      vscale = sqrt(mscale / rscale)
+      Escale = ke_family !mscale * vscale**2
+      vscale = sqrt(2 * Escale / mscale) !sqrt(mscale / rscale)
       tscale = rscale / vscale
       Lscale = mscale * rscale * vscale
-      Escale = mscale * vscale**2
 
       xcom(:) = xcom(:) / rscale
       vcom(:) = vcom(:) / vscale
@@ -173,7 +173,15 @@ subroutine set_scale_factors()
 
       m_frag = m_frag / mscale
       rad_frag = rad_frag / rscale
+      rot_frag = rot_frag * tscale
       Qloss = Qloss / Escale
+      ke_family = ke_family / Escale
+      Etot_before = Etot_before / Escale
+      pe_before = pe_before / Escale
+      ke_spin_before = ke_spin_before / Escale
+      ke_orbit_before = ke_orbit_before / Escale
+      Ltot_before = Ltot_before / Lscale
+      Lmag_before = Lmag_before / Lscale
 
       return
    end subroutine set_scale_factors