From fa50454b5753f2732640df0f631a00c5138d8f87 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 5 Oct 2021 23:28:24 -0400 Subject: [PATCH] Fixed bug in OMP directive --- src/encounter/encounter.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/encounter/encounter.f90 b/src/encounter/encounter.f90 index e2830445b..02e78cfcd 100644 --- a/src/encounter/encounter.f90 +++ b/src/encounter/encounter.f90 @@ -14,7 +14,7 @@ module subroutine encounter_check_all_flat_plpl(nplplm, k_plpl, x, v, renc, dt, integer(I4B), dimension(:,:), intent(in) :: k_plpl !! List of all pl-pl encounters real(DP), dimension(:,:), intent(in) :: x !! Position vectors of massive bodies real(DP), dimension(:,:), intent(in) :: v !! Velocity vectors of massive bodies - real(DP), dimension(:), intent(in) :: renc !! Hill's radii of massive bodies + real(DP), dimension(:), intent(in) :: renc !! Hill's radii of massive bodies real(DP), intent(in) :: dt !! Step size logical, dimension(:), intent(out) :: lencounter !! Logical array indicating which pair is in an encounter state logical, dimension(:), intent(out) :: loc_lvdotr !! Logical array indicating the sign of v .dot. x for each encounter @@ -25,7 +25,7 @@ module subroutine encounter_check_all_flat_plpl(nplplm, k_plpl, x, v, renc, dt, !$omp parallel do simd default(private) schedule(static)& !$omp shared(nplplm, k_plpl, x, v, renc, dt, lencounter, loc_lvdotr) & - !$omp lastprivate(xr, yr, zr, vxr, vyr, vzr, renc1, renc2) + !$omp lastprivate(xr, yr, zr, vxr, vyr, vzr, renc12) do k = 1_I8B, nplplm i = k_plpl(1, k) j = k_plpl(2, k)