Changeset 2414

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)

Location:

LMDZ5/branches/testing/libf/phylmd

Files:

• conf_phys_m.F90 (modified) (4 diffs)
• fisrtilp.F90 (modified) (2 diffs)
• fisrtilp.h (modified) (2 diffs)

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

@@ -202 +202 @@
     LOGICAL, SAVE :: ok_conserv_q_omp
     INTEGER, SAVE :: iflag_fisrtilp_qsat_omp
+    INTEGER, SAVE :: iflag_bergeron_omp
     LOGICAL,SAVE :: ok_strato_omp
     LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp
@@ -754 +755 @@
 
     !Config  Key  = iflag_fisrtilp_qsat
-    !Config  Desc = Flag de convection
-    !Config  Def  = 1
-    !Config  Help = Flag  pour la convection les options suivantes existent :
-    !Config         -1 pour Kinetic energy correction
-    !Config         1  conservation kinetic and enthalpy
+    !Config  Desc = Flag de fisrtilp
+    !Config  Def  = 0
+    !Config  Help = Flag  pour la pluie grande-échelle les options suivantes existent :
+    !Config         >1 nb iterations pour converger dans le calcul de qsat
     iflag_fisrtilp_qsat_omp = 0
     CALL getin('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat_omp)
+
+    !Config  Key  = iflag_bergeron
+    !Config  Desc = Flag de fisrtilp
+    !Config  Def  = 0
+    !Config  Help = Flag  pour la pluie grande-échelle les options suivantes existent :
+    !Config         0 pas d effet Bergeron
+    !Config         1 effet Bergeron pour T<0
+    iflag_bergeron_omp = 0
+    CALL getin('iflag_bergeron',iflag_bergeron_omp)
 
     !
@@ -1958 +1967 @@
     ok_conserv_q = ok_conserv_q_omp
     iflag_fisrtilp_qsat = iflag_fisrtilp_qsat_omp
+    iflag_bergeron = iflag_bergeron_omp
 
     epmax = epmax_omp
@@ -2228 +2238 @@
     write(lunout,*)'ok_conserv_q=',ok_conserv_q
     write(lunout,*)'iflag_fisrtilp_qsat=',iflag_fisrtilp_qsat
+    write(lunout,*)'iflag_bergeron=',iflag_bergeron
     write(lunout,*)' epmax = ', epmax
     write(lunout,*)' ok_adj_ema = ', ok_adj_ema

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                      

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

@@ -15 +15 @@
       INTEGER iflag_pdf
       INTEGER iflag_fisrtilp_qsat
+      INTEGER iflag_bergeron
 
       common/comfisrtilp/                                               &
@@ -26 +27 @@
      &     ,reevap_ice                                                  &
      &     ,iflag_fisrtilp_qsat                                         &
+     &     ,iflag_bergeron                                              &
      &     ,iflag_pdf