diff --git a/src/MNH/advecuvw_rk.f90 b/src/MNH/advecuvw_rk.f90 index b8d00b9fbf3dfcbc46f131d5a18a86b0d7eb26ce..403eb32770b116afca596af5dd44d35e9fbd9976 100644 --- a/src/MNH/advecuvw_rk.f90 +++ b/src/MNH/advecuvw_rk.f90 @@ -416,13 +416,15 @@ IF ( GIBM ) THEN !$acc end kernels END IF ! -!$acc kernels present(PRUS_ADV,PRVS_ADV,PRWS_ADV) present(ZUT,ZVT,ZWT) present(PU,PV,PW) IF (GIBM .AND. CIBM_ADV=='FREEZE') THEN +!$acc kernels WHERE (XIBM_LS(:,:,:,2).GT.-ZIBM_EPSI) ZIBM(:,:,:,1) = 0. WHERE (XIBM_LS(:,:,:,3).GT.-ZIBM_EPSI) ZIBM(:,:,:,2) = 0. WHERE (XIBM_LS(:,:,:,4).GT.-ZIBM_EPSI) ZIBM(:,:,:,3) = 0. +!$acc end kernels ENDIF ! +!$acc kernels present(PRUS_ADV,PRVS_ADV,PRWS_ADV) present(ZUT,ZVT,ZWT) present(PU,PV,PW) PRUS_ADV = 0. PRVS_ADV = 0. PRWS_ADV = 0. @@ -492,13 +494,13 @@ RKLOOP: DO JS = 1, ISPL !* 4. Advection with WENO ! -------------------------- ! -!$acc kernels IF (GIBM .AND. CIBM_ADV=='LOWORD') THEN +!$acc kernels ZIBM(:,:,:,1)=ZRUS(:,:,:,JS) ZIBM(:,:,:,2)=ZRVS(:,:,:,JS) ZIBM(:,:,:,3)=ZRWS(:,:,:,JS) - ENDIF !$acc end kernels + ENDIF ! IF (HUVW_ADV_SCHEME=='WENO_K') THEN CALL ADVECUVW_WENO_K (HLBCX, HLBCY, KWENO_ORDER, ZUT, ZVT, ZWT, &