diff --git a/src/MNH/ice4_sedimentation_split.f90 b/src/MNH/ice4_sedimentation_split.f90
index 518ba63fb62cc9c4d9eeb5e568f9d6236e1e2e37..c5fb8ddeafa8102ab60a81db00c4083059751863 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