Fixed units in the new edge case GMass calc

MintonGroup · Feb 25, 2024 · c486412 · c486412
1 parent 86122b1
commit c486412
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/swiftest/init_cond.py b/swiftest/init_cond.py
@@ -66,7 +66,7 @@ def get_Gmass(raw_response):
                     M = M.split('=')[1].strip().split(' ')[0].strip()
                     M = float(M)  * mult * unit_conv
                     try:
-                        return M * swiftest.GC
+                        return M * swiftest.GC * 1e-9 # Return units of km**3 / s**2 for consistency
                     except:
                         return None
             return None
@@ -537,4 +537,4 @@ def vec2xr(param: Dict, **kwargs: Any):
                             )
         ds = xr.combine_by_coords([ds, clm_xr])
 
-    return ds
+    return ds