diff --git a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb
index d9be0df4d..29dcf43aa 100644
--- a/examples/rmvs_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb
+++ b/examples/rmvs_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb
@@ -81,8 +81,8 @@
     {
      "data": {
       "text/plain": [
-       "[<matplotlib.lines.Line2D at 0x2af2269e12d0>,\n",
-       " <matplotlib.lines.Line2D at 0x2af2269eeb10>]"
+       "[<matplotlib.lines.Line2D at 0x2b47a0e74650>,\n",
+       " <matplotlib.lines.Line2D at 0x2b47a0e71250>]"
       ]
      },
      "execution_count": 6,
diff --git a/src/rmvs/rmvs_step.f90 b/src/rmvs/rmvs_step.f90
index 113b4d02f..74e6958c9 100644
--- a/src/rmvs/rmvs_step.f90
+++ b/src/rmvs/rmvs_step.f90
@@ -442,6 +442,7 @@ subroutine rmvs_make_planetocentric(param, cb, pl, tp)
                      tpenci%cb_heliocentric = cb
                      tpenci%ipleP = i
                      tpenci%lmask(:) = .true.
+                     tpenci%status(:) = ACTIVE
                      ! Grab all the encountering test particles and convert them to a planetocentric frame
                      tpenci%id(:) = pack(tp%id(:), encmask(:)) 
                      do j = 1, NDIM