Changeset 5647 for LMDZ6/trunk

Timestamp:

May 10, 2025, 9:25:12 PM (4 weeks ago)

Author:

evignon

Message:

petite correction pour la phase mixte dans les thermiques, nouvelle param Lea Etienne

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90

@@ -149 +149 @@
   REAL, DIMENSION(klon,klev+1),    INTENT(IN)   :: tke             ! turbulent kinetic energy [m2/s2]
   REAL, DIMENSION(klon,klev+1),    INTENT(IN)   :: tke_dissip      ! TKE dissipation [m2/s3]
-  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: entr_therm      ! thermal plume entrainment rate [kg/s/m2] ! per mesh surface unit
-  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: detr_therm      ! thermal plume detrainment rate [kg/s/m2] ! per mesh surface unit
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: entr_therm      ! thermal plume entrainment rate * dz [kg/s/m2] 
+  REAL, DIMENSION(klon,klev),      INTENT(IN)   :: detr_therm      ! thermal plume detrainment rate * dz [kg/s/m2]
 
 
@@ -312 +312 @@
   REAL, DIMENSION(klon) :: znebprecip, znebprecipclr, znebprecipcld
   REAL, DIMENSION(klon) :: tot_zneb
-  REAL, DIMENSION(klon) :: zdistcltop, ztemp_cltop
+  REAL, DIMENSION(klon) :: zdistcltop, ztemp_cltop, zdeltaz
   REAL, DIMENSION(klon) :: zqliq, zqice, zqvapcl, zqliqth, zqiceth, zqvapclth, sursat_e, invtau_e ! for icefrac_lscp_turb
 
@@ -658 +658 @@
          Tbefthm1(:)=Tbefth(:)       
         ENDIF
-        zqth=qta(:,k)
+        zqth(:)=qta(:,k)
+        zdeltaz(:)=(paprs(:,k)-paprs(:,k+1))/RG/zp(:)*RD*zt(:)
 
         ! Treatment of stratiform clouds (lognormale or ice-sursat) or all clouds (including cloudth
@@ -849 +850 @@
                         ENDIF 
                         sursat_e(i)=cloudth_senv(i,k)/zqsi(i)
-                        invtau_e(i)=gamma_mixth*MAX(entr_therm(i,k)-detr_therm(i,k),0.)*RD*Tbef(i)/zp(i)
+                        invtau_e(i)=gamma_mixth*MAX(entr_therm(i,k)-detr_therm(i,k),0.)*RD*Tbef(i)/zp(i)/zdeltaz(i)
                      ENDDO
                      CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbefth, zp, paprs(:,k), paprs(:,k+1), wth(:,k), zqi_ini,  &
@@ -860 +861 @@
                         IF (k .GT. 1) THEN
                            ! evaluate the mixing sursaturation using saturation deficit at level below
+                           ! as air pacels detraining into clouds have not (less) seen yet entrainement from above
                            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)
+                           ! mixing is assumed to scales with intensity of net detrainment/entrainment rate (D/dz-E/dz) / rho
+                           invtau_e(i)=gamma_mixth*MAX(detr_therm(i,k)-entr_therm(i,k),0.)*RD*Tbef(i)/zp(i)/zdeltaz(i)
                         ELSE
                            sursat_e(i)=0.