Skip to content
Snippets Groups Projects
Commit 321c3bf2 authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

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)

(cherry picked from commit ad799fd0)
parent ed98eb21
No related branches found
No related tags found
No related merge requests found
...@@ -748,8 +748,10 @@ DO WHILE (ANY(ZREMAINT>0.)) ...@@ -748,8 +748,10 @@ DO WHILE (ANY(ZREMAINT>0.))
ENDIF ENDIF
!$acc kernels !$acc kernels
ZMAX_TSTEP(:,:) = ZREMAINT(:,:) ZMAX_TSTEP(:,:) = ZREMAINT(:,:)
! acc loop independent private(JI,JJ,JK,JL) !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)
!$mnh_do_concurrent( JL = 1 : ISEDIM ) !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) JI=I1(JL)
JJ=I2(JL) JJ=I2(JL)
JK=I3(JL) JK=I3(JL)
...@@ -759,7 +761,7 @@ DO WHILE (ANY(ZREMAINT>0.)) ...@@ -759,7 +761,7 @@ DO WHILE (ANY(ZREMAINT>0.))
PRXT(JI, JJ, JK) * PDZZ(JI, JJ, JK) / ZWSED(JI, JJ, JK)) PRXT(JI, JJ, JK) * PDZZ(JI, JJ, JK) / ZWSED(JI, JJ, JK))
!acc end atomic !acc end atomic
ENDIF ENDIF
!$mnh_end_do() END DO
ZREMAINT(:,:) = ZREMAINT(:,:) - ZMAX_TSTEP(:,:) ZREMAINT(:,:) = ZREMAINT(:,:) - ZMAX_TSTEP(:,:)
!$acc end kernels !$acc end kernels
!$acc kernels !$acc kernels
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment