From c48d9593b61937fa42761812c83470eb8364d5e9 Mon Sep 17 00:00:00 2001
From: David Minton <daminton@purdue.edu>
Date: Tue, 18 May 2021 13:40:25 -0400
Subject: [PATCH] Corrected unit conversion factors to account for implied G in
 mass terms

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

diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90
index 9262c8186..762ccf0d1 100644
--- a/src/symba/symba_frag_pos.f90
+++ b/src/symba/symba_frag_pos.f90
@@ -225,9 +225,9 @@ subroutine set_scale_factors()
       vcom(:) = (mass(1) * v(:,1) + mass(2) * v(:,2)) / mtot
 
       ! Set scale factors
-      mscale = mtot
+      mscale = mtot !! Because of the implied G, mass is actually G*mass with units of distance**3 / time**2
       rscale = norm2(x(:,2) - x(:,1))
-      vscale = norm2(v(:,2) - v(:,1))
+      vscale = sqrt(mscale / rscale)
       tscale = rscale / vscale
       Lscale = mscale * rscale * vscale
       Escale = mscale * vscale**2