From 5ed7a59af0e005fd1bfc1842f1b6aaf0aca3fe8e Mon Sep 17 00:00:00 2001
From: Hanno Rein <hanno@hanno-rein.de>
Date: Thu, 13 Jul 2023 18:55:37 -0400
Subject: [PATCH] Set G=1 to avoid rounding errors when using horizons

---
 rebound/simulation.py | 1 +
 1 file changed, 1 insertion(+)

diff --git a/rebound/simulation.py b/rebound/simulation.py
index d58ccb276..755555e1d 100644
--- 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)