diff --git a/src/MNH/advection_metsv.f90 b/src/MNH/advection_metsv.f90 index 1169bcf63295ce60f47f87a58187d2e1b171a8c6..c7c864704e2172e8ec09784e60cc8faaaf9512c8 100644 --- a/src/MNH/advection_metsv.f90 +++ b/src/MNH/advection_metsv.f90 @@ -844,13 +844,17 @@ DO JSPL=1,KSPLIT !$acc kernels !$acc loop seq DO JR = 1, KRR - ZR(:,:,:,JR) = ZR(:,:,:,JR) + ( ZRRS_PPM(:,:,:,JR) + ZRRS_OTHER(:,:,:,JR) + PRRS_CLD(:,:,:,JR) ) & - * ZTSTEP_PPM / PRHODJ(:,:,:) + DO CONCURRENT ( JI=1:JIU,JJ=1:JJU,JK=1:JKU) + ZR(JI,JJ,JK,JR) = ZR(JI,JJ,JK,JR) + ( ZRRS_PPM(JI,JJ,JK,JR) + ZRRS_OTHER(JI,JJ,JK,JR) + PRRS_CLD(JI,JJ,JK,JR) ) & + * ZTSTEP_PPM / PRHODJ(JI,JJ,JK) + END DO !CONCURRENT END DO !$acc loop seq DO JSV = 1, KSV - ZSV(:,:,:,JSV) = ZSV(:,:,:,JSV) + ( ZRSVS_PPM(:,:,:,JSV) + ZRSVS_OTHER(:,:,:,JSV) + & - PRSVS_CLD(:,:,:,JSV) ) * ZTSTEP_PPM / PRHODJ(:,:,:) + DO CONCURRENT ( JI=1:JIU,JJ=1:JJU,JK=1:JKU) + ZSV(JI,JJ,JK,JSV) = ZSV(JI,JJ,JK,JSV) + ( ZRSVS_PPM(JI,JJ,JK,JSV) + ZRSVS_OTHER(JI,JJ,JK,JSV) + & + PRSVS_CLD(JI,JJ,JK,JSV) ) * ZTSTEP_PPM / PRHODJ(JI,JJ,JK) + END DO !CONCURRENT END DO !$acc end kernels END IF