Skip to content
Snippets Groups Projects
Commit 91abc302 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Quentin 24/02/2023: bugfix, DEAR computation of PLM must be done only on IKTB:IKTE here

parent cfc0ae97
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/common/turb/turb.F90
+ 1
− 1
View file @ 91abc302
...@@ -1731,7 +1731,7 @@ IF ( TURBN%CTURBDIM /= '1DIM' ) THEN ...@@ -1731,7 +1731,7 @@ IF ( TURBN%CTURBDIM /= '1DIM' ) THEN
END IF END IF
! 1D turbulence scheme ! 1D turbulence scheme
!$mnh_expand_array(JIJ=IIJB:IIJE,JK=IKTB:IKTE) !$mnh_expand_array(JIJ=IIJB:IIJE,JK=IKTB:IKTE)
PLM(IIJB:IIJE,IKTB:IKTE) = PZZ(IIJB:IIJE,IKL+IKTB:IKTE+IKL) - PZZ(IIJB:IIJE,IKTB:IKTE) PLM(IIJB:IIJE,IKTB:IKTE) = PZZ(IIJB:IIJE,IKTB+IKL:IKTE+IKL) - PZZ(IIJB:IIJE,IKTB:IKTE)
!$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=IKTB:IKTE) !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=IKTB:IKTE)
!$mnh_expand_array(JIJ=IIJB:IIJE) !$mnh_expand_array(JIJ=IIJB:IIJE)
PLM(IIJB:IIJE,IKU) = PLM(IIJB:IIJE,IKE) PLM(IIJB:IIJE,IKU) = PLM(IIJB:IIJE,IKE)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme