Ignore:
Timestamp:
Jul 20, 2024, 12:00:23 PM (12 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/cosp/lidar_simulator.F90

    r5082 r5087  
    378378!     opt. thick of each layer
    379379      tau_mol(:,1:nlev) = alpha_mol(:,1:nlev) &
    380          & *(zheight(:,2:nlev+1)-zheight(:,1:nlev))
     380   *(zheight(:,2:nlev+1)-zheight(:,1:nlev))
    381381!     opt. thick from TOA
    382382      DO k = nlev-1, 1, -1
     
    390390!       opt. thick of each layer
    391391        tau_part(:,:,i) = tau_part(:,:,i) &
    392            & * (zheight(:,2:nlev+1)-zheight(:,1:nlev) )
     392   * (zheight(:,2:nlev+1)-zheight(:,1:nlev) )
    393393!       opt. thick from TOA
    394394        DO k = nlev-1, 1, -1
     
    400400!      Upper layer
    401401       pmol(:,nlev) = beta_mol(:,nlev) / (2.*tau_mol(:,nlev)) &
    402             & * (1.-exp(-2.0*tau_mol(:,nlev)))
     402   * (1.-exp(-2.0*tau_mol(:,nlev)))
    403403!      Other layers
    404404       DO k= nlev-1, 1, -1
     
    406406        WHERE (tau_mol_lay(:)>0.)
    407407          pmol(:,k) = beta_mol(:,k) * EXP(-2.0*tau_mol(:,k+1)) / (2.*tau_mol_lay(:)) &
    408             & * (1.-exp(-2.0*tau_mol_lay(:)))
     408   * (1.-exp(-2.0*tau_mol_lay(:)))
    409409        ELSEWHERE
    410410!         This must never happend, but just in case, to avoid div. by 0
     
    429429!     Upper layer
    430430      pnorm(:,nlev) = betatot(:,nlev) / (2.*tautot(:,nlev)) &
    431             & * (1.-exp(-2.0*tautot(:,nlev)))
     431   * (1.-exp(-2.0*tautot(:,nlev)))
    432432
    433433!     Other layers
     
    436436        WHERE (tautot_lay(:)>0.)
    437437          pnorm(:,k) = betatot(:,k) * EXP(-2.0*tautot(:,k+1)) / (2.*tautot_lay(:)) &
    438                & * (1.-EXP(-2.0*tautot_lay(:)))
     438   * (1.-EXP(-2.0*tautot_lay(:)))
    439439        ELSEWHERE
    440440!         This must never happend, but just in case, to avoid div. by 0
     
    468468!     Upper layer
    469469      pnorm_ice(:,nlev) = betatot_ice(:,nlev) / (2.*tautot_ice(:,nlev)) &
    470             & * (1.-exp(-2.0*tautot_ice(:,nlev)))
     470   * (1.-exp(-2.0*tautot_ice(:,nlev)))
    471471
    472472      DO k= nlev-1, 1, -1
     
    474474        WHERE (tautot_lay_ice(:)>0.)
    475475         pnorm_ice(:,k)=betatot_ice(:,k)*EXP(-2.0*tautot_ice(:,k+1))/(2.*tautot_lay_ice(:)) &
    476                & * (1.-EXP(-2.0*tautot_lay_ice(:)))
     476   * (1.-EXP(-2.0*tautot_lay_ice(:)))
    477477        ELSEWHERE
    478478         pnorm_ice(:,k)=betatot_ice(:,k)*EXP(-2.0*tautot_ice(:,k+1))
     
    483483!     Upper layer
    484484      pnorm_liq(:,nlev) = betatot_liq(:,nlev) / (2.*tautot_liq(:,nlev)) &
    485             & * (1.-exp(-2.0*tautot_liq(:,nlev)))
     485   * (1.-exp(-2.0*tautot_liq(:,nlev)))
    486486
    487487      DO k= nlev-1, 1, -1
     
    489489        WHERE (tautot_lay_liq(:)>0.)
    490490          pnorm_liq(:,k)=betatot_liq(:,k)*EXP(-2.0*tautot_liq(:,k+1))/(2.*tautot_lay_liq(:)) &
    491                & * (1.-EXP(-2.0*tautot_lay_liq(:)))
     491   * (1.-EXP(-2.0*tautot_lay_liq(:)))
    492492        ELSEWHERE
    493493          pnorm_liq(:,k)=betatot_liq(:,k)*EXP(-2.0*tautot_liq(:,k+1))
     
    510510!     Upper layer
    511511      beta_perp_ice(:,nlev) = pnorm_perp_ice(:,nlev) * (2.*tautot_ice(:,nlev)) &
    512             & / (1.-exp(-2.0*tautot_ice(:,nlev)))
     512   / (1.-exp(-2.0*tautot_ice(:,nlev)))
    513513
    514514      DO k= nlev-1, 1, -1
     
    516516        WHERE (tautot_lay_ice(:)>0.)
    517517         beta_perp_ice(:,k) = pnorm_perp_ice(:,k)/ EXP(-2.0*tautot_ice(:,k+1)) * (2.*tautot_lay_ice(:)) &
    518             & / (1.-exp(-2.0*tautot_lay_ice(:)))
     518   / (1.-exp(-2.0*tautot_lay_ice(:)))
    519519
    520520        ELSEWHERE
     
    526526!     Upper layer
    527527      beta_perp_liq(:,nlev) = pnorm_perp_liq(:,nlev) * (2.*tautot_liq(:,nlev)) &
    528             & / (1.-exp(-2.0*tautot_liq(:,nlev)))
     528   / (1.-exp(-2.0*tautot_liq(:,nlev)))
    529529
    530530      DO k= nlev-1, 1, -1
     
    532532        WHERE (tautot_lay_liq(:)>0.)
    533533         beta_perp_liq(:,k) = pnorm_perp_liq(:,k)/ max(seuil,EXP(-2.0*tautot_liq(:,k+1))) &
    534             & * (2.*tautot_lay_liq(:)) / (1.-exp(-2.0*tautot_lay_liq(:)))
     534   * (2.*tautot_lay_liq(:)) / (1.-exp(-2.0*tautot_lay_liq(:)))
    535535
    536536        ELSEWHERE
     
    550550          pnorm_perp_tot(:,nlev) = &
    551551              (beta_perp_ice(:,nlev)+beta_perp_liq(:,nlev)-(beta_mol(:,nlev)/(1+1/0.0284))) / (2.*tautot(:,nlev)) &
    552               & * (1.-exp(-2.0*tautot(:,nlev)))
     552   * (1.-exp(-2.0*tautot(:,nlev)))
    553553    ELSEWHERE
    554554    pnorm_perp_tot(:,nlev) = 0.
     
    570570                          (beta_perp_ice(:,k)+beta_perp_liq(:,k)-(beta_mol(:,k)/(1+1/0.0284))) * &
    571571                          EXP(-2.0*tautot(:,k+1)) / (2.*tautot_lay(:)) &
    572                           & * (1.-EXP(-2.0*tautot_lay(:)))
     572   * (1.-EXP(-2.0*tautot_lay(:)))
    573573                    ELSEWHERE
    574574          !         This must never happen, but just in case, to avoid div. by 0
Note: See TracChangeset for help on using the changeset viewer.