Changeset 5496 for LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

Timestamp:

Jan 20, 2025, 10:37:37 PM (12 days ago)

Author:

evignon

Message:

nouvelle option pour le calcul du rayon optique des cristaux de glace nuageuse
qui inclue une dependance en temperature ET en contenu en glace (contrairement a la formule
utilisee aujourdhui qui n'a qu'une dependance en T) La formule derivee ici est relativement simple
avec deux coefficients de tuning dont les plages de valeurs acceptables sont donnes en commentaires.
Elle se compare tres bien avec la formule " de reference " de Sun and Rikkus 1999 et Sun 2001

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_cloud_optics_prop_ini.f90

@@ -21 +21 @@
   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, PARAMETER :: thres_tau=0.3, thres_neb=0.001
@@ -38 +39 @@
 !$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) 
 
@@ -98 +100 @@
     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)