Changeset 4413


Ignore:
Timestamp:
Feb 1, 2023, 5:47:09 PM (22 months ago)
Author:
evignon
Message:

pour les mini-projets: controle des coefficients pour le detrainement
et l'entrainement des downdrafts dans les .def

Location:
LMDZ6/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90

    r4381 r4413  
    88   integer        :: iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure
    99   integer        :: iflag_thermals_down
     10   real           :: entrdn_coef, detrdn_coef
    1011
    1112!$OMP THREADPRIVATE(dvdq,dqimpl,prt_level,lunout)
     
    1415!$OMP THREADPRIVATE(iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure)
    1516!$OMP THREADPRIVATE(iflag_thermals_down)
     17!$OMP THREADPRIVATE(entrdn_coef, detrdn_coef)
    1618
    1719
     
    9193
    9294
     95   entrdn_coef = 0.5
     96   CALL getin_p('entrdn_coef',entrdn_coef)
     97
     98   detrdn_coef = 0.5
     99   CALL getin_p('detrdn_coef',detrdn_coef)
     100
     101
    93102     CALL getin_p('thermals_fact_epsilon',fact_epsilon)
    94103     CALL getin_p('thermals_betalpha',betalpha)
  • LMDZ6/trunk/libf/phylmd/thermcell_main.F90

    r4396 r4413  
    1919      USE thermcell_ini_mod, ONLY: iflag_thermals_closure,iflag_thermals_ed,tau_thermals,r_aspect_thermals
    2020      USE thermcell_ini_mod, ONLY: iflag_thermals_down
     21      USE thermcell_ini_mod, ONLY: entrdn_coef, detrdn_coef
    2122      USE thermcell_ini_mod, ONLY: RD,RG
    2223
     
    134135      real, dimension(ngrid,nlay) :: zqsat,csc
    135136      real, dimension(ngrid) :: zcon,zcon2,alim_star_tot,f
     137      real, dimension(ngrid,nlay) :: entrdn,detrdn
    136138
    137139      character (len=20) :: modname='thermcell_main'
     
    502504        print*,'AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA'
    503505        print*,'WARNING !!! routine thermcell_down en cours de developpement'
    504         ! on veut transporter la temperature potentielle, l'eau totale, qté de mouvement
    505         CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,zthl,zdthladj)
    506         CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,po,pdoadj)
    507         CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,zu,pduadj)
    508         CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,0.5*detr0,0.5*entr0,masse,zv,pdvadj)
     506        entrdn=entrdn_coef*detr0
     507        detrdn=detrdn_coef*entr0
     508        ! we want to transport potential temperature, total water and momentum
     509        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zthl,zdthladj)
     510        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,po,pdoadj)
     511        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zu,pduadj)
     512        CALL thermcell_updown_dq(ngrid,nlay,ptimestep,lmax,entr0,detr0,entrdn,detrdn,masse,zv,pdvadj)
    509513      ELSE
    510514      !--------------------------------------------------------------
Note: See TracChangeset for help on using the changeset viewer.