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

Timestamp:

Jul 24, 2023, 4:46:18 PM (16 months ago)

Author:

evignon

Message:

modification de Lea pour inclure une dependance de la phase a la temperature du sommet du nuage

File:

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

LMDZ6/trunk/libf/phylmd/newmicro.F90

@@ -3 +3 @@
 SUBROUTINE newmicro(flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, paprs, pplay, t, pqlwp, picefra, pclc, &
     pcltau, pclemi, pch, pcl, pcm, pct, pctlwp, xflwp, xfiwp, xflwc, xfiwc, &
-    mass_solu_aero, mass_solu_aero_pi, pcldtaupi, latitude_deg,distcltop, re, fl, reliq, reice, &
+    mass_solu_aero, mass_solu_aero_pi, pcldtaupi, latitude_deg, distcltop, temp_cltop, re, fl, reliq, reice, &
     reliq_pi, reice_pi)
 
@@ -43 +43 @@
   ! latitude_deg-input latitude in degrees
   ! distcltop ---input- distance from cloud top
+  ! temp_cltop  -input- temperature of cloud top
 
   ! re------output-R-Cloud droplet effective radius multiplied by fl [um]
@@ -100 +101 @@
 
   REAL distcltop(klon,klev)
+  REAL temp_cltop(klon,klev)
   LOGICAL lo
 
@@ -244 +246 @@
 
       IF (ok_new_lscp) THEN
-          CALL icefrac_lscp(klon,t(:,k),iflag_ice_thermo,distcltop(:,k),zfice(:,k),dzfice(:,k))
+          CALL icefrac_lscp(klon,t(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k),zfice(:,k),dzfice(:,k))
       ELSE
           CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k))