From 1e032efadca6dd7378ca002d1cb25842928d114f Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Mon, 15 Jan 2024 15:40:29 +0100
Subject: [PATCH] Philippe 15/01/2024: bugfix: wrong loop indices if MNH_BITREP

---
 src/MNH/bl89.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/MNH/bl89.f90 b/src/MNH/bl89.f90
index 6aa37569a..2dc62d15f 100644
--- a/src/MNH/bl89.f90
+++ b/src/MNH/bl89.f90
@@ -319,8 +319,8 @@ ZBL89EXP = LOG(16.)/(4.*LOG(XKARMAN)+LOG(XCED)-3.*LOG(XCMFS))
 #ifdef MNH_COMPILER_NVHPC
 !$acc loop independent collapse(2)
 #endif
-do concurrent( ji = 1 : iiu, jj = 1 : iju )
-  zsqrt_tke(ji, jj) = Br_pow( ztkem(ji,jj), 0.5 )
+do concurrent( ji = 1 : iiu * iju, jk = 1 : ikt )
+  zsqrt_tke(ji, jk) = Br_pow( ztkem(ji,jk), 0.5 )
 end do
 !ZBL89EXP is defined here because (and not in ini_cturb) because XCED is defined in read_exseg (depending on BL89/RM17)
 ZBL89EXP = Br_log( 16. ) / ( 4. * Br_log( XKARMAN )+ Br_log( XCED ) - 3. * Br_log( XCMFS) )
-- 
GitLab