Changeset 5453 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_condensation.f90

Timestamp:

Dec 23, 2024, 8:19:39 PM (37 hours ago)

Author:

aborella

Message:

Added water emissions and IO routines for contrails.
Work to be done: contrails initial cross section and radiative transfer

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_condensation.f90 (modified) (5 diffs)

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

@@ -157 +157 @@
 !
 REAL,     INTENT(INOUT), DIMENSION(klon) :: rcont_seri    ! ratio of contrails fraction to total cloud fraction [-]
-REAL,     INTENT(IN),    DIMENSION(klon) :: flight_dist   ! 
-REAL,     INTENT(IN),    DIMENSION(klon) :: flight_h2o    ! 
+REAL,     INTENT(IN),    DIMENSION(klon) :: flight_dist   ! aviation distance flown within the mesh [m/s/mesh]
+REAL,     INTENT(IN),    DIMENSION(klon) :: flight_h2o    ! aviation H2O emitted within the mesh [kgH2O/s/mesh]
 !
 !  Output
@@ -238 +238 @@
 !
 ! for cell properties
-REAL :: rho, rhodz, dz
-REAL :: V_cell, M_cell
+REAL :: rho, rhodz, dz, V_cell
 
 qzero(:) = 0.
@@ -336 +335 @@
       !--Cell volume [m3]
       V_cell = dz * cell_area(i)
-      !--Cell dry air mass [kg]
-      M_cell = rhodz * cell_area(i)
 
 
@@ -504 +501 @@
         ENDIF
         pdf_e3 = k0_pdf_lscp + kappa_pdf_lscp * MAX( temp_nowater - temp(i), 0. )
-        pdf_alpha = EXP( rhl_clr / 100. ) * pdf_e3
-        pdf_alpha = MIN(10., pdf_alpha)
+        pdf_alpha = EXP( MIN(1000., rhl_clr) / 100. ) * pdf_e3
         
         IF ( ok_warm_cloud ) THEN
@@ -790 +786 @@
             dtime, pplay(i), temp(i), qsat(i), qsatl(i), gamma_cond(i), &
             rcont_seri(i), flight_dist(i), cldfra(i), qvc(i), &
-            V_cell, M_cell, pdf_loc, pdf_scale, pdf_alpha, &
+            dz, V_cell, pdf_loc, pdf_scale, pdf_alpha, &
             Tcritcont(i), qcritcont(i), potcontfraP(i), potcontfraNP(i), contfra(i), &
             dcf_avi(i), dqvc_avi(i), dqi_avi(i) &