From 3c9df9f06188fca01db051328ce7693313c8efd7 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Mon, 5 Dec 2022 16:49:26 +0100
Subject: [PATCH] Quentin 05/12/2022: bugfix LES Ocean validated
---
src/common/turb/mode_turb_ver_thermo_flux.F90 | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/src/common/turb/mode_turb_ver_thermo_flux.F90 b/src/common/turb/mode_turb_ver_thermo_flux.F90
index f6b21804e..c88dbec8b 100644
--- a/src/common/turb/mode_turb_ver_thermo_flux.F90
+++ b/src/common/turb/mode_turb_ver_thermo_flux.F90
@@ -594,9 +594,7 @@ ELSE ! atmosp bottom
! atmos top
#ifdef REPRO48
#else
- ZF(IIJB:IIJE,IKE)=0.
- !TODO merge : the following solution must be kept :
- !ZF((IIJB:IIJE,IKE+1)=0.
+ ZF(IIJB:IIJE,IKE+1)=0.
#endif
END IF
!
@@ -989,7 +987,7 @@ IF (KRR /= 0) THEN
! atmos top
#ifdef REPRO48
#else
- ZF(IIJB:IIJE,IKE)=0.
+ ZF(IIJB:IIJE,IKE+1)=0.
#endif
END IF
! Compute the split conservative potential temperature at t+deltat
--
GitLab