Changeset 3245 for LMDZ6/trunk/libf/phylmd/newmicro.F90

Timestamp:

Mar 9, 2018, 5:09:21 PM (6 years ago)

Author:

fhourdin

Message:

Seuil minimum nombre de noyaux de condensation cdnc_min
passe a physiq.def

File:

• LMDZ6/trunk/libf/phylmd/newmicro.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/newmicro.F90

@@ -13 +13 @@
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
+  USE ioipsl_getin_p_mod, ONLY : getin_p
+
   IMPLICIT NONE
   ! ======================================================================
@@ -161 +163 @@
   REAL reliq_pi(klon, klev), reice_pi(klon, klev)
 
+  REAL,SAVE :: cdnc_min=-1.
+  REAL,SAVE :: cdnc_min_m3
+  !$OMP THREADPRIVATE(cdnc_min,cdnc_min_m3)
+
   ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! FH : 2011/05/24
@@ -172 +178 @@
   ! Pour retrouver les résultats numériques de la version d'origine,
   ! on impose 0.71 quand on est proche de 0.71
+
+  if (first) THEN
+      call getin_p('cdnc_min',cdnc_min)
+      cdnc_min_m3=cdnc_min*1E6
+      IF (cdnc_min_m3<0.) cdnc_min_m3=20.E6 ! astuce pour retrocompatibilite
+  ENDIF
 
   d_rei_dt = (rei_max-rei_min)/81.4
@@ -239 +251 @@
         cdnc(i, k) = 10.**(bl95_b0+bl95_b1*log(max(mass_solu_aero(i,k), &
           1.E-4))/log(10.))*1.E6 !-m-3
-        cdnc(i, k) = min(1000.E6, max(20.E6,cdnc(i,k)))
+        cdnc(i, k) = min(1000.E6, max(cdnc_min_m3,cdnc(i,k)))
 
         ! --pre-industrial case
         cdnc_pi(i, k) = 10.**(bl95_b0+bl95_b1*log(max(mass_solu_aero_pi(i,k), &
           1.E-4))/log(10.))*1.E6 !-m-3
-        cdnc_pi(i, k) = min(1000.E6, max(20.E6,cdnc_pi(i,k)))
+        cdnc_pi(i, k) = min(1000.E6, max(cdnc_min_m3,cdnc_pi(i,k)))
 
         ! --present-day case
@@ -645 +657 @@
   END IF !ok_cdnc
 
+  first=.false. !to be sure
+
   RETURN