Changeset 4823

Timestamp:

Feb 15, 2024, 9:58:39 PM (3 months ago)

Author:

evignon

Message:

correction des seuils du melting et freezing dans lmdz_lscp_poprecip

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90

@@ -691 +691 @@
                    * capa_snowflake / RLMLT * coef_ventil &
                    * MAX(0., temp_wetbulb - RTT) * dtime
-    ENDIF
+      dqsclrmelt = MAX( dqsclrmelt , -snowclr(i) / hum_to_flux(i))
+      ELSE
+      dqsclrmelt = 0.
+   ENDIF
 
     !--In cloudy air
@@ -704 +707 @@
                    * capa_snowflake / RLMLT * coef_ventil &
                    * MAX(0., temp_wetbulb - RTT) * dtime
+      dqscldmelt = MAX( dqscldmelt , -snowcld(i) / hum_to_flux(i))
+    ELSE
+      dqscldmelt = 0.
     ENDIF
     
@@ -765 +771 @@
     tau_freez = 1. / ( gamma_freez &
               * EXP( - alpha_freez * ( temp(i) - temp_nowater ) / ( RTT - temp_nowater ) ) )
-    !--Exact solution of dqrain/dt = -qrain/tau_freez
-    dqrtotfreez = raincld(i) / hum_to_flux(i) * ( EXP( - dtime / tau_freez ) - 1. )
-
-    !--Barrier
-    !--It is a MAX because everything is negative
-    dqrtotfreez = MAX(dqrtotfreez, dqrfreez_max)
-
-    !--The partition between clear and cloudy air
-    !--is proportionnal to the rain fluxes in clear / cloudy air
-    dqrclrfreez = dqrtotfreez * rainclr(i) / ( rainclr(i) + raincld(i) )
-    dqrcldfreez = dqrtotfreez - dqrclrfreez
+
+    !--In clear air
+    IF ( rainclr(i) .GT. 0. ) THEN
+       !--Exact solution of dqrain/dt = -qrain/tau_freez
+       dqrclrfreez = rainclr(i) / hum_to_flux(i) * ( EXP( - dtime / tau_freez ) - 1. )
+       dqrclrfreez = MAX(dqrclrfreez, -rainclr(i) / hum_to_flux(i))
+    ELSE
+       dqrclrfreez = 0.
+    ENDIF
+
+    !--In cloudy air
+    IF ( raincld(i) .GT. 0. ) THEN
+       !--Exact solution of dqrain/dt = -qrain/tau_freez
+       dqrcldfreez = raincld(i) / hum_to_flux(i) * ( EXP( - dtime / tau_freez ) - 1. )
+       dqrcldfreez = MAX(dqrcldfreez, -raincld(i) / hum_to_flux(i))
+    ELSE
+       dqrcldfreez = 0.
+    ENDIF
+
+    !--Barriers
+    !--It is activated if the total is LOWER than the max
+    !--because everything is negative
+    dqrtotfreez = dqrclrfreez + dqrcldfreez
+    IF ( dqrtotfreez .LT. dqrfreez_max ) THEN
+      !--We redistribute the max freezed rain keeping
+      !--the clear/cloud partition of the freezing rain
+      dqrclrfreez = dqrfreez_max * dqrclrfreez / dqrtotfreez
+      dqrcldfreez = dqrfreez_max * dqrcldfreez / dqrtotfreez
+      dqrtotfreez = dqrfreez_max
+    ENDIF
+
 
     !--Add tendencies