From 4e1e44b3e85bd95dafb904f30a17cb9b8d036fe3 Mon Sep 17 00:00:00 2001 From: Carlisle Wishard Date: Thu, 26 Aug 2021 16:26:55 -0400 Subject: [PATCH] changed the names of velocity labels to match horizons database --- python/swiftest/swiftest/init_cond.py | 13 +++++++------ 1 file changed, 7 insertions(+), 6 deletions(-) diff --git a/python/swiftest/swiftest/init_cond.py b/python/swiftest/swiftest/init_cond.py index a159ba9c3..0d7807524 100644 --- a/python/swiftest/swiftest/init_cond.py +++ b/python/swiftest/swiftest/init_cond.py @@ -210,9 +210,9 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): p1.append(pldata[key].vectors()['x'][0] * DCONV) p2.append(pldata[key].vectors()['y'][0] * DCONV) p3.append(pldata[key].vectors()['z'][0] * DCONV) - p4.append(pldata[key].vectors()['vhx'][0] * VCONV) - p5.append(pldata[key].vectors()['vhy'][0] * VCONV) - p6.append(pldata[key].vectors()['vhz'][0] * VCONV) + p4.append(pldata[key].vectors()['vx'][0] * VCONV) + p5.append(pldata[key].vectors()['vy'][0] * VCONV) + p6.append(pldata[key].vectors()['vz'][0] * VCONV) p7.append(pldata[key].elements()['a'][0] * DCONV) p8.append(pldata[key].elements()['e'][0]) p9.append(pldata[key].elements()['incl'][0]) @@ -229,10 +229,11 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds): p7.append(pldata[key].vectors()['x'][0] * DCONV) p8.append(pldata[key].vectors()['y'][0] * DCONV) p9.append(pldata[key].vectors()['z'][0] * DCONV) - p10.append(pldata[key].vectors()['vhx'][0] * VCONV) - p11.append(pldata[key].vectors()['vhy'][0] * VCONV) - p12.append(pldata[key].vectors()['vhz'][0] * VCONV) + p10.append(pldata[key].vectors()['vx'][0] * VCONV) + p11.append(pldata[key].vectors()['vy'][0] * VCONV) + p12.append(pldata[key].vectors()['vz'][0] * VCONV) pvec = np.vstack([p1, p2, p3, p4, p5, p6]) + if ispl: Rpl.append(planetradius[key] * DCONV) GMpl.append(GMcb[0] / MSun_over_Mpl[key])