Skip to content
Snippets Groups Projects
Commit f33d1c6f 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 8503006d
No related branches found
No related tags found
No related merge requests found
......@@ -1731,7 +1731,7 @@ IF ( TURBN%CTURBDIM /= '1DIM' ) THEN
END IF
! 1D turbulence scheme
!$mnh_expand_array(JIJ=IIJB:IIJE,JK=IKTB:IKTE)
PLM(:,:) = PZZ(:,IKTB+IKL:IKTE+IKL) - PZZ(:,:)
PLM(:,IKTB:IKTE) = PZZ(:,IKTB+IKL:IKTE+IKL) - PZZ(:,IKTB:IKTE)
!$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=IKTB:IKTE)
!$mnh_expand_array(JIJ=IIJB:IIJE)
PLM(:,IKU) = PLM(:,IKE)
......
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