Changeset 5790 for LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

Timestamp:

Jul 28, 2025, 6:44:28 PM (4 months ago)

Author:

aborella

Message:

Major modifs to treatment of contrails (from 2 classes to 2 moments) + diagnostics. Increased numerical efficiency

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop_ini.f90 (modified) (3 diffs)

LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

@@ -24 +24 @@
   REAL, PROTECTED :: rei_coef, rei_min_temp
   REAL, PROTECTED :: zepsec
-  REAL, PROTECTED :: rei_linear_contrails=7.5E-6 ! [m] effective radius of ice crystals in linear contrails
-  REAL, PROTECTED :: rei_cirrus_contrails=10.E-6 ! [m] effective radius of ice crystals in contrails cirrus
+  REAL, PROTECTED :: eff2vol_radius_contrails=0.7
+  REAL, PROTECTED :: rho_ice=920. ! Ice crystal density (assuming spherical geometry) [kg/m3]
   REAL, PARAMETER :: thres_tau=0.3, thres_neb=0.001
   REAL, PARAMETER :: prmhc=440.*100. ! Pressure between medium and high level cloud in Pa
   REAL, PARAMETER :: prlmc=680.*100. ! Pressure between low and medium level cloud in Pa
   REAL, PARAMETER :: coef_froi=0.09, coef_chau =0.13
-  REAL, PROTECTED :: seuil_neb=0.001
+  REAL, PARAMETER :: seuil_neb=0.001
 ! if iflag_t_glace=0, old values are used for liquid/ice partitionning:
   REAL, PARAMETER :: t_glace_min_old = 258.
@@ -44 +44 @@
 !$OMP THREADPRIVATE(rei_coef, rei_min_temp)
 !$OMP THREADPRIVATE(zepsec) 
-!$OMP THREADPRIVATE(rei_linear_contrails, rei_cirrus_contrails, seuil_neb)
+!$OMP THREADPRIVATE(eff2vol_radius_contrails, rho_ice)
 
   
@@ -110 +110 @@
     rei_min_temp = 175.
     CALL getin_p('rei_min_temp',rei_min_temp)
-    CALL getin_p('rei_linear_contrails', rei_linear_contrails)
-    write(lunout,*)'rei_linear_contrails=',rei_linear_contrails
-    CALL getin_p('rei_cirrus_contrails', rei_cirrus_contrails)
-    write(lunout,*)'rei_cirrus_contrails=',rei_cirrus_contrails
-    CALL getin_p('seuil_neb_rad', seuil_neb)
-    write(lunout,*)'seuil_neb_rad=',seuil_neb
+    CALL getin_p('eff2vol_radius_contrails', eff2vol_radius_contrails)
+    write(lunout,*)'eff2vol_radius_contrails=',eff2vol_radius_contrails