Ignore:
Timestamp:
Jul 20, 2024, 12:00:23 PM (4 months ago)
Author:
abarral
Message:

remove fixed-form \s+& remaining in .f90,.F90

File:
1 edited

Legend:

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

    r5082 r5087  
    33
    44      SUBROUTINE thermcell_height(ngrid,nlay,lalim,lmin,linter,lcong,lintercong,lmix,  &
    5      &           zw2,zlev,lmax,zmax,zmax0,zmix,wmax,zcong)
     5             zw2,zlev,lmax,zmax,zmax0,zmix,wmax,zcong)
    66      IMPLICIT NONE
    77
     
    105105! calcul de zlevinter
    106106          zlevinter(ig)=(zlev(ig,lmax(ig)+1)-zlev(ig,lmax(ig)))*  &
    107      &    linter(ig)+zlev(ig,lmax(ig))-lmax(ig)*(zlev(ig,lmax(ig)+1)  &
    108      &    -zlev(ig,lmax(ig)))
     107      linter(ig)+zlev(ig,lmax(ig))-lmax(ig)*(zlev(ig,lmax(ig)+1)  &
     108      -zlev(ig,lmax(ig)))
    109109!pour le cas ou on prend tjs lmin=1
    110110!       zmax(ig)=max(zmax(ig),zlevinter(ig)-zlev(ig,lmin(ig)))
     
    117117! calcul de zlevintercong
    118118         zlevintercong(ig)=(zlev(ig,lcong(ig)+1)-zlev(ig,lcong(ig)))*  &
    119      &    lintercong(ig)+zlev(ig,lcong(ig))-lcong(ig)*(zlev(ig,lcong(ig)+1)  &
    120      &    -zlev(ig,lcong(ig)))
     119      lintercong(ig)+zlev(ig,lcong(ig))-lcong(ig)*(zlev(ig,lcong(ig)+1)  &
     120      -zlev(ig,lcong(ig)))
    121121         zcong(ig)=zlevintercong(ig)-zlev(ig,1)
    122122!         print*,"calcul zcong",lcong(ig),lintercong(ig),zlevintercong(ig),zcong(ig)
     
    131131! test
    132132              if (((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
    133      &        *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
    134      &        -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
    135      &        *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig)))))>1e-10)  &
    136      &        then
     133          *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
     134          -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
     135          *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig)))))>1e-10)  &
     136          then
    137137!             
    138138            zmix(ig)=((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
    139      &        *((zlev(ig,lmix(ig)))**2-(zlev(ig,lmix(ig)+1))**2)  &
    140      &        -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
    141      &        *((zlev(ig,lmix(ig)-1))**2-(zlev(ig,lmix(ig)))**2))  &
    142      &        /(2.*((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
    143      &        *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
    144      &        -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
    145      &        *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))))))
     139          *((zlev(ig,lmix(ig)))**2-(zlev(ig,lmix(ig)+1))**2)  &
     140          -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
     141          *((zlev(ig,lmix(ig)-1))**2-(zlev(ig,lmix(ig)))**2))  &
     142          /(2.*((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
     143          *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
     144          -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
     145          *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))))))
    146146              else
    147147              zmix(ig)=zlev(ig,lmix(ig))
     
    162162         do l=1,nlay
    163163            if (zmix(ig)>=zlev(ig,l).and.  &
    164      &          zmix(ig)<zlev(ig,l+1)) then
     164            zmix(ig)<zlev(ig,l+1)) then
    165165              lmix(ig)=l
    166166             endif
Note: See TracChangeset for help on using the changeset viewer.