From a9f58877e886533ab924fe403a58df4c073d6704 Mon Sep 17 00:00:00 2001 From: David Minton Date: Mon, 7 Nov 2022 16:35:45 -0500 Subject: [PATCH] Removed the longdouble conversion from constant numbers as it can screw with xarray --- python/swiftest/swiftest/constants.py | 21 ++++++++++++--------- 1 file changed, 12 insertions(+), 9 deletions(-) diff --git a/python/swiftest/swiftest/constants.py b/python/swiftest/swiftest/constants.py index 5bd5bc7a4..fe3253c86 100644 --- a/python/swiftest/swiftest/constants.py +++ b/python/swiftest/swiftest/constants.py @@ -13,15 +13,18 @@ import astropy.constants as const # Constants in SI units -GC = np.longdouble(const.G.value) -AU2M = np.longdouble(const.au.value) -GMSunSI = np.longdouble(const.GM_sun.value) -MSun = np.longdouble(const.M_sun.value) -RSun = np.longdouble(const.R_sun.value) +GC = const.G.value +AU2M = const.au.value +GMSun = const.GM_sun.value +MSun = const.M_sun.value +RSun = const.R_sun.value +MEarth = const.M_earth.value +REarth = const.R_earth.value +GMEarth = const.GM_earth.value JD2S = 86400 -YR2S = np.longdouble(365.25 * JD2S) -einsteinC = np.longdouble(299792458.0) +YR2S = 365.25 * JD2S +einsteinC = 299792458.0 # Solar oblatenes values: From Mecheri et al. (2004), using Corbard (b) 2002 values (Table II) -J2Sun = np.longdouble(2.198e-7) -J4Sun = np.longdouble(-4.805e-9) +J2Sun = 2.198e-7 +J4Sun = -4.805e-9