diff --git a/src/pw/ps_implicit_methods.F b/src/pw/ps_implicit_methods.F
index f5122280a1..bd03bbc891 100644
--- a/src/pw/ps_implicit_methods.F
+++ b/src/pw/ps_implicit_methods.F
@@ -1150,6 +1150,11 @@ CONTAINS
 
          ALLOCATE (ps_implicit_env%idx_1dto3d(data_size))
          l = 1
+         ! Suppress OpenMP (at least the Intel compiler has an issue here)
+         ! An automatic OpenMP parallelization of this loop might be tricky
+         ! because of the l incrementation
+!$OMP PARALLEL IF(.FALSE.)
+!$OMP DO
          DO k = lb3, ub3
             DO j = lb2, ub2
                DO i = lb1, ub1
@@ -1160,6 +1165,8 @@ CONTAINS
                END DO
             END DO
          END DO
+!$OMP END DO
+!$OMP END PARALLEL
 
          n_contacts = SIZE(ps_implicit_env%contacts)
          n_tiles_tot = 0
@@ -1281,7 +1288,7 @@ CONTAINS
          l = 1
          ! Suppress OpenMP (at least the Intel compiler has an issue here)
          ! An automatic OpenMP parallelization of this loop might be tricky
-         ! because of the incrementation of l
+         ! because of the l incrementation
 !$OMP PARALLEL IF(.FALSE.)
 !$OMP DO
          DO k = lb3, ub3