Removed the longdouble conversion from constant numbers as it can screw with xarray

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