From 8cd8cd58d5aceb5f76227684d342ef51fe1d51d0 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Beno=C3=AEt=20Vi=C3=A9?= <benoit.vie@meteo.fr>
Date: Mon, 29 Aug 2022 10:44:07 +0200
Subject: [PATCH] Bugfix lambda computation for rain in case NMOM_R=1
---
src/MNH/ini_lima_warm.f90 | 2 +-
src/MNH/lima_rain_evaporation.f90 | 3 +--
2 files changed, 2 insertions(+), 3 deletions(-)
diff --git a/src/MNH/ini_lima_warm.f90 b/src/MNH/ini_lima_warm.f90
index 5ef1977aa..38664ac88 100644
--- a/src/MNH/ini_lima_warm.f90
+++ b/src/MNH/ini_lima_warm.f90
@@ -166,7 +166,7 @@ XLBEXC = 1.0/XBC
IF (NMOM_R.EQ.1) THEN
XCCR = 8.E6
XCXR = -1.
- XLBEXR = 1.0/XBR
+ XLBEXR = 1.0/(XCXR-XBR)
XLBR = ( XAR*XCCR*MOMG(XALPHAR,XNUR,XBR) )**(-XLBEXR)
ELSE
XLBR = XAR*ZGAMR(2)
diff --git a/src/MNH/lima_rain_evaporation.f90 b/src/MNH/lima_rain_evaporation.f90
index f4c0df7d0..c7211f2fc 100644
--- a/src/MNH/lima_rain_evaporation.f90
+++ b/src/MNH/lima_rain_evaporation.f90
@@ -150,8 +150,7 @@ ELSE
!
! Compute the function G(T)
!
- ZZW2(:) = 1. / ( XRHOLW*((((PLV(:)/PT(:))**2)/(XTHCO*XRV)) + & ! G
- (XRV*PT(:))/(XDIVA*PEVSAT(:))))
+ ZZW2(:) = 1. / ( XRHOLW*((((PLV(:)/PT(:))**2)/(XTHCO*XRV)) + (XRV*PT(:))/(XDIVA*PEVSAT(:)))) !G
!
! Compute the evaporation tendency
!
--
GitLab