From 3817d727160d70682561098b3dac79cb06e2601c Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Beno=C3=AEt=20Vi=C3=A9?= <benoit.vie@meteo.fr>
Date: Fri, 26 Aug 2022 11:51:45 +0200
Subject: [PATCH] bugfix low values of ri everywhere
---
src/MNH/ini_lima.f90 | 2 +-
src/MNH/sources_neg_correct.f90 | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/MNH/ini_lima.f90 b/src/MNH/ini_lima.f90
index d90f1e160..74235232d 100644
--- a/src/MNH/ini_lima.f90
+++ b/src/MNH/ini_lima.f90
@@ -147,7 +147,7 @@ ALLOCATE( XCTMIN(7) )
XCTMIN(1) = 1.0 ! Not used
XCTMIN(2) = 1.0E-3 ! Nc
XCTMIN(3) = 1.0E-3 ! Nr
-XCTMIN(4) = 1.0E-3 ! Ni
+XCTMIN(4) = 1.0E-0 ! Ni
XCTMIN(5) = 1.0E-3 ! Not used
XCTMIN(6) = 1.0E-3 ! Not used
XCTMIN(7) = 1.0E-3 ! Not used
diff --git a/src/MNH/sources_neg_correct.f90 b/src/MNH/sources_neg_correct.f90
index aa5978a54..d54bdac40 100644
--- a/src/MNH/sources_neg_correct.f90
+++ b/src/MNH/sources_neg_correct.f90
@@ -275,7 +275,7 @@ CLOUD: select case ( hcloud )
! Correction where ri<0 or Ni<0
if ( lcold_lima ) then
zmask(:,:,:)=(prrs(:, :, :, 4) < xrtmin_lima(4) / ptstep)
- if (nmom_i.ge.2) zmask(:,:,:)=(zmask(:,:,:) .or. prsvs(:, :, :, nsv_lima_ni) < 0. )
+ if (nmom_i.ge.2) zmask(:,:,:)=(zmask(:,:,:) .or. prsvs(:, :, :, nsv_lima_ni) < xctmin_lima(4) / ptstep)
where ( zmask(:,:,:) )
prrs(:, :, :, 1) = prrs(:, :, :, 1) + prrs(:, :, :, 4)
prths(:, :, :) = prths(:, :, :) - prrs(:, :, :, 4) * zls(:, :, :) / &
--
GitLab