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

Timestamp:

Feb 12, 2025, 10:18:07 AM (3 days ago)

Author:

aborella

Message:

Added new formulation for effective radius for ice crystals in natural clouds + removed a print

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

@@ -22 +22 @@
   REAL, PROTECTED :: rad_chau1, rad_chau2, rad_froid
   REAL, PROTECTED :: rei_max, rei_min
+  REAL, PROTECTED :: rei_coef, rei_min_temp
   REAL, PROTECTED :: zepsec
   REAL, PROTECTED :: re_ice_crystals_contrails=7.5E-6 ! [m] effective radius of ice crystals in contrails
@@ -40 +41 @@
 !$OMP THREADPRIVATE(cdnc_min, cdnc_min_m3, rpi, rg, rd)
 !$OMP THREADPRIVATE(rad_chau1, rad_chau2, rad_froid, rei_max, rei_min)
+!$OMP THREADPRIVATE(rei_coef, rei_min_temp)
 !$OMP THREADPRIVATE(zepsec) 
 !$OMP THREADPRIVATE(re_ice_crystals_contrails) 
@@ -104 +106 @@
     rei_max = 61.29
     CALL getin_p('rei_max',rei_max)
+    rei_coef = 2.4
+    CALL getin_p('rei_coef',rei_coef)
+    rei_min_temp = 175.
+    CALL getin_p('rei_min_temp',rei_min_temp)
     CALL getin_p('re_ice_crystals_contrails', re_ice_crystals_contrails)