Fixed bug where status flag was not set properly for test particles i…

…n planetocentric systems
MintonGroup · Aug 2, 2021 · 26ed1c9 · 26ed1c9
1 parent 1e7f96c
commit 26ed1c9
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 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
@@ -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
@@ -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