Fixed array range problem that was causing an exception.

MintonGroup · Aug 23, 2021 · c3022d3 · c3022d3
1 parent 321da65
commit c3022d3
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/src/kick/kick.f90 b/src/kick/kick.f90
@@ -41,7 +41,7 @@ module subroutine kick_getacch_int_pl(self)
             if (lmask(i) .and. lmask(j)) call kick_getacch_int_one_pl(rji2, dx, dy, dz, Gmass(i), Gmass(j), ahi(1,i), ahi(2,i), ahi(3,i), ahj(1,j), ahj(2,j), ahj(3,j))
          end do
          !$omp end parallel do
-         ah(:,:) = ah(:,:) + ahi(:,:) + ahj(:,:)
+         ah(:,1:self%nbody) = ah(:,1:self%nbody) + ahi(:,1:self%nbody) + ahj(:,1:self%nbody)
       end associate
 
       return

diff --git a/src/symba/symba_kick.f90 b/src/symba/symba_kick.f90
@@ -42,7 +42,7 @@ module subroutine symba_kick_getacch_int_pl(self)
             if ((rji2 > rlim2) .and. lmask(i) .and. lmask(j)) call kick_getacch_int_one_pl(rji2, dx, dy, dz, Gmass(i), Gmass(j), ahi(1,i), ahi(2,i), ahi(3,i), ahj(1,j), ahj(2,j), ahj(3,j))
          end do
          !$omp end parallel do
-         ah(:,:) = ah(:,:) + ahi(:,:) + ahj(:,:)
+         ah(:,1:self%nbody) = ah(:,1:self%nbody) + ahi(:,1:self%nbody) + ahj(:,1:self%nbody)
       end associate
 
       return
@@ -96,7 +96,7 @@ module subroutine symba_kick_getacch_pl(self, system, param, t, lbeg)
                if (rji2 > rlim2) call kick_getacch_int_one_pl(rji2, dx, dy, dz, Gmass(i), Gmass(j), ahi(1,i), ahi(2,i), ahi(3,i), ahj(1,j), ahj(2,j), ahj(3,j))
             end do
             !$omp end parallel do
-            ah(:,:) = ah(:,:) - ahi(:,:) - ahj(:,:)
+            ah(:,1:self%nbody) = ah(:,1:self%nbody) - ahi(:,1:self%nbody) - ahj(:,1:self%nbody)
          end associate
       end select