- Timestamp:
- Jul 24, 2024, 4:23:34 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_closure.F90
r5116 r5117 21 21 22 22 ! --- arguments ------------------------------------------ 23 integer, intent(in) :: ngrid,nlay24 real, intent(in) :: r_aspect,ptimestep25 real, intent(in), dimension(ngrid,nlay) :: alim_star,rho,zlev26 integer, intent(in), dimension(ngrid) :: lalim27 real, intent(in), dimension(ngrid) :: zmax,wmax23 INTEGER, INTENT(IN) :: ngrid,nlay 24 REAL, INTENT(IN) :: r_aspect,ptimestep 25 REAL, INTENT(IN), DIMENSION(ngrid,nlay) :: alim_star,rho,zlev 26 INTEGER, INTENT(IN), DIMENSION(ngrid) :: lalim 27 REAL, INTENT(IN), DIMENSION(ngrid) :: zmax,wmax 28 28 29 real, intent(out), dimension(ngrid) :: f29 REAL, INTENT(OUT), DIMENSION(ngrid) :: f 30 30 31 31 32 32 ! --- local ------------------------------------------ 33 real, dimension(ngrid) :: zdenom,alim_star2,alim_star_tot33 REAL, DIMENSION(ngrid) :: zdenom,alim_star2,alim_star_tot 34 34 INTEGER llmax 35 35 INTEGER ig,k … … 45 45 llmax=1 46 46 DO ig=1,ngrid 47 if(lalim(ig)>llmax) llmax=lalim(ig)47 IF (lalim(ig)>llmax) llmax=lalim(ig) 48 48 END DO 49 49 … … 53 53 DO k=1,llmax-1 54 54 do ig=1,ngrid 55 if(k<lalim(ig)) THEN55 IF (k<lalim(ig)) THEN 56 56 alim_star2(ig)=alim_star2(ig)+alim_star(ig,k)**2 & 57 57 & /(rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k))) … … 63 63 64 64 DO ig=1,ngrid 65 if(alim_star2(ig)>1.e-10) THEN65 IF (alim_star2(ig)>1.e-10) THEN 66 66 f(ig)=wmax(ig)*alim_star_tot(ig)/ & 67 67 & (max(500.,zmax(ig))*r_aspect*alim_star2(ig))
Note: See TracChangeset
for help on using the changeset viewer.