From f33d1c6fc813ed7560452956521168ed1c8758ee Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 24 Feb 2023 16:28:27 +0100
Subject: [PATCH] Quentin 24/02/2023: bugfix, DEAR computation of PLM must be
 done only on IKTB:IKTE here

---
 src/PHYEX/turb/turb.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/PHYEX/turb/turb.f90 b/src/PHYEX/turb/turb.f90
index e3f302639..6da972021 100644
--- a/src/PHYEX/turb/turb.f90
+++ b/src/PHYEX/turb/turb.f90
@@ -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)
-- 
GitLab