Changeset 5633 for LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90

Timestamp:

Apr 28, 2025, 2:46:35 PM (6 weeks ago)

Author:

evignon

Message:

small changes related to mixed phase clouds treatment
Etienne and Lea

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90

@@ -289 +289 @@
   REAL, DIMENSION(klon) :: zdqsdT_raw
   REAL, DIMENSION(klon) :: gammasat,dgammasatdt                   ! coefficient to make cold condensation at the correct RH and derivative wrt T
-  REAL, DIMENSION(klon) :: Tbef,Tbefth,qlibef,DT                  ! temperature, humidity and temp. variation during condensation iteration
+  REAL, DIMENSION(klon) :: Tbef,Tbefth,Tbefthm1,qlibef,DT                  ! temperature, humidity and temp. variation during condensation iteration
   REAL :: num,denom
   REAL :: cste
@@ -536 +536 @@
     ENDIF ! (ok_poprecip)
     
-    ! Calculation of qsat, L/Cp*dqsat/dT and ncoreczq counter
+    ! Calculation of qsat,L/cp*dqsat/dT and ncoreczq counter
     !-------------------------------------------------------
 
@@ -653 +653 @@
         qlibef(:)=0.
         Tbefth(:)=tla(:,k)*pspsk(:,k)
+        IF (k .GT. 1) THEN
+         Tbefthm1(:)=tla(:,k-1)*pspsk(:,k-1)
+        ELSE
+         Tbefthm1(:)=Tbefth(:)       
+        ENDIF
         zqth=qta(:,k)
 
@@ -817 +822 @@
                         invtau_e(i) = 0.
                      ENDDO
-                     CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbef, zp, paprs(:,k), paprs(:,k+1), wls(:,k), zqi_ini, ziflcld, qincloud, &
-                     zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, zficeenv, dzficeenv,                   &
-                     cldfraliq(:,k),sigma2_icefracturb(:,k),mean_icefracturb(:,k))                     
+                     CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbef, zp, paprs(:,k), paprs(:,k+1), wls(:,k), zqi_ini,    &
+                     ziflcld, znebprecipcld, qincloud, zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, & 
+                     zficeenv, dzficeenv, cldfraliq(:,k),sigma2_icefracturb(:,k),mean_icefracturb(:,k))                     
                      DO i=1,klon
                         IF (pticefracturb(i)) THEN
@@ -846 +851 @@
                         invtau_e(i)=gamma_mixth*MAX(entr_therm(i,k)-detr_therm(i,k),0.)*RD*Tbef(i)/zp(i)
                      ENDDO
-                     CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbefth, zp, paprs(:,k), paprs(:,k+1), wth(:,k), zqi_ini, zeroklon, qincloud, &
-                     zcf, zeroklon, zeroklon, sursat_e, invtau_e, zqliqth, zqvapclth, zqiceth, zficeth, dzficeth,                      &
-                     cldfraliqth(:,k), sigma2_icefracturbth(:,k), mean_icefracturbth(:,k))
+                     CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbefth, zp, paprs(:,k), paprs(:,k+1), wth(:,k), zqi_ini,  &
+                     zeroklon, znebprecipcld, qincloud, zcf, zeroklon, zeroklon, sursat_e, invtau_e, zqliqth, zqvapclth, zqiceth, &
+                     zficeth, dzficeth,cldfraliqth(:,k), sigma2_icefracturbth(:,k), mean_icefracturbth(:,k))
                      !Environment
                      DO i=1,klon
                         qincloud(i)=zqn(i)/(1.-fraca(i,k))
                         zcf(i)=MIN(MAX(rneb(i,k)-rnebth(i,k), 0.),1.)*(1.-fraca(i,k))
-                        sursat_e(i)=cloudth_sth(i,k)/zqsith(i)
-                        invtau_e(i)=gamma_mixth*MAX(detr_therm(i,k)-entr_therm(i,k),0.)*RD*Tbef(i)/zp(i)
+                        IF (k .GT. 1) THEN
+                           ! evaluate the mixing sursaturation using saturation deficit at level below
+                           sursat_e(i)=cloudth_sth(i,k-1)/(zqsith(i)+zdqsith(i)*RCPD/RLSTT*(Tbefthm1(i)-Tbefth(i)))
+                           invtau_e(i)=gamma_mixth*MAX(detr_therm(i,k)-entr_therm(i,k),0.)*RD*Tbef(i)/zp(i)
+                        ELSE
+                           sursat_e(i)=0.
+                           invtau_e(i)=0.
+                        ENDIF
                      ENDDO
-                     CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbef, zp, paprs(:,k), paprs(:,k+1), wls(:,k), zqi_ini, ziflcld, qincloud, &
-                     zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, zfice, dzfice,                      &
-                     cldfraliq(:,k),sigma2_icefracturb(:,k), mean_icefracturb(:,k)) 
+                     CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbef, zp, paprs(:,k), paprs(:,k+1), wls(:,k), zqi_ini,    &
+                     ziflcld, znebprecipcld, qincloud, zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, &
+                     zfice, dzfice, cldfraliq(:,k),sigma2_icefracturb(:,k), mean_icefracturb(:,k)) 
   
                     ! adjust zfice to account for condensates in thermals'fraction