Changeset 4413
- Timestamp:
- Feb 1, 2023, 5:47:09 PM (22 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/thermcell_ini_mod.F90
r4381 r4413 8 8 integer :: iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure 9 9 integer :: iflag_thermals_down 10 real :: entrdn_coef, detrdn_coef 10 11 11 12 !$OMP THREADPRIVATE(dvdq,dqimpl,prt_level,lunout) … … 14 15 !$OMP THREADPRIVATE(iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure) 15 16 !$OMP THREADPRIVATE(iflag_thermals_down) 17 !$OMP THREADPRIVATE(entrdn_coef, detrdn_coef) 16 18 17 19 … … 91 93 92 94 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 93 102 CALL getin_p('thermals_fact_epsilon',fact_epsilon) 94 103 CALL getin_p('thermals_betalpha',betalpha) -
LMDZ6/trunk/libf/phylmd/thermcell_main.F90
r4396 r4413 19 19 USE thermcell_ini_mod, ONLY: iflag_thermals_closure,iflag_thermals_ed,tau_thermals,r_aspect_thermals 20 20 USE thermcell_ini_mod, ONLY: iflag_thermals_down 21 USE thermcell_ini_mod, ONLY: entrdn_coef, detrdn_coef 21 22 USE thermcell_ini_mod, ONLY: RD,RG 22 23 … … 134 135 real, dimension(ngrid,nlay) :: zqsat,csc 135 136 real, dimension(ngrid) :: zcon,zcon2,alim_star_tot,f 137 real, dimension(ngrid,nlay) :: entrdn,detrdn 136 138 137 139 character (len=20) :: modname='thermcell_main' … … 502 504 print*,'AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA' 503 505 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) 509 513 ELSE 510 514 !--------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.