From 15f00b51c3ca112d9de047786bb3b67615391e6d Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Wed, 7 Jul 2021 14:30:08 +0200
Subject: [PATCH] Quentin 07/07/2021: bugfix LIMA prognostic supersaturation

---
 src/MNH/prognos_lima.f90 | 1 +
 1 file changed, 1 insertion(+)

diff --git a/src/MNH/prognos_lima.f90 b/src/MNH/prognos_lima.f90
index 967dd166a..648348502 100644
--- a/src/MNH/prognos_lima.f90
+++ b/src/MNH/prognos_lima.f90
@@ -297,6 +297,7 @@ ENDWHERE
 !Calcul de S+dS
 PS0(:)=PS0(:)+((ZA1(:)-(((ZB(:)*(PS0(:)+1.0)+1.0)*ZDZRC(:))/ZRVSAT1(:)))*PTSTEP)
 !
+PS0=MAX(PS0,-0.98)
 !Ajustement tel que rv=(s+1)*rvs
 ZTL(:)=PTT(:)-(PLV(:)/PCPH(:))*PRC(:)
 ZRT(:)=PRC(:)+PRV(:)
-- 
GitLab