Changeset 6077 for LMDZ6/trunk/libf

Timestamp:

Feb 16, 2026, 2:44:02 PM (4 days ago)

Author:

evignon

Message:

changemeent condition (simplification) appel nouvelle param de phase

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp_main.f90 (modified) (2 diffs)
• lmdz_lscp_phase.f90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90

@@ -746 +746 @@
                ! only in stratiform clouds in the mixed phase regime (Raillard et al. 2025)
                ! it overwrites temperature-dependent phase partitioning except for convective boundary layer clouds
-               pticefracturb(:) = (fraca(:, k) < min_frac_th_cld) .AND. (Tbef(:) .GT. temp_nowater) .AND. (Tbef(:) .LT. RTT)
+               pticefracturb(:) = (fraca(:, k) < min_frac_th_cld) 
                DO i = 1, klon
                   qincloud(i) = zqn(i)
@@ -768 +768 @@
                ! hence we assume that no snow falls in thermals.
                ! we activate only in the mixed phase regime not to distrub the cirrus param at very cold T
-               pticefracturb(:) = (Tbef(:) .GT. temp_nowater) .AND. (Tbef(:) .LT. RTT)
+               pticefracturb(:) = .TRUE.
                !Thermals
                DO i = 1, klon

LMDZ6/trunk/libf/phylmd/lmdz_lscp_phase.f90

@@ -245 +245 @@
          ! because cldfra is intent in, but can be locally modified due to test
          cldfra1D = cldfra(i)
-         ! activate param for concerned grid points and for cloudy conditions
-         IF ((pticefracturb(i)) .AND. (cldfra(i) .GT. 0.)) THEN
-            IF (cldfra(i) .GE. 1.0) THEN
-               cldfra1D = 1.0
-            END IF
-
-            ! T>0°C, no ice allowed
-            IF (temp(i) .GE. RTT) THEN
-               qvap_cld(i) = qsatl(i)*cldfra1D
-               qliq(i) = MAX(0.0, qtot_incl(i) - qsatl(i))*cldfra1D
-               qice(i) = 0.
-               cldfraliq(i) = 1.
-               icefrac(i) = 0.
-               dicefracdT(i) = 0.
-
-               ! T<-38°C, no liquid allowed
-            ELSE IF (temp(i) .LE. temp_nowater) THEN
-               qvap_cld(i) = qsati(i)*cldfra1D
-               qliq(i) = 0.
-               qice(i) = MAX(0.0, qtot_incl(i) - qsati(i))*cldfra1D
-               cldfraliq(i) = 0.
-               icefrac(i) = 1.
-               dicefracdT(i) = 0.
-
-               !---------------------------------------------------------
-               !--             MIXED PHASE TEMPERATURE REGIME
-               !---------------------------------------------------------
-               !--In the mixed phase regime (-38°C< T <0°C) we distinguish
-               !--3 possible subcases.
-               !--1.  No pre-existing ice
-               !--2A. Pre-existing ice and no turbulence
-               !--2B. Pre-existing ice and turbulence
-            ELSE
-
-               ! gamma_snwretro controls the contribution of snowflakes to the negative feedback
-               ! note that for reasons related to inetarctions with the condensation iteration in lscp_main
-               ! we consider here the mean snowflake concentration in the mesh (not the in-cloud concentration)
-               ! when poprecip is active, it will be worth testing considering the incloud fraction, dividing
-               ! by znebprecipcld
-               ! qiceini_incl  = qice_ini(i) / cldfra1D + &
-               !              gamma_snwretro * snowcld(i) * RG * dtime / ( paprsdn(i) - paprsup(i) )
-               ! assuming constant snowfall velocity
+         ! T>0°C, no ice allowed
+         IF (temp(i) .GE. RTT) THEN
+            qvap_cld(i) = qsatl(i)*cldfra1D
+            qliq(i) = MAX(0.0, qtot_incl(i) - qsatl(i))*cldfra1D
+            qice(i) = 0.
+            cldfraliq(i) = 1.
+            icefrac(i) = 0.
+            dicefracdT(i) = 0.
+
+            ! T<-38°C, no liquid allowed
+         ELSE IF (temp(i) .LE. temp_nowater) THEN
+            qvap_cld(i) = qsati(i)*cldfra1D
+            qliq(i) = 0.
+            qice(i) = MAX(0.0, qtot_incl(i) - qsati(i))*cldfra1D
+            cldfraliq(i) = 0.
+            icefrac(i) = 1.
+            dicefracdT(i) = 0.
+          ELSE 
+            !---------------------------------------------------------
+            !--             MIXED PHASE TEMPERATURE REGIME
+            !---------------------------------------------------------
+            !--In the mixed phase regime (-38°C< T <0°C) we distinguish
+            !--3 possible subcases.
+            !--1.  No pre-existing ice
+            !--2A. Pre-existing ice and no turbulence
+            !--2B. Pre-existing ice and turbulence
+
+            ! gamma_snwretro controls the contribution of snowflakes to the negative feedback
+            ! note that for reasons related to inetarctions with the condensation iteration in lscp_main
+            ! we consider here the mean snowflake concentration in the mesh (not the in-cloud concentration)
+            ! when poprecip is active, it will be worth testing considering the incloud fraction, dividing
+            ! by znebprecipcld
+            ! qiceini_incl  = qice_ini(i) / cldfra1D + &
+            !              gamma_snwretro * snowcld(i) * RG * dtime / ( paprsdn(i) - paprsup(i) )
+            ! assuming constant snowfall velocity
+
+            ! activate param for concerned grid points and for cloudy conditions
+            IF ((pticefracturb(i)) .AND. (cldfra(i) .GT. 0.)) THEN
+               IF (cldfra(i) .GE. 1.0) THEN
+                 cldfra1D = 1.0
+               END IF
+
                qiceini_incl = qice_ini(i)/cldfra1D + gamma_snwretro*snowcld(i)/pplay(i)*RD*temp(i)/snow_fallspeed
 
@@ -466 +466 @@
                END IF ! End qini
 
-            END IF ! ! MPC temperature
-
-         END IF ! pticefracturb and cldfra
+            END IF ! pticefracturb and cldfra
+
+         END IF ! MPC temp
 
       END DO ! klon