Changeset 2414


Ignore:
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:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/testing/libf/phylmd/conf_phys_m.F90

    r2408 r2414  
    202202    LOGICAL, SAVE :: ok_conserv_q_omp
    203203    INTEGER, SAVE :: iflag_fisrtilp_qsat_omp
     204    INTEGER, SAVE :: iflag_bergeron_omp
    204205    LOGICAL,SAVE :: ok_strato_omp
    205206    LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp
     
    754755
    755756    !Config  Key  = iflag_fisrtilp_qsat
    756     !Config  Desc = Flag de convection
    757     !Config  Def  = 1
    758     !Config  Help = Flag  pour la convection les options suivantes existent :
    759     !Config         -1 pour Kinetic energy correction
    760     !Config         1  conservation kinetic and enthalpy
     757    !Config  Desc = Flag de fisrtilp
     758    !Config  Def  = 0
     759    !Config  Help = Flag  pour la pluie grande-échelle les options suivantes existent :
     760    !Config         >1 nb iterations pour converger dans le calcul de qsat
    761761    iflag_fisrtilp_qsat_omp = 0
    762762    CALL getin('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat_omp)
     763
     764    !Config  Key  = iflag_bergeron
     765    !Config  Desc = Flag de fisrtilp
     766    !Config  Def  = 0
     767    !Config  Help = Flag  pour la pluie grande-échelle les options suivantes existent :
     768    !Config         0 pas d effet Bergeron
     769    !Config         1 effet Bergeron pour T<0
     770    iflag_bergeron_omp = 0
     771    CALL getin('iflag_bergeron',iflag_bergeron_omp)
    763772
    764773    !
     
    19581967    ok_conserv_q = ok_conserv_q_omp
    19591968    iflag_fisrtilp_qsat = iflag_fisrtilp_qsat_omp
     1969    iflag_bergeron = iflag_bergeron_omp
    19601970
    19611971    epmax = epmax_omp
     
    22282238    write(lunout,*)'ok_conserv_q=',ok_conserv_q
    22292239    write(lunout,*)'iflag_fisrtilp_qsat=',iflag_fisrtilp_qsat
     2240    write(lunout,*)'iflag_bergeron=',iflag_bergeron
    22302241    write(lunout,*)' epmax = ', epmax
    22312242    write(lunout,*)' ok_adj_ema = ', ok_adj_ema
  • LMDZ5/branches/testing/libf/phylmd/fisrtilp.F90

    r2408 r2414  
    123123  PARAMETER (ztfondue=278.15)
    124124  REAL dzfice(klon)
     125  REAL zsolid
     126 
    125127  !
    126128  LOGICAL appel1er
     
    932934               *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
    933935           zifl(i) = zifl(i)+ MAX(zcond(i)*zfice(i)-zoliqi(i),0.0) &
    934                     *(paprs(i,k)-paprs(i,k+1))/(RG*dtime) 
     936                    *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
    935937     !      zrfl(i) = zrfl(i)+  zpluie                         &
    936938     !          *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
    937939     !      zifl(i) = zifl(i)+  zice                    &
    938      !               *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)                                   
     940     !               *(paprs(i,k)-paprs(i,k+1))/(RG*dtime)
     941         
     942!CR : on prend en compte l'effet Bergeron dans les flux de precipitation
     943           if ((iflag_bergeron.eq.1).and.(zt(i).LT.273.15)) then
     944              zsolid = zrfl(i)
     945              zifl(i) = zifl(i)+zrfl(i)
     946              zrfl(i) = 0.
     947              zt(i)=zt(i)+zsolid*(RG*dtime)/(paprs(i,k)-paprs(i,k+1)) &
     948                      *(RLSTT-RLVTT)/RCPD/(1.0+RVTMP2*zq(i))
     949           endif
     950!RC                                   
    939951
    940952         ENDIF                     
  • LMDZ5/branches/testing/libf/phylmd/fisrtilp.h

    r1910 r2414  
    1515      INTEGER iflag_pdf
    1616      INTEGER iflag_fisrtilp_qsat
     17      INTEGER iflag_bergeron
    1718
    1819      common/comfisrtilp/                                               &
     
    2627     &     ,reevap_ice                                                  &
    2728     &     ,iflag_fisrtilp_qsat                                         &
     29     &     ,iflag_bergeron                                              &
    2830     &     ,iflag_pdf       
    2931
Note: See TracChangeset for help on using the changeset viewer.