Philippe 26/06/2023: ice4_sedimentation_split: remove a DO CONCURRENT because...

Philippe 26/06/2023: ice4_sedimentation_split: remove a DO CONCURRENT because of a reduction in loop (wrong reults with NVHPC 23.5)

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