Changeset 4913 for LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

Timestamp:

Apr 21, 2024, 1:29:21 PM (4 weeks ago)

Author:

evignon

Message:

amelioration de poprecip avec:

• formule du rayon des cristaux
• qvap du ciel clair pour l'évap des precip

Lea et Audran

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

@@ -267 +267 @@
   REAL, DIMENSION(klon) :: zlh_solid
   REAL, DIMENSION(klon) :: ztupnew
+  REAL, DIMENSION(klon) :: zqvapclr, zqupnew ! for poprecip evap / subl
   REAL :: zm_solid         ! for liquid -> solid conversion
   REAL, DIMENSION(klon) :: zrflclr, zrflcld
@@ -387 +388 @@
 dqsfreez(:,:) = 0.
 dqsmelt(:,:)  = 0.
+zqupnew(:)    = 0.
+zqvapclr(:)   = 0.
 
 
@@ -413 +416 @@
         zq(i)=qt(i,k)
         IF (.not. iftop) THEN
-           ztupnew(i)=temp(i,k+1)+d_t(i,k+1)
+           ztupnew(i)  = temp(i,k+1) + d_t(i,k+1)
+           zqupnew(i)  = qt(i,k+1) + d_q(i,k+1) + d_ql(i,k+1) + d_qi(i,k+1)
+           !--zqs(i) is the saturation specific humidity in the layer above
+           zqvapclr(i) = MAX(0., qt(i,k+1) + d_q(i,k+1) - rneb(i,k+1) * zqs(i))
         ENDIF
         !c_iso init of iso
@@ -424 +430 @@
             CALL poprecip_precld(klon, dtime, iftop, paprs(:,k), paprs(:,k+1), pplay(:,k), &
                               zt, ztupnew, zq, zmqc, znebprecipclr, znebprecipcld, &
+                              zqvapclr, zqupnew, &
                               zrfl, zrflclr, zrflcld, &
                               zifl, ziflclr, ziflcld, &