Fixed mismatch between DataSet labels and JPL/Horizons-retrieved ephermides

MintonGroup · Aug 27, 2021 · 4d53903 · 4d53903
1 parent 2ff7ebe
commit 4d53903
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/python/swiftest/swiftest/init_cond.py b/python/swiftest/swiftest/init_cond.py
@@ -231,7 +231,7 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
             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])
+        pvec = np.vstack([p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12])
 
         if ispl:
             Rpl.append(planetradius[key] * DCONV)
@@ -259,7 +259,6 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
                 pvec = np.vstack([pvec, Ip1, Ip2, Ip3, rotx, roty, rotz])
         else:
             plab = tlab.copy()
-        pvec = np.vstack([pvec, p7, p8, p9, p10, p11, p12])
 
         if idval is None:
             plid = np.array([planetid[key]], dtype=int) + 1