Changeset 5174 for LMDZ6/branches/Amaury_dev

Timestamp:

Sep 9, 2024, 11:03:45 AM (3 months ago)

Author:

abarral

Message:

Merge r5169, r5170, r5171

Location:

LMDZ6/branches/Amaury_dev

Files:

• . (modified) (1 prop)
• libf/phylmd/lmdz_lscp_tools.F90 (modified) (6 diffs)
• libf/phylmd/pbl_surface_mod.F90 (modified) (1 diff)

LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_lscp_tools.F90

@@ -233 +233 @@
     !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     ! Compute the liquid, ice and vapour content (+ice fraction) based
-    ! on turbulence (see Fields 2014, Furtado 2016)
+    ! on turbulence (see Fields 2014, Furtado 2016, Raillard 2025)
+  ! L.Raillard (30/08/24)
     !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
@@ -363 +364 @@
             ! Enough TKE
           ELSE
+            print*,"MOUCHOIRACTIVE"
             !---------------------------------------------------------
             !--               ICE SUPERSATURATION PDF
@@ -390 +392 @@
             IF (qiceini_incl > eps) THEN
               nb_crystals = 1.e3 * 5.94e-5 * (RTT - temp(i))**3.33 * naero5**(0.0264 * (RTT - temp(i)) + 0.0033)
-              lambda_PSD = ((RPI * rho_ice * nb_crystals * 24.) / (6. * (qiceini_incl + gamma_snwretro * qsnowcld_incl))) ** (1. / 3.)
+              lambda_PSD = ((RPI * rho_ice * nb_crystals ) / (rho_air * (qiceini_incl + gamma_snwretro * qsnowcld_incl))) ** (1. / 3.)
               N0_PSD = nb_crystals * lambda_PSD
-              moment1_PSD = N0_PSD / 2. / lambda_PSD**2
+              moment1_PSD = N0_PSD /  lambda_PSD**2
             ELSE
               moment1_PSD = 0.
@@ -443 +445 @@
             mean_pdf = sursat_env * 1. / tau_mixingenv / (invtau_phaserelax + 1. / tau_mixingenv)
             cldfraliq(i) = 0.5 * (1. - erf((sursat_iceliq - mean_pdf) / (SQRT(2. * sigma2_pdf))))
-            !IF (cldfraliq(i) .GT. liqfra_max) THEN
-            !    cldfraliq(i) = liqfra_max
-            !ENDIF
+            IF (cldfraliq(i) .GT. liqfra_max) THEN
+                cldfraliq(i) = liqfra_max
+            ENDIF
 
             qliq_incl = qsati(i) * SQRT(sigma2_pdf) / SQRT(2. * RPI) * EXP(-1. * (sursat_iceliq - mean_pdf)**2. / (2. * sigma2_pdf))  &
@@ -452 +454 @@
             sigma2_icefracturb(i) = sigma2_pdf
             mean_icefracturb(i) = mean_pdf
-            !------------ ICE AMOUNT AND WATER CONSERVATION  ------------
+            !------------ SPECIFIC VAPOR CONTENT AND WATER CONSERVATION  ------------
 
             IF ((qliq_incl <= eps) .OR. (cldfraliq(i) <= eps)) THEN
@@ -459 +461 @@
             END IF
 
-            !--Choice for in-cloud vapor :
-            !--1.Weighted mean between qvap in MPC parts and in ice-only parts
-            !--2.Always at ice saturation
+            !--Specific humidity is the max between qsati and the weighted mean between
+              !--qv in MPC patches and qv in ice-only parts. We assume that MPC parts are
+            !--always at qsatl and ice-only parts slightly subsaturated (qsati*sursat_iceext+1)
+              !--The whole cloud can therefore be supersaturated but never subsaturated.
             qvap_incl = MAX(qsati(i), (1. - cldfraliq(i)) * (sursat_iceext + 1.) * qsati(i) + cldfraliq(i) * qsatl(i))
 

LMDZ6/branches/Amaury_dev/libf/phylmd/pbl_surface_mod.F90

@@ -2701 +2701 @@
           yt1_new = (1. / RCPD) * (AcoefH(j) + BcoefH(j) * y_flux_t1(j) * dtime)
           ytsurf_new(j) = yt1_new - y_flux_t1(j) / (Kech_h(j) * RCPD)
-        ENDDO
+        ! for cases forced in flux and for which forcing in Ts is needed
+            ! to prevent the latter to reach unrealistic value (even if not used,
+            ! Ts is calculated and hgardfou can appear during the calculation
+            ! of surface saturation humidity for example
+            if (ok_forc_tsurf) ytsurf_new(j)=tg
+          ENDDO
 
         DO j = 1, knon