symba_frag_pos set Escale to be ke_family and adjust all scale factors

MintonGroup · Jun 1, 2021 · e3d1cfb · e3d1cfb
1 parent 58d65af
commit e3d1cfb
Showing 1 changed file with 10 additions and 2 deletions.
diff --git 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