Skip to content
Snippets Groups Projects
Commit c88d6543 authored by ESCOBAR Juan's avatar ESCOBAR Juan
Browse files

Juan 05/08/2020:tridiag_w.f90, PGI/NVHPC BUG, DOCONCURRENT from 2 -> 3 nestloop

parent 42e1d392
No related branches found
No related tags found
No related merge requests found
......@@ -293,14 +293,14 @@ END DO !CONCURRENT
!$acc end kernels
!
!$acc kernels ! async
DO CONCURRENT ( JI=1:JIU,JJ=1:JJU)
ZY(JI,JJ,IKB+1:IKE-1) = ZMZM_RHODJ(JI,JJ,IKB+1:IKE-1)*PVARM(JI,JJ,IKB+1:IKE-1)/PTSTEP &
- PRHODJ(JI,JJ,IKB+1:IKE-1 ) * PF(JI,JJ,IKB+1:IKE-1 )/PMZF_DZZ(JI,JJ,IKB+1:IKE-1 ) &
+ PRHODJ(JI,JJ,IKB:IKE-2) * PF(JI,JJ,IKB:IKE-2)/PMZF_DZZ(JI,JJ,IKB:IKE-2) &
+ ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,IKB+1:IKE-1 ) * PVARM(JI,JJ,IKB+2:IKE) &
- ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,IKB+1:IKE-1 ) * PVARM(JI,JJ,IKB+1:IKE-1 ) &
- ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,IKB:IKE-2) * PVARM(JI,JJ,IKB+1:IKE-1 ) &
+ ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,IKB:IKE-2) * PVARM(JI,JJ,IKB:IKE-2)
DO CONCURRENT ( JI=1:JIU,JJ=1:JJU,JK=IKB+1:IKE-1)
ZY(JI,JJ,JK) = ZMZM_RHODJ(JI,JJ,JK)*PVARM(JI,JJ,JK)/PTSTEP &
- PRHODJ(JI,JJ,JK ) * PF(JI,JJ,JK )/PMZF_DZZ(JI,JJ,JK ) &
+ PRHODJ(JI,JJ,JK-1) * PF(JI,JJ,JK-1)/PMZF_DZZ(JI,JJ,JK-1) &
+ ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,JK ) * PVARM(JI,JJ,JK+1) &
- ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,JK ) * PVARM(JI,JJ,JK ) &
- ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,JK-1) * PVARM(JI,JJ,JK ) &
+ ZRHODJ_DFDDWDZ_O_DZ2(JI,JJ,JK-1) * PVARM(JI,JJ,JK-1)
END DO !CONCURRENT
!$acc end 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