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])