Changeset 2414 for LMDZ5/branches/testing/libf/phylmd/fisrtilp.F90

Timestamp:

Dec 23, 2015, 5:05:06 PM (9 years ago)

Author:

crio

Message:

Prise en compte de l'effet Bergeron dans les flux de pluie grande-echelle
Activer en mettant iflag_bergeron=1 dans physiq.def (valeur par défaut 0)
Representation of the Bergeron effect in the large-scale precipitation fluxes
Activate by iflag_bergeron=1 in physiq.def (standard value 0)

File:

• LMDZ5/branches/testing/libf/phylmd/fisrtilp.F90 (modified) (2 diffs)

LMDZ5/branches/testing/libf/phylmd/fisrtilp.F90

@@ -123 +123 @@
   PARAMETER (ztfondue=278.15)
   REAL dzfice(klon)
+  REAL zsolid
+  
   !
   LOGICAL appel1er
@@ -932 +934 @@
                *(paprs(i,k)-paprs(i,k+1))/(RG*dtime) 
            zifl(i) = zifl(i)+ MAX(zcond(i)*zfice(i)-zoliqi(i),0.0) &
-                    *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)  
+                    *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
      !      zrfl(i) = zrfl(i)+  zpluie                         &
      !          *(paprs(i,k)-paprs(i,k+1))/(RG*dtime) 
      !      zifl(i) = zifl(i)+  zice                    &
-     !               *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)                                    
+     !               *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
+         
+!CR : on prend en compte l'effet Bergeron dans les flux de precipitation
+           if ((iflag_bergeron.eq.1).and.(zt(i).LT.273.15)) then
+              zsolid = zrfl(i)
+              zifl(i) = zifl(i)+zrfl(i)
+              zrfl(i) = 0.
+              zt(i)=zt(i)+zsolid*(RG*dtime)/(paprs(i,k)-paprs(i,k+1)) &
+                      *(RLSTT-RLVTT)/RCPD/(1.0+RVTMP2*zq(i))
+           endif
+!RC                                    
 
          ENDIF