diff --git a/src/common/turb/mode_tridiag.F90 b/src/common/turb/mode_tridiag.F90 index e12da996eb8acb466cb8a9ffd0396aeaa12a8866..6baf20df6d83d87e4299b8253b412cf79329ad66 100644 --- a/src/common/turb/mode_tridiag.F90 +++ b/src/common/turb/mode_tridiag.F90 @@ -235,11 +235,11 @@ IF ( PIMPL > 1.E-10 ) THEN ! ELSE ! - !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE) DO JK=IKTB,IKTE + !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE) PVARP(IIB:IIE,IJB:IJE,JK) = ZY(IIB:IIE,IJB:IJE,JK) + !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE) END DO - !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE) ! END IF ! diff --git a/src/common/turb/mode_tridiag_thermo.F90 b/src/common/turb/mode_tridiag_thermo.F90 index e0ed28e1307cd639220a93fc4dbe39dc2a42c412..bd131e170f013a1a904a52d2a61b7ed332bb44f3 100644 --- a/src/common/turb/mode_tridiag_thermo.F90 +++ b/src/common/turb/mode_tridiag_thermo.F90 @@ -196,8 +196,8 @@ ZY(IIB:IIE,IJB:IJE,IKB) = PRHODJ(IIB:IIE,IJB:IJE,IKB)*PVARM(IIB:IIE,IJB:IJE,IKB) - ZRHODJ_DFDDTDZ_O_DZ2(IIB:IIE,IJB:IJE,IKB+D%NKL) * PIMPL * PVARM(IIB:IIE,IJB:IJE,IKB ) !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE) ! -!$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE) DO JK=IKTB+1,IKTE-1 + !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE) ZY(IIB:IIE,IJB:IJE,JK) = PRHODJ(IIB:IIE,IJB:IJE,JK)*PVARM(IIB:IIE,IJB:IJE,JK)/PTSTEP & - ZMZM_RHODJ(IIB:IIE,IJB:IJE,JK+D%NKL) * PF(IIB:IIE,IJB:IJE,JK+D%NKL)/PDZZ(IIB:IIE,IJB:IJE,JK+D%NKL) & + ZMZM_RHODJ(IIB:IIE,IJB:IJE,JK ) * PF(IIB:IIE,IJB:IJE,JK )/PDZZ(IIB:IIE,IJB:IJE,JK ) & @@ -205,8 +205,8 @@ DO JK=IKTB+1,IKTE-1 - ZRHODJ_DFDDTDZ_O_DZ2(IIB:IIE,IJB:IJE,JK+D%NKL) * PIMPL * PVARM(IIB:IIE,IJB:IJE,JK ) & - ZRHODJ_DFDDTDZ_O_DZ2(IIB:IIE,IJB:IJE,JK ) * PIMPL * PVARM(IIB:IIE,IJB:IJE,JK ) & + ZRHODJ_DFDDTDZ_O_DZ2(IIB:IIE,IJB:IJE,JK ) * PIMPL * PVARM(IIB:IIE,IJB:IJE,JK-D%NKL) + !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE) END DO -!$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE) ! !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE) ZY(IIB:IIE,IJB:IJE,IKE) = PRHODJ(IIB:IIE,IJB:IJE,IKE)*PVARM(IIB:IIE,IJB:IJE,IKE)/PTSTEP & diff --git a/src/common/turb/turb.F90 b/src/common/turb/turb.F90 index 42d7bbd87230af1662a647d6948138820dd992bb..a7933fa968c9c770910e2937bad32693a56882de 100644 --- a/src/common/turb/turb.F90 +++ b/src/common/turb/turb.F90 @@ -754,13 +754,13 @@ SELECT CASE (HTURBLEN) ZALPHA=0.5**(-1.5) ! - !$mnh_expand_array(JI=1:D%NIT,JJ=1:D%NJT) DO JK=IKTB,IKTE + !$mnh_expand_array(JI=1:D%NIT,JJ=1:D%NJT) ZLM(:,:,JK) = ( 0.5*(PZZ(:,:,JK)+PZZ(:,:,JK+D%NKL)) - & & PZZ(:,:,D%NKA+JPVEXT_TURB*D%NKL) ) * PDIRCOSZW(:,:) ZLM(:,:,JK) = ZALPHA * ZLM(:,:,JK) * ZL0 / ( ZL0 + ZALPHA*ZLM(:,:,JK) ) + !$mnh_end_expand_array(JI=1:D%NIT,JJ=1:D%NJT) END DO - !$mnh_end_expand_array(JI=1:D%NIT,JJ=1:D%NJT) ! ZLM(:,:,IKTB-1) = ZLM(:,:,IKTB) ZLM(:,:,IKTE+1) = ZLM(:,:,IKTE)