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

Timestamp:

Jan 17, 2025, 6:12:48 PM (13 days ago)

Author:

aborella

Message:

Added dependance of rei to contrail fraction

File:

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

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

@@ -11 +11 @@
   LOGICAL, PROTECTED :: ok_cdnc
   LOGICAL, PROTECTED :: ok_icefra_lscp, ok_new_lscp
+  LOGICAL, PROTECTED :: ok_plane_contrail
   REAL, PROTECTED :: bl95_b0, bl95_b1 ! Parameter in B&L 95-Formula
   REAL, ALLOCATABLE :: latitude_deg(:)
@@ -22 +23 @@
   REAL, PROTECTED :: rei_max, rei_min
   REAL, PROTECTED :: zepsec
+  REAL, PROTECTED :: re_ice_crystals_contrails=7.5E-6 ! [m] effective radius of ice crystals in contrails
   REAL, PARAMETER :: thres_tau=0.3, thres_neb=0.001
   REAL, PARAMETER :: prmhc=440.*100. ! Pressure between medium and high level cloud in Pa
@@ -39 +41 @@
 !$OMP THREADPRIVATE(rad_chau1, rad_chau2, rad_froid, rei_max, rei_min)
 !$OMP THREADPRIVATE(zepsec) 
+!$OMP THREADPRIVATE(re_ice_crystals_contrails) 
 
   
@@ -46 +49 @@
        & ok_cdnc_in, bl95_b0_in, &
        & bl95_b1_in, latitude_deg_in, rpi_in, rg_in, rd_in, zepsec_in, novlp_in, &
-       & iflag_ice_thermo_in, ok_new_lscp_in)
+       & iflag_ice_thermo_in, ok_new_lscp_in, &
+       & ok_plane_contrail_in)
 
     USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -56 +60 @@
     INTEGER, INTENT(IN) :: novlp_in, iflag_ice_thermo_in
     LOGICAL, INTENT(IN) :: ok_cdnc_in, ok_new_lscp_in
+    LOGICAL, INTENT(IN) :: ok_plane_contrail_in
     REAL, INTENT(IN) :: bl95_b0_in, bl95_b1_in
     REAL, INTENT(IN) :: latitude_deg_in(klon)
@@ -77 +82 @@
     iflag_ice_thermo = iflag_ice_thermo_in
     ok_new_lscp = ok_new_lscp_in
+    ok_plane_contrail = ok_plane_contrail_in
     
     call getin_p('cdnc_min',cdnc_min)
@@ -98 +104 @@
     rei_max = 61.29
     CALL getin_p('rei_max',rei_max)
+    CALL getin_p('re_ice_crystals_contrails', re_ice_crystals_contrails)