Ignore:
Timestamp:
Jul 22, 2024, 9:29:09 PM (2 months ago)
Author:
abarral
Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:
1 edited

Legend:

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

    r5082 r5099  
    11! $Id$
    2 !
     2
    33SUBROUTINE readaerosol_optic(debut, flag_aerosol, itap, rjourvrai, &
    44     pdtphys, pplay, paprs, t_seri, rhcl, presnivs, &
     
    1010! 1) recevie the aerosols(already read and interpolated) corresponding to flag_aerosol
    1111! 2) calculate the optical properties for the aerosols
    12 !
    13  
     12
    1413  USE dimphy
    1514  USE aero_mod
     
    7574!****************************************************************************************
    7675! 1) Get aerosol mass
    77 !   
     76
    7877!****************************************************************************************
    7978! Read and interpolate sulfate
     
    141140  ENDIF
    142141
    143 !
    144142! Store all aerosols in one variable
    145 !
     143
    146144  m_allaer(:,:,id_ASBCM_phy)  = bcsol(:,:)        ! ASBCM
    147145  m_allaer(:,:,id_ASPOMM_phy) = pomsol(:,:)       ! ASPOMM
     
    172170  m_allaer_pi(:,:,id_CSNO3M_phy) = 0.0
    173171  m_allaer_pi(:,:,id_CINO3M_phy) = 0.0
    174 !
     172
    175173! Calculate the total mass of all soluble aersosols
    176 !
     174
    177175  mass_solu_aero(:,:)    = sulfate(:,:)    + bcsol(:,:)    + pomsol(:,:) + ssacu(:,:)
    178176  mass_solu_aero_pi(:,:) = sulfate_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + ssacu_pi(:,:)
     
    180178!****************************************************************************************
    181179! 2) Calculate optical properties for the aerosols
    182 !
     180
    183181!****************************************************************************************
    184182  DO k = 1, klev
Note: See TracChangeset for help on using the changeset viewer.