Changeset 4888

Timestamp:

Apr 2, 2024, 9:29:46 PM (8 months ago)

Author:

evignon

Message:

ajout d'une condition sur les fractions precipitantes pour la fonte et le freezing dans poprecip

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90

@@ -361 +361 @@
 !--Initialisation of variables
 
+
 qzero(:)    = 0.
 
@@ -382 +383 @@
   dqlauto = 0.
   dqlrim  = 0.
-  dqifreez= 0.
 
   !--dhum_to_dflux: coef to convert a delta in specific quantity to a flux 
@@ -710 +710 @@
     air_thermal_conduct=(5.69+0.017*(temp(i)-RTT))*1.e-3*4.184                ! thermal conductivity of the air, SI
 
+
     !--In clear air
-    IF ( snowclr(i) .GT. 0. ) THEN
+    IF ( ( snowclr(i) .GT. 0. ) .AND.  ( precipfracclr(i) .GT. 0. ) ) THEN
       !--Calculated according to
       !-- flux = velocity_snowflakes * nb_snowflakes * volume_snowflakes * rho_snow
@@ -719 +720 @@
                    * capa_snowflake / RLMLT * coef_ventil &
                    * MAX(0., temp_wetbulb - RTT) * dtime
-
+      
       !--Barrier to limit the melting flux to the clr snow flux in the mesh
       dqsclrmelt = MAX( dqsclrmelt , -snowclr(i) / dhum_to_dflux(i))
     ENDIF
 
+
     !--In cloudy air
-    IF ( snowcld(i) .GT. 0. ) THEN
+    IF ( ( snowcld(i) .GT. 0. ) .AND. ( precipfraccld(i) .GT. 0. ) ) THEN
       !--Calculated according to
       !-- flux = velocity_snowflakes * nb_snowflakes * volume_snowflakes * rho_snow
@@ -736 +738 @@
       !--Barrier to limit the melting flux to the cld snow flux in the mesh
       dqscldmelt = MAX(dqscldmelt , - snowcld(i) / dhum_to_dflux(i))
-    ENDIF
+   ENDIF
     
+
     !--Barrier on temperature. If the total melting flux leads to a
     !--positive temperature, it is limited to keep temperature above 0 degC.
@@ -758 +761 @@
     snowclr(i) = MAX(0., snowclr(i) + dqsclrmelt * dhum_to_dflux(i))
     snowcld(i) = MAX(0., snowcld(i) + dqscldmelt * dhum_to_dflux(i))
-
 
     !--Temperature adjustment with the release of latent
@@ -807 +809 @@
     dqrtotfreez_step1 = 0.
 
-    IF ( ( qice(i) .GT. 0. ) .AND. ( raincld(i) .GT. 0. ) ) THEN
+    IF ( ( qice(i) .GT. 0. ) .AND. ( raincld(i) .GT. 0. ) .AND. ( precipfraccld(i) .GT. 0. ) ) THEN
       dqrclrfreez = 0.
       dqrcldfreez = 0.
@@ -913 +915 @@
 
   ENDIF
+
+
 
   !--If the local flux of rain+snow in clear/cloudy air is lower than rain_int_min,
@@ -943 +947 @@
                  + snowcld(i) / ( pplay(i) / RD / temp(i) ) / snow_fallspeed_cld )
 
+
 ENDDO ! loop on klon
 
+
 END SUBROUTINE poprecip_postcld