diff --git a/src/MNH/tools.f90 b/src/MNH/tools.f90
index 846462cbbcc4e3a6b61c4724bb643ba9874501a6..a811b054567e80016073f18b470e4c00af38510a 100644
--- a/src/MNH/tools.f90
+++ b/src/MNH/tools.f90
@@ -188,7 +188,7 @@ ic = 0
 ! different runs of this subroutine BUT final result should be the same
 !Comment the following line + atomic directives to have consistent values for debugging
 !Warning: huge impact on performance
-!$acc parallel loop private(idx) independent
+!$acc parallel loop private(idx) copy(ic) independent
   do ji = 1, size( ltab, 1 )
     if ( ltab(ji ) ) then
 !$acc atomic capture
@@ -198,7 +198,7 @@ ic = 0
       i1(idx) = ji
     end if
   end do
-! acc end kernels
+! acc end parallel
 
 else
    
@@ -259,7 +259,7 @@ ic = 0
 ! different runs of this subroutine BUT final result should be the same
 !Comment the following line + atomic directives to have consistent values for debugging
 !Warning: huge impact on performance
-!$acc parallel loop collapse(2) private(idx) independent
+!$acc parallel loop collapse(2) private(idx) copy(ic) independent
   do jj = 1, size( ltab, 2 )
     do ji = 1, size( ltab, 1 )
       if ( ltab(ji, jj ) ) then
@@ -272,7 +272,7 @@ ic = 0
       end if
     end do
   end do
-! acc end kernels
+! acc end parallel
 
 else
 
@@ -337,7 +337,7 @@ ic = 0
 ! different runs of this subroutine BUT final result should be the same
 !Comment the following line + atomic directives to have consistent values for debugging
 !Warning: huge impact on performance
-!$acc parallel loop collapse(3) private(idx) independent
+!$acc parallel loop collapse(3) private(idx) copy(ic) independent
   do jk = 1, size( ltab, 3 )
     do jj = 1, size( ltab, 2 )
        do ji = 1, size( ltab, 1 )
@@ -353,7 +353,7 @@ ic = 0
       end do
     end do
   end do
-! acc end kernels
+! acc end parellel
 
 else