Changeset 5790 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_precip.f90

Timestamp:

Jul 28, 2025, 6:44:28 PM (19 hours ago)

Author:

aborella

Message:

Major modifs to treatment of contrails (from 2 classes to 2 moments) + diagnostics. Increased numerical efficiency

File:

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

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

@@ -18 +18 @@
 SUBROUTINE histprecip_precld( &
            klon, dtime, iftop, paprsdn, paprsup, pplay, zt, ztupnew, zq, &
-           zmqc, zneb, znebprecip, znebprecipclr, flsed, flsed_lincont, flsed_circont, &
+           zmqc, zneb, znebprecip, znebprecipclr, flsed, flsed_cont, &
            zrfl, zrflclr, zrflcld, zifl, ziflclr, ziflcld, dqreva, dqssub &
            )
@@ -45 +45 @@
 REAL,    INTENT(INOUT), DIMENSION(klon) :: zmqc           !--specific humidity in the precipitation falling from the upper layer [kg/kg]
 REAL,    INTENT(IN),    DIMENSION(klon) :: flsed          !--sedimentated ice flux [kg/s/m2]
-REAL,    INTENT(IN),    DIMENSION(klon) :: flsed_lincont  !--linear contrails sedimentated ice flux [kg/s/m2]
-REAL,    INTENT(IN),    DIMENSION(klon) :: flsed_circont  !--contrail cirrus sedimentated ice flux [kg/s/m2]
+REAL,    INTENT(IN),    DIMENSION(klon) :: flsed_cont     !--contrails sedimentated ice flux [kg/s/m2]
 
 REAL,    INTENT(IN),    DIMENSION(klon) :: zneb           !--cloud fraction IN THE LAYER ABOVE [-]
@@ -134 +133 @@
     !--If the sedimentation of ice crystals is activated, the falling ice is sublimated and
     !--added to the total water content of the gridbox
-    IF ( (flsed(i) + flsed_lincont(i) + flsed_circont(i)) .GT. 0. ) THEN
-      qice_sedim = (flsed(i) + flsed_lincont(i) + flsed_circont(i)) &
+    IF ( (flsed(i) + flsed_cont(i)) .GT. 0. ) THEN
+      qice_sedim = (flsed(i) + flsed_cont(i)) &
           / ( paprsdn(i) - paprsup(i) ) * RG * dtime
 
@@ -757 +756 @@
            klon, dtime, iftop, paprsdn, paprsup, pplay, temp, tempupnew, qvap, &
            qprecip, precipfracclr, precipfraccld, qvapclrup, qtotupnew, &
-           flsed, flsed_lincont, flsed_circont, &
+           flsed, flsed_cont, &
            cldfra, qvc, qliq, qice, &
            rain, rainclr, raincld, snow, snowclr, snowcld, &
@@ -793 +792 @@
 REAL,    INTENT(IN),    DIMENSION(klon) :: qtotupnew      !--total specific humidity IN THE LAYER ABOVE [kg/kg]
 REAL,    INTENT(IN),    DIMENSION(klon) :: flsed          !--sedimentated ice water flux [kg/s/m2]
-REAL,    INTENT(IN),    DIMENSION(klon) :: flsed_lincont  !--sedimentated ice water flux [kg/s/m2]
-REAL,    INTENT(IN),    DIMENSION(klon) :: flsed_circont  !--sedimentated ice water flux [kg/s/m2]
+REAL,    INTENT(IN),    DIMENSION(klon) :: flsed_cont     !--sedimentated ice water flux [kg/s/m2]
 
 REAL,    INTENT(INOUT), DIMENSION(klon) :: cldfra         !--cloud fraction at the beginning of lscp - used only if the cloud properties are advected [-]
@@ -896 +894 @@
   cpw = RCPD * RVTMP2
   DO i = 1, klon
-    IF ( (flsed(i) + flsed_lincont(i) + flsed_circont(i)) .GT. 0. ) THEN
-      qice_sedim = (flsed(i) + flsed_lincont(i) + flsed_circont(i)) / dhum_to_dflux(i)
+    IF ( (flsed(i) + flsed_cont(i)) .GT. 0. ) THEN
+      qice_sedim = (flsed(i) + flsed_cont(i)) / dhum_to_dflux(i)
       !--No condensed water so cp=cp(vapor+dry air)
       !-- RVTMP2=rcpv/rcpd-1
@@ -1005 +1003 @@
     !--If the sedimentation of ice crystals is activated, the falling ice is sublimated and
     !--added to the total water content of the gridbox
-    IF ( (flsed(i) + flsed_lincont(i) + flsed_circont(i)) .GT. 0. ) THEN
-      qice_sedim = (flsed(i) + flsed_lincont(i) + flsed_circont(i)) / dhum_to_dflux(i)
+    IF ( (flsed(i) + flsed_cont(i)) .GT. 0. ) THEN
+      qice_sedim = (flsed(i) + flsed_cont(i)) / dhum_to_dflux(i)
       !--Vapor is updated after evaporation/sublimation (it is increased)
       qvap(i) = qvap(i) + qice_sedim