From 21246adb62dde710520d21d50e60d8f451074bbc Mon Sep 17 00:00:00 2001
From: Pierre Tulet <pierre.tulet@aero.obs-mip.fr>
Date: Mon, 22 May 2023 17:46:27 +0200
Subject: [PATCH] P.Tulet 22/02/23: correction pour eviter d'avoir une
 concentraotin lognormale nulle. correction to avoid having a zero in the
 lognormal concentration.

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

diff --git a/src/MNH/ch_aer_reallfin.f90 b/src/MNH/ch_aer_reallfin.f90
index 505f3acfa..101835c6c 100644
--- a/src/MNH/ch_aer_reallfin.f90
+++ b/src/MNH/ch_aer_reallfin.f90
@@ -134,7 +134,7 @@ ZRHOI(:) = 1.8e3
 ZRHOI(JP_AER_H2O) = 1.0e3   ! water
 ZRHOI(JP_AER_DST) = XDENSITY_DUST
 !
-PSV(:,:,:,:)  = MAX(PSV(:,:,:,:), 0.)
+PSV(:,:,:,:)  = MAX(PSV(:,:,:,:), 1E-20)
 PCO(:,:,:)    = MAX(PCO(:,:,:), 10.E-9*PRHODREF(:,:,:))
 !
 ! Special treatment for BC and OC (link to CO gaseous concentration)
@@ -310,6 +310,7 @@ DO JN=1,JPMODE
 !
   ZRG(:,:,:,JN)=(ZM(:,:,:,NM3(JN))/ZN(:,:,:,JN))**(1./3.)*EXP(-1.5*ZSIGMA(:,:,:))
 !
+
 ENDDO
 !
 !conversion into ppv
-- 
GitLab