Juan 04/11/2022:tools.f90, BYpass CCE > 14.X , "acc atomic" needed now...

Juan 04/11/2022:tools.f90, BYpass CCE > 14.X , "acc atomic" needed now "copy(ic)" to get corretly the ic value

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