Correct indexing issue in OMP/Thomas

xcompact3d · Aug 13, 2024 · a04874f · a04874f
1 parent 015b0a1
commit a04874f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/omp/exec_thom.f90 b/src/omp/exec_thom.f90
@@ -131,7 +131,7 @@ subroutine der_univ_thom_per(du, u, tdsops)
             du(i, j, k) = sum(coeffs(1:4) * u(i, jm(1:4), k))
             du(i, j, k) = du(i, j, k) + coeffs(5) * u(i, j, k)
             du(i, j, k) = du(i, j, k) + sum(coeffs(6:9) * u(i, jp(1:4), k))
-            du(i, j, k) = du(i, j, k) - du(i, jm(4), k) * thom_s(j)
+            du(i, j, k) = du(i, j, k) - du(i, jm(1), k) * thom_s(j)
           end do
           !$omp end simd
         end do