Fixed bug in which bodies with more than one encounter have their acc…

…elrations zeroed out between evaluations of each pair
MintonGroup · Aug 11, 2021 · 873f5d1 · 873f5d1
1 parent 00fb68e
commit 873f5d1
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/src/symba/symba_kick.f90 b/src/symba/symba_kick.f90
@@ -133,13 +133,13 @@ module subroutine symba_kick_pltpenc(self, system, dt, irec, sgn)
                else
                   irecl = irec
                end if
+               if (isplpl) then
+                  pl%ah(:,ind1(1:self%nenc)) = 0.0_DP
+                  pl%ah(:,ind2(1:self%nenc)) = 0.0_DP
+               else
+                  tp%ah(:,ind2(1:self%nenc)) = 0.0_DP
+               end if
                do k = 1, self%nenc
-                  if (isplpl) then
-                     pl%ah(:,ind1(k)) = 0.0_DP
-                     pl%ah(:,ind2(k)) = 0.0_DP
-                  else
-                     tp%ah(:,ind2(k)) = 0.0_DP
-                  end if
                   if (isplpl) then
                      lgoodlevel = (pl%levelg(ind1(k)) >= irm1) .and. (pl%levelg(ind2(k)) >= irm1)
                   else