Changeset 5488 for LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop.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.f90 (modified) (5 diffs)

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

@@ -9 +9 @@
     reliq_pi, reice_pi, scdnc, cldncl, reffclwtop, lcc, reffclws, &
     reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra,  & 
-    icefrac_optics, dNovrN, ptconv,rnebcon, ccwcon)
+    icefrac_optics, dNovrN, ptconv, rnebcon, ccwcon, rcontrail)
 
   USE lmdz_cloud_optics_prop_ini , ONLY : flag_aerosol, ok_cdnc
@@ -28 +28 @@
   USE lmdz_cloud_optics_prop_ini , ONLY : ok_icefra_lscp, rei_max, rei_min
   USE lmdz_cloud_optics_prop_ini , ONLY : zepsec, novlp, iflag_ice_thermo, ok_new_lscp
+  USE lmdz_cloud_optics_prop_ini , ONLY : ok_plane_contrail, re_ice_crystals_contrails
   
 
@@ -69 +70 @@
 
   LOGICAL, INTENT(IN) :: ptconv(klon, klev) ! flag for grid points affected by deep convection
+
+  REAL, INTENT(IN) :: rcontrail(klon, klev) ! ratio of contrails to total cloud fraction, used only if ok_plane_contrail=y [-]
 
   ! inout:
@@ -333 +336 @@
 
           IF (zfiwp_var==0. .OR. rei<=0.) rei = 1.
+          IF ( ok_plane_contrail ) THEN
+            !--If contrails are activated, rei is a weighted average between the natural
+            !--rei and the contrails rei, with the weights being the fraction of natural
+            !--vs contrail cirrus in the gridbox
+            rei = rei * ( 1. - rcontrail(i,k) ) + re_ice_crystals_contrails * rcontrail(i,k)
+          ENDIF
           pcldtaupi(i, k) = 3.0/2.0*zflwp_var/rad_chaud_pi(i, k) + &
             zfiwp_var*(3.448E-03+2.431/rei)
@@ -442 +451 @@
         IF (zflwp_var==0.) rel = 1.
         IF (zfiwp_var==0. .OR. rei<=0.) rei = 1.
+        IF ( ok_plane_contrail ) THEN
+          !--If contrails are activated, rei is a weighted average between the natural
+          !--rei and the contrails rei, with the weights being the fraction of natural
+          !--vs contrail cirrus in the gridbox
+          rei = rei * ( 1. - rcontrail(i,k) ) + re_ice_crystals_contrails * rcontrail(i,k)
+        ENDIF
         pcltau(i, k) = 3.0/2.0*(zflwp_var/rel) + zfiwp_var*(3.448E-03+2.431/ &
           rei)