Fixed bug resulting from accidentally removing the status flag settin…

…g for encounters
MintonGroup · Sep 17, 2021 · 9bec149 · 9bec149
1 parent 8ae5215
commit 9bec149
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 1 deletion.
diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90
@@ -291,7 +291,6 @@ module function symba_collision_check_encounter(self, system, param, t, dt, irec
       class is (symba_pl)
          select type(tp => system%tp)
          class is (symba_tp)
-            nenc = self%nenc
             nenc = self%nenc
             allocate(lmask(nenc))
             lmask(:) = ((self%status(1:nenc) == ACTIVE) .and. (pl%levelg(self%index1(1:nenc)) >= irec))

diff --git a/src/symba/symba_encounter_check.f90 b/src/symba/symba_encounter_check.f90
@@ -188,6 +188,7 @@ module function symba_encounter_check_pl(self, param, system, dt, irec) result(l
                j = plplenc_list%index2(k)
                plplenc_list%id1(k) = pl%id(i)
                plplenc_list%id2(k) = pl%id(j)
+               plplenc_list%status(k) = ACTIVE
                plplenc_list%level(k) = irec
                pl%lencounter(i) = .true.
                pl%lencounter(j) = .true.