From 9bec1497ac7b6ea0090e50939194d524ed688097 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 17 Sep 2021 18:21:24 -0400 Subject: [PATCH] Fixed bug resulting from accidentally removing the status flag setting for encounters --- src/symba/symba_collision.f90 | 1 - src/symba/symba_encounter_check.f90 | 1 + 2 files changed, 1 insertion(+), 1 deletion(-) diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90 index a3990cabd..864eaa723 100644 --- 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 index 743ef677e..3ba6b4adf 100644 --- 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.