Changeset 5933

Timestamp:

Dec 16, 2025, 10:03:43 AM (4 weeks ago)

Author:

evignon

Message:

petite modif de la param des nuages de phase suite aux reponses aux reviews du papier de Lea

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_lscp_tools.f90

@@ -310 +310 @@
    REAL :: rho_air
    REAL :: psati                                                          !--saturation vapor pressure wrt ice             [Pa]
-
+   REAL :: sigmaw2                                                        !--variance of vertical turbulent velocity       [m2/s2]
                                                                        
     REAL :: tempvig1, tempvig2
@@ -334 +334 @@
     DO i=1,klon
      rho_air  = pplay(i) / temp(i) / RD
+     ! assuming turbulence isotropy, tke=3/2*sigmaw2
+     sigmaw2=2./3*tke(i)
      ! because cldfra is intent in, but can be locally modified due to test 
      cldfra1D = cldfra(i)
@@ -384 +386 @@
            !--cloud else fully iced cloud
            IF ( (qiceini_incl .LT. eps) .AND. (invtau_e(i) .LT. eps) ) THEN
-              IF ( (wvel(i) .GT. eps) .OR. (tke(i) .GT. eps) ) THEN
+              IF ( (wvel(i)+sqrt(sigmaw2) .GT. eps) .OR. (tke(i) .GT. eps) ) THEN
                  qvap_cld(i)   = qsatl(i) * cldfra1D
                  qliq(i)       = MAX(0.,qtot_incl(i)-qsatl(i)) * cldfra1D
@@ -484 +486 @@
                  !--Tau_dissipturb is the time needed for turbulence to decay 
                  !--due to viscosity
-                 tau_dissipturb = gamma_taud * 2. * 2./3. * tke(i) / tke_dissip(i) / C0
+                 tau_dissipturb = gamma_taud * 2. * sigmaw2 / tke_dissip(i) / C0
 
                  !--------------------- PDF COMPUTATIONS ---------------------
@@ -495 +497 @@
 
                  liqfra_max = MAX(0., (MIN (1.,( qtot_incl(i) - (qice_ini(i) / cldfra1D) - qsati(i) * (1 + sursat_iceext ) ) / ( qsatl(i) - qsati(i) ) ) ) )
-                 sigma2_pdf = 1./2. * ( ai**2 ) *  2./3. * tke(i) * tau_dissipturb / (invtau_phaserelax + invtau_e(i))
+                 sigma2_pdf = 1./2. * ( ai**2 ) *  sigmaw2 * tau_dissipturb / (invtau_phaserelax + invtau_e(i))
                  ! sursat ranges between -1 and 1, so we prevent sigma2 so exceed 1
                  cldfraliq(i) = 0.5 * (1. - erf( ( sursat_iceliq - sursat_equ) / (SQRT(2.* sigma2_pdf) ) ) )