diff --git a/src/MNH/ice4_sedimentation_split.f90 b/src/MNH/ice4_sedimentation_split.f90
index 52d13f1ccfbd381ce191630ba541beabfac78411..bd7f66ddc3ded80335f8972a3f269327312db23a 100644
--- a/src/MNH/ice4_sedimentation_split.f90
+++ b/src/MNH/ice4_sedimentation_split.f90
@@ -748,8 +748,10 @@ DO WHILE (ANY(ZREMAINT>0.))
   ENDIF
 !$acc kernels
   ZMAX_TSTEP(:,:) = ZREMAINT(:,:)
-! acc loop independent private(JI,JJ,JK,JL)
-  !$mnh_do_concurrent( JL = 1 : ISEDIM )
+!PW: does not work with NVHPC 23.5 if !mnh_do_concurrent because it adds a "acc loop independent" and in that case, results are wrong (reduction operation)
+!PW: DO CONCURRENT is probably not correct either (because of the reduction)
+! acc loop independent private(JI,JJ,JK)
+  DO JL = 1, ISEDIM
     JI=I1(JL)
     JJ=I2(JL)
     JK=I3(JL)
@@ -759,7 +761,7 @@ DO WHILE (ANY(ZREMAINT>0.))
                            PRXT(JI, JJ, JK) * PDZZ(JI, JJ, JK) / ZWSED(JI, JJ, JK))
 !acc end atomic
     ENDIF
-  !$mnh_end_do()
+  END DO
   ZREMAINT(:,:) = ZREMAINT(:,:) - ZMAX_TSTEP(:,:)
 !$acc end kernels
 !$acc kernels