Changeset 5602 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp.f90

Timestamp:

Apr 3, 2025, 4:53:58 PM (3 months ago)

Author:

aborella

Message:

Removed deep convection clouds from prognostic cloud properties

File:

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

LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp.f90 ¶

@@ -8 +8 @@
 SUBROUTINE lscp(klon,klev,dtime,missing_val,            &
      paprs, pplay, omega, temp, qt, ql_seri, qi_seri,   &
-     ptconv, ratqs, sigma_qtherm,                       &
+     ptconv, cldfracv, qcondcv, ratqs, sigma_qtherm,    &
      d_t, d_q, d_ql, d_qi, rneb, rneblsvol,             &
      pfraclr, pfracld,                                  &
@@ -153 +153 @@
                                                                    ! CR: if iflag_ice_thermo=2, only convection is active
   LOGICAL, DIMENSION(klon,klev),   INTENT(IN)   :: ptconv          ! grid points where deep convection scheme is active
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: cldfracv        ! cloud fraction from deep convection [-]
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: qcondcv         ! in-cloud condensed specific humidity from deep convection [kg/kg]
 
   !Inputs associated with thermal plumes
@@ -741 +743 @@
 
                     CALL condensation_ice_supersat( &
-                        klon, dtime, missing_val, &
-                        pplay(:,k), paprs(:,k), paprs(:,k+1), &
+                        klon, dtime, pplay(:,k), paprs(:,k), paprs(:,k+1), &
+                        cldfracv(:,k), qcondcv(:,k), &
                         cf_seri(:,k), qvc_seri(:,k), qliq_in, qice_in, &
                         shear, tke_dissip(:,k), cell_area, stratomask(:,k), &
@@ -949 +951 @@
 
     IF (ok_plane_contrail) THEN
-      !--Contrails do not precipitate. We remove then from the variables temporarily
+      !!--Contrails do not precipitate. We remove then from the variables temporarily
+      !DO i = 1, klon
+      !  rneb(i,k) = rneb(i,k) - contfra(i)
+      !  zoliqi(i) = zoliqi(i) - ( qcont(i) - zqs(i) * contfra(i) )
+      !ENDDO
+      !--Contrails precipitate as natural clouds. We save the partition of ice
+      !--between natural clouds and contrails
+      !--NB. we use qcont as a temporary variable to save this partition
       DO i = 1, klon
-        rneb(i,k) = rneb(i,k) - contfra(i)
-        zoliqi(i) = zoliqi(i) - ( qcont(i) - zqs(i) * contfra(i) )
+        IF ( zoliqi(i) .GT. 0. ) THEN
+          qcont(i) = ( qcont(i) - zqs(i) * contfra(i) ) / zoliqi(i)
+        ELSE
+          qcont(i) = 0.
+        ENDIF
       ENDDO
     ENDIF
@@ -990 +1002 @@
     
     IF (ok_plane_contrail) THEN
-      !--Contrails are reintroduced in the variables
+      !!--Contrails are reintroduced in the variables
+      !DO i = 1, klon
+      !  rneb(i,k) = rneb(i,k) + contfra(i)
+      !  zoliqi(i) = zoliqi(i) + ( qcont(i) - zqs(i) * contfra(i) )
+      !ENDDO
+      !--Contrails fraction is left unchanged, but contrails water has changed
       DO i = 1, klon
-        rneb(i,k) = rneb(i,k) + contfra(i)
-        zoliqi(i) = zoliqi(i) + ( qcont(i) - zqs(i) * contfra(i) )
+        IF ( zoliqi(i) .LE. 0. ) THEN
+          contfra(i) = 0.
+          qcont(i) = 0.
+        ELSE
+          qcont(i) = zqs(i) * contfra(i) + zoliqi(i) * qcont(i)
+        ENDIF
       ENDDO
     ENDIF