Removed bad array index in potential energy summation.

MintonGroup · Jan 27, 2022 · f4fb1ae · f4fb1ae
1 parent 50fed7a
commit f4fb1ae
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/util/util_get_energy_momentum.f90 b/src/util/util_get_energy_momentum.f90
@@ -199,7 +199,7 @@ subroutine util_get_energy_potential_triangular(npl, lmask, GMcb, Gmass, mass, x
          do concurrent(j = i+1:npl, lmask(i) .and. lmask(j))
             pepl(j) = - (Gmass(i) * mass(j)) / norm2(xb(:, i) - xb(:, j))
          end do
-         pe = pe + sum(pepl(i+1:npl), lmask(i) .and. lmask(j))
+         pe = pe + sum(pepl(i+1:npl), lmask(i))
       end do
       !$omp end parallel do
       pe = pe + sum(pecb(1:npl), lmask(1:npl))