Changeset 6021 for LMDZ6/trunk/libf/phylmd/lmdz_call_lscp.f90

Timestamp:

Jan 13, 2026, 11:42:56 AM (5 days ago)

Author:

evignon

Message:

petites modifs pour:

• adapter le diag de rhliq et rhice (revision papier mixed-phase dans ACP)
• mettre les precip totales en inout de call_lscp

File:

• LMDZ6/trunk/libf/phylmd/lmdz_call_lscp.f90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_call_lscp.f90

@@ -19 +19 @@
                         sigma2_icefracturb, sigma2_icefracturbth, &
                         mean_icefracturb, mean_icefracturbth, &
-                        radocond, radicefrac, rain, snow, &
+                        radocond, radicefrac, rain_lsc, snow_lsc, &
+                        rain_fall, snow_fall, &
                         frac_impa, frac_nucl, beta, &
                         pfrac_impa, pfrac_nucl, pfrac_1nucl, &
@@ -124 +125 @@
       REAL, DIMENSION(klon, klev), INTENT(INOUT)   :: pfrac_nucl   ! product of nucleation scavenging coeff.
       REAL, DIMENSION(klon, klev), INTENT(INOUT)   :: pfrac_1nucl  ! product of nucleation scavenging coeff. (alpha=1)
+      REAL, DIMENSION(klon), INTENT(INOUT)   :: rain_fall    ! total surface rainfall [kg/m2/s]
+      REAL, DIMENSION(klon), INTENT(INOUT)   :: snow_fall    ! total surface snowfall [kg/m2/s]
 
       ! INPUT/OUTPUT condensation and ice supersaturation
@@ -158 +161 @@
       REAL, DIMENSION(klon, klev), INTENT(OUT)  :: radicefrac       ! ice fraction of condensed water for radiation scheme
       REAL, DIMENSION(klon, klev), INTENT(OUT)  :: rhcl             ! clear-sky relative humidity [-]
-      REAL, DIMENSION(klon), INTENT(OUT)  :: rain             ! surface large-scale rainfall [kg/s/m2]
-      REAL, DIMENSION(klon), INTENT(OUT)  :: snow             ! surface large-scale snowfall [kg/s/m2]
+      REAL, DIMENSION(klon), INTENT(OUT)  :: rain_lsc             ! surface large-scale rainfall [kg/s/m2]
+      REAL, DIMENSION(klon), INTENT(OUT)  :: snow_lsc             ! surface large-scale snowfall [kg/s/m2]
       REAL, DIMENSION(klon, klev + 1), INTENT(OUT)  :: prfl             ! large-scale rainfall flux in the column [kg/s/m2]
       REAL, DIMENSION(klon, klev + 1), INTENT(OUT)  :: psfl             ! large-scale snowfall flux in the column [kg/s/m2]
@@ -276 +279 @@
                    sigma2_icefracturb, sigma2_icefracturbth, &
                    mean_icefracturb, mean_icefracturbth, &
-                   radocond, radicefrac, rain, snow, &
+                   radocond, radicefrac, rain_lsc, snow_lsc, &
                    frac_impa, frac_nucl, beta, &
                    prfl, psfl, rhcl, &
@@ -303 +306 @@
                        temp, qt, ptconv, ratqs, sigma_qtherm, &
                        d_t, d_q, d_ql, d_qi, rneb, rneblsvol, radocond, &
-                       rain, snow, &
+                       rain_lsc, snow_lsc, &
                        pfrac_impa, pfrac_nucl, pfrac_1nucl, &
                        frac_impa, frac_nucl, beta, &
@@ -318 +321 @@
 
       ! rain and snow are set to 0 when negative
-      WHERE (rain < 0) rain = 0.
-      WHERE (snow < 0) snow = 0.
+      WHERE (rain_lsc < 0) rain_lsc = 0.
+      WHERE (snow_lsc < 0) snow_lsc = 0.
 
       ! so-called 'numerical rain' is computed when qlnew=ql+dql>qlmax and qinew=qi+dqi>qimax
@@ -346 +349 @@
       END IF
 
+      ! Total precipitation 
+      DO i = 1, klon
+         rain_fall(i) = rain_fall(i) + rain_lsc(i)
+         snow_fall(i) = snow_fall(i) + snow_lsc(i)
+      ENDDO
+
+
 !===============================================================================
 ! Add tendencies