- Timestamp:
- Jul 22, 2024, 9:29:09 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_optic.F90
r5082 r5099 1 1 ! $Id$ 2 ! 2 3 3 SUBROUTINE readaerosol_optic(debut, flag_aerosol, itap, rjourvrai, & 4 4 pdtphys, pplay, paprs, t_seri, rhcl, presnivs, & … … 10 10 ! 1) recevie the aerosols(already read and interpolated) corresponding to flag_aerosol 11 11 ! 2) calculate the optical properties for the aerosols 12 ! 13 12 14 13 USE dimphy 15 14 USE aero_mod … … 75 74 !**************************************************************************************** 76 75 ! 1) Get aerosol mass 77 ! 76 78 77 !**************************************************************************************** 79 78 ! Read and interpolate sulfate … … 141 140 ENDIF 142 141 143 !144 142 ! Store all aerosols in one variable 145 ! 143 146 144 m_allaer(:,:,id_ASBCM_phy) = bcsol(:,:) ! ASBCM 147 145 m_allaer(:,:,id_ASPOMM_phy) = pomsol(:,:) ! ASPOMM … … 172 170 m_allaer_pi(:,:,id_CSNO3M_phy) = 0.0 173 171 m_allaer_pi(:,:,id_CINO3M_phy) = 0.0 174 ! 172 175 173 ! Calculate the total mass of all soluble aersosols 176 ! 174 177 175 mass_solu_aero(:,:) = sulfate(:,:) + bcsol(:,:) + pomsol(:,:) + ssacu(:,:) 178 176 mass_solu_aero_pi(:,:) = sulfate_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + ssacu_pi(:,:) … … 180 178 !**************************************************************************************** 181 179 ! 2) Calculate optical properties for the aerosols 182 ! 180 183 181 !**************************************************************************************** 184 182 DO k = 1, klev
Note: See TracChangeset
for help on using the changeset viewer.