- Timestamp:
- Jun 17, 2022, 4:24:49 PM (2 years ago)
- Location:
- LMDZ6/branches/LMDZ-ECRAD
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-ECRAD
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/thermcell_closure.F90
r2311 r4171 3 3 ! 4 4 SUBROUTINE thermcell_closure(ngrid,nlay,r_aspect,ptimestep,rho, & 5 & zlev,lalim,alim_star, f_star,zmax,wmax,f,lev_out)5 & zlev,lalim,alim_star,zmax,wmax,f) 6 6 7 7 !------------------------------------------------------------------------- … … 17 17 IMPLICIT NONE 18 18 19 #include "thermcell.h" 20 INTEGER ngrid,nlay 19 ! --- arguments ------------------------------------------ 20 integer, intent(in) :: ngrid,nlay 21 real, intent(in) :: r_aspect,ptimestep 22 real, intent(in), dimension(ngrid,nlay) :: alim_star,rho,zlev 23 integer, intent(in), dimension(ngrid) :: lalim 24 real, intent(in), dimension(ngrid) :: zmax,wmax 25 26 real, intent(out), dimension(ngrid) :: f 27 28 29 ! --- local ------------------------------------------ 30 real, dimension(ngrid) :: zdenom,alim_star2,alim_star_tot 31 INTEGER llmax 21 32 INTEGER ig,k 22 REAL r_aspect,ptimestep23 integer lev_out ! niveau pour les print24 25 INTEGER lalim(ngrid)26 REAL alim_star(ngrid,nlay)27 REAL f_star(ngrid,nlay+1)28 REAL rho(ngrid,nlay)29 REAL zlev(ngrid,nlay)30 REAL zmax(ngrid)31 REAL wmax(ngrid)32 REAL zdenom(ngrid)33 REAL alim_star2(ngrid)34 REAL f(ngrid)35 36 REAL alim_star_tot(ngrid)37 INTEGER llmax38 33 39 34 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 72 67 73 68 74 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!75 ! TESTS POUR UNE NOUVELLE FERMETURE DANS LAQUELLE ALIM_STAR NE SERAIT76 ! PAS NORMALISE77 ! f(ig)=f(ig)*f_star(ig,2)/(f_star(ig,lalim(ig)))78 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!79 69 80 return70 RETURN 81 71 end
Note: See TracChangeset
for help on using the changeset viewer.