From 26ed1c99f3158d986d4ed1bfe2c2534505d9d62a Mon Sep 17 00:00:00 2001
From: David A Minton <daminton@purdue.edu>
Date: Mon, 2 Aug 2021 12:42:43 -0400
Subject: [PATCH] Fixed bug where status flag was not set properly for test
 particles in planetocentric systems

---
 .../1pl_1tp_encounter/swiftest_vs_swifter.ipynb               | 4 ++--
 src/rmvs/rmvs_step.f90                                        | 1 +
 2 files changed, 3 insertions(+), 2 deletions(-)

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