Set G=1 to avoid rounding errors when using horizons

hannorein · Jul 13, 2023 · 5ed7a59 · 5ed7a59
1 parent 7c6a4fb
commit 5ed7a59
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/rebound/simulation.py b/rebound/simulation.py
@@ -1662,6 +1662,7 @@ def add(self, particle=None, **kwargs):
             elif isinstance(particle,str):
                 if self.python_unit_l == 0 or self.python_unit_m == 0 or self.python_unit_t == 0:
                     self.units = ('AU', 'yr2pi', 'Msun')
+                    self.G = 1.0
                 builtindatasets = ["solar system", "outer solar system"]
                 if particle.lower() == "solar system":          # built in test dataset
                     data.add_solar_system(self)