From 699f1a3993c77b5d691abe60abec7f4ab1981458 Mon Sep 17 00:00:00 2001
From: Gaelle DELAUTIER <gaelle.delautier@meteo.fr>
Date: Tue, 15 May 2018 14:19:38 +0200
Subject: [PATCH]  Q.Rodier 15/5/2018 : new constant XRM17

---
 src/MNH/ini_cturb.f90  | 7 ++++++-
 src/MNH/modd_cturb.f90 | 1 +
 2 files changed, 7 insertions(+), 1 deletion(-)

diff --git a/src/MNH/ini_cturb.f90 b/src/MNH/ini_cturb.f90
index a7a37098f..3e7bc3275 100644
--- a/src/MNH/ini_cturb.f90
+++ b/src/MNH/ini_cturb.f90
@@ -81,10 +81,11 @@ IMPLICIT NONE
 !         1.1 Constant for dissipation of Tke
 !
 !XCED  = 0.70
-XCED  = 0.85
+XCED  = 0.84
 !       Redelsperger-Sommeria (1981) = 0.70
 !       Schmidt-Schumann      (1989) = 0.845
 !       Cheng-Canuto-Howard   (2002) = 0.845
+!       Rodier, Masson, Couvreux, Paci (2017) = 0.34
 !
 !
 !         1.2 Constant for wind pressure-correlations
@@ -162,6 +163,10 @@ XALPSBL = 4.63
 !       Stull 1988                  = 4.75
 !
 !
+!         1.11  Value related to the shear term in mixing length computation
+!
+XRM17 = 0.5  ! Rodier et al 2017
+!
 !
 !         2. Derivated constants
 !            -------------------
diff --git a/src/MNH/modd_cturb.f90 b/src/MNH/modd_cturb.f90
index ed58fd962..6992865f2 100644
--- a/src/MNH/modd_cturb.f90
+++ b/src/MNH/modd_cturb.f90
@@ -68,6 +68,7 @@ REAL,SAVE :: XCDD         ! ct. for the destruction term in the dissipation eq.
 REAL,SAVE :: XCDT         ! ct. for the transport term in the dissipation eq.
 !
 REAL,SAVE :: XTKEMIN      ! mimimum value for the TKE
+REAL,SAVE :: XRM17        ! Rodier et al 2017 constant in shear term for mixing length
 !
 REAL,SAVE :: XLINI        ! initial value for BL mixing length
 REAL,SAVE :: XLINF        ! to prevent division by zero in the BL algorithm
-- 
GitLab