From 02638304b53c02a83df0373f7c095db8667c0108 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <rodierq@sxphynh>
Date: Mon, 24 Jan 2022 11:15:36 +0100
Subject: [PATCH] M. Mandement 24/01/2022:bugfix : init of theta_l in mixtures
 was too low (0.1K) ==> 300K This resolves the bug trigged in
 th_r_from_thl_rt_1d.f90 at the LOG(ZT) when some HALO points has T values <
 0. TODO: the computation in compute_entr_detr should be done only on updraft
 points

---
 src/MNH/compute_entr_detr.f90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/MNH/compute_entr_detr.f90 b/src/MNH/compute_entr_detr.f90
index 80a9d68db..004f706a1 100644
--- a/src/MNH/compute_entr_detr.f90
+++ b/src/MNH/compute_entr_detr.f90
@@ -122,6 +122,7 @@ END MODULE MODI_COMPUTE_ENTR_DETR
 !!      R.Honnert Oct 2016 : Update with AROME
 !  P. Wautelet 08/02/2019: bugfix: compute ZEPSI_CLOUD only once and only when it is needed
 !  P. Wautelet 10/02/2021: bugfix: initialized PPART_DRY everywhere
+!  M. Mandement 24/01/2022:bugfix: init of theta_l in mixtures was too low (0.1K) 
 !! --------------------------------------------------------------------------
 !
 !*      0. DECLARATIONS
@@ -232,7 +233,7 @@ INTEGER :: JI,JLOOP
   ZFRAC_ICE(:)=PFRAC_ICE(:) ! to not modify fraction of ice
  
   ZPRE(:)=PPRE_MINUS_HALF(:)
-  ZMIXTHL(:)=0.1
+  ZMIXTHL(:)=300.0
   ZMIXRT(:)=0.1
 
   !Initialize PPART_DRY everywhere to prevent access to non-initialized values
-- 
GitLab