Changeset 5576 for LMDZ6

Timestamp:

Mar 15, 2025, 12:46:52 PM (5 months ago)

Author:

aborella

Message:

Fix for radius of contrails ice crystals (it had the wrong unit)

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_cloud_optics_prop.f90 (modified) (2 diffs)

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

@@ -362 +362 @@
             !--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)
+            !--Beware, re_ice_crystals_contrails is in m, while rei is in microns
+            rei = rei * ( 1. - rcontrail(i,k) ) + re_ice_crystals_contrails * 1.E6 * rcontrail(i,k)
           ENDIF
           pcldtaupi(i, k) = 3.0/2.0*zflwp_var/rad_chaud_pi(i, k) + &
@@ -497 +498 @@
           !--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)
+          !--Beware, re_ice_crystals_contrails is in m, while rei is in microns
+          rei = rei * ( 1. - rcontrail(i,k) ) + re_ice_crystals_contrails * 1.E6 * rcontrail(i,k)
         ENDIF
         pcltau(i, k) = 3.0/2.0*(zflwp_var/rel) + zfiwp_var*(3.448E-03+2.431/ &