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

Timestamp:

Apr 22, 2025, 6:03:13 PM (2 months ago)

Author:

aborella

Message:

Correction to coupling with deep convection + correction to cloud dissipation

File:

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

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

@@ -1204 +1204 @@
       pcf_seri(:,k) = perscontfra(:)
       qva_seri(:,k) = zqs(:) * contfra(:)
-      qia_seri(:,k) = qcont(:) - zqs(:) * contfra(:)
       DO i = 1, klon
+        IF ( zoliqi(i) .GT. 0. ) THEN
+          !--The quantity of ice in linear contrails has been reduced by precipitation
+          !--with the same factor as the rest of the clouds
+          qia_seri(i,k) = ( qcont(i) - zqs(i) * contfra(i) ) / zoliqi(i) * radocond(i,k)
+        ELSE
+          qia_seri(i,k) = 0.
+        ENDIF
         IF ( ( rneb(i,k) - cfa_seri(i,k) ) .GT. eps ) THEN
+          !--The in-cloud quantity of ice in persistent contrail cirrus is the same as
+          !--the one from all natural cirrus clouds
           qice_perscont(i,k) = ( radocond(i,k) - qia_seri(i,k) ) &
               * perscontfra(i) / ( rneb(i,k) - cfa_seri(i,k) )
@@ -1219 +1227 @@
       DO i = 1, klon
 
+        !--We save the cloud properties that will be advected
         cf_seri(i,k) = rneb(i,k)
-
+        qvc_seri(i,k) = qvc(i)
+
+        !--Deep convection clouds properties are removed from radiative properties
+        !--outputed from lscp (NB. rneb and radocond are only used for the radiative
+        !--properties and are NOT prognostics)
+        !--We must have iflag_coupl == 5 for this coupling to work
+        IF ( ptconv(i,k) .AND. pt_pron_clds(i) ) THEN
+          rneb(i,k) = rneb(i,k) - cfcon(i,k)
+          radocond(i,k) = radocond(i,k) - qccon(i,k) * cfcon(i,k)
+        ENDIF
+
+        !--If everything was precipitated, the remaining empty cloud is dissipated
+        !--and everything is transfered to the subsaturated clear sky region
+        !--NB. we do not change rneb, as it is a diagnostic only
         IF ( zoliq(i) .LE. 0. ) THEN
-          !--If everything was precipitated, the remaining empty cloud is dissipated
-          !--and everything is transfered to the subsaturated clear sky region
-          !--NB. we do not change rneb, as it is a diagnostic only
           cf_seri(i,k) = 0.
+          qvc_seri(i,k) = 0.
           qvc(i) = 0.
         ENDIF
-
-        qvc_seri(i,k) = qvc(i)
 
         !--Diagnostics