From c3022d3121b40fae4d7b77d851257ec388035dcd Mon Sep 17 00:00:00 2001 From: David A Minton Date: Mon, 23 Aug 2021 07:04:48 -0400 Subject: [PATCH] Fixed array range problem that was causing an exception. --- src/kick/kick.f90 | 2 +- src/symba/symba_kick.f90 | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/src/kick/kick.f90 b/src/kick/kick.f90 index 49da7b14c..c351d0656 100644 --- 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 index 68ce0086a..d282f9c53 100644 --- 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