Ignore:
Timestamp:
Jul 24, 2024, 4:23:34 PM (4 months ago)
Author:
abarral
Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_closure.F90

    r5116 r5117  
    2121
    2222! --- arguments ------------------------------------------
    23 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
     23INTEGER, INTENT(IN) :: ngrid,nlay
     24REAL, INTENT(IN) :: r_aspect,ptimestep
     25REAL, INTENT(IN), DIMENSION(ngrid,nlay) :: alim_star,rho,zlev
     26INTEGER, INTENT(IN), DIMENSION(ngrid) :: lalim
     27REAL, INTENT(IN), DIMENSION(ngrid) :: zmax,wmax
    2828
    29 real, intent(out), dimension(ngrid) :: f
     29REAL, INTENT(OUT), DIMENSION(ngrid) :: f
    3030
    3131
    3232! --- local ------------------------------------------
    33 real, dimension(ngrid) :: zdenom,alim_star2,alim_star_tot
     33REAL, DIMENSION(ngrid) :: zdenom,alim_star2,alim_star_tot
    3434INTEGER llmax
    3535INTEGER ig,k       
     
    4545llmax=1
    4646DO ig=1,ngrid
    47    if (lalim(ig)>llmax) llmax=lalim(ig)
     47   IF (lalim(ig)>llmax) llmax=lalim(ig)
    4848END DO
    4949
     
    5353DO k=1,llmax-1
    5454   do ig=1,ngrid
    55       if (k<lalim(ig)) THEN
     55      IF (k<lalim(ig)) THEN
    5656         alim_star2(ig)=alim_star2(ig)+alim_star(ig,k)**2  &
    5757&                    /(rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k)))
     
    6363
    6464DO ig=1,ngrid
    65    if (alim_star2(ig)>1.e-10) THEN
     65   IF (alim_star2(ig)>1.e-10) THEN
    6666      f(ig)=wmax(ig)*alim_star_tot(ig)/  &
    6767&     (max(500.,zmax(ig))*r_aspect*alim_star2(ig))
Note: See TracChangeset for help on using the changeset viewer.