From 5d7780a97a971336a83714bd39feb645e2c4ad59 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Tue, 16 Aug 2022 12:45:58 +0200
Subject: [PATCH] Quentin 16/08/2022: bugfix bit-repro packing
turb_ver_dyn_flux
---
src/common/turb/mode_turb_ver_dyn_flux.F90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/common/turb/mode_turb_ver_dyn_flux.F90 b/src/common/turb/mode_turb_ver_dyn_flux.F90
index c4f3cc9ca..680129b12 100644
--- a/src/common/turb/mode_turb_ver_dyn_flux.F90
+++ b/src/common/turb/mode_turb_ver_dyn_flux.F90
@@ -1060,7 +1060,7 @@ IF(HTURBDIM=='3DIM') THEN
ZWORK31D(IIJB:IIJE ) = - ZFLXZ(IIJB:IIJE,IKB+D%NKL) &
* ( ZWORK1(IIJB:IIJE,IKB+D%NKL) - ZWORK51D(IIJB:IIJE ) &
* PDZY(IIJB:IIJE,IKB+D%NKL) ) &
- / (0.5*(PDYY(IIJB:IIJE,IKB+D%NKL)+PDYY(IIJB:IIJE,IKE)))
+ / (0.5*(PDYY(IIJB:IIJE,IKB+D%NKL)+PDYY(IIJB:IIJE,IKB)))
!$mnh_end_expand_array(JIJ=IIJB:IIJE)
CALL MYF2D_PHY(D,ZWORK31D,ZWORK41D)
ZA(IIJB:IIJE,IKB) = ZWORK41D(IIJB:IIJE)
--
GitLab