From 9f6275d109fc775a14fdfc82548441244734d578 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Beno=C3=AEt=20Vi=C3=A9?= <benoit.vie@meteo.fr>
Date: Thu, 5 May 2022 16:53:07 +0200
Subject: [PATCH] ongoing work on ice nucleation for NMOM_I=1

---
 src/MNH/ini_lima_cold_mixed.f90   | 13 +++++++++++++
 src/MNH/lima_nucleation_procs.f90 |  7 ++++---
 2 files changed, 17 insertions(+), 3 deletions(-)

diff --git a/src/MNH/ini_lima_cold_mixed.f90 b/src/MNH/ini_lima_cold_mixed.f90
index e6daea87e..756164cfb 100644
--- a/src/MNH/ini_lima_cold_mixed.f90
+++ b/src/MNH/ini_lima_cold_mixed.f90
@@ -56,6 +56,7 @@ USE MODD_PARAM_LIMA
 USE MODD_PARAM_LIMA_WARM
 USE MODD_PARAM_LIMA_COLD
 USE MODD_PARAM_LIMA_MIXED
+use MODD_RAIN_ICE_PARAM,  only: XALPHA1, XALPHA2, XBETA1, XBETA2, IMNU0=>XMNU0, XNU10, XNU20
 USE MODD_REF
 !
 use mode_msg
@@ -146,6 +147,18 @@ REAL :: ZRHOIW ! ice density
 !
 ILUOUT0 = TLUOUT0%NLU
 !
+! Ice nucleation for NMOM_I=1
+!
+!
+XNU10 = 50.*ZFACT_NUCL
+XALPHA1 = 4.5
+XBETA1 = 0.6
+!
+XNU20 = 1000.*ZFACT_NUCL
+XALPHA2 = 12.96
+XBETA2 = 0.639
+!
+IMNU0 = 6.88E-13
 !
 !*       1.     CHARACTERISTICS OF THE SPECIES
 !   	        ------------------------------
diff --git a/src/MNH/lima_nucleation_procs.f90 b/src/MNH/lima_nucleation_procs.f90
index 91fb9deaf..b28c91cb6 100644
--- a/src/MNH/lima_nucleation_procs.f90
+++ b/src/MNH/lima_nucleation_procs.f90
@@ -313,9 +313,10 @@ IF (LCOLD .AND. LNUCL .AND. NMOM_I.EQ.1) THEN
 !  Z_TH_HIND=ZTHS*PTSTEP-PTHT
 !  Z_RI_HIND=ZRIS*PTSTEP-PRIT
 !  Z_CI_HIND=ZCIT-PCIT
-!  PCIT=ZCIT
-!  PRIT=ZRIS*PTSTEP
-!  PTHT=ZTHS*PTSTEP
+  PCIT=ZCIT
+  PRIT=ZRIS*PTSTEP
+  PTHT=ZTHS*PTSTEP
+  PRVT=ZRVS*PTSTEP
 !  Z_TH_HINC=0.
 !  Z_RC_HINC=0.
 !  Z_CC_HINC=0.
-- 
GitLab