Changeset 2738 for LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

Timestamp:

Jan 2, 2017, 11:24:30 PM (8 years ago)

Author:

oboucher

Message:

Swapping the order of CSSO4 and ASSO4 aerosols (and fixing an issue on aerindex). Preparing the ground for nitrate aerosols (coarse soluble, accumulation soluble, coarse insoluble). Modifying the LW aeropt routine so that it is compatible with both INCA and climatological aerosols (for dust only). Adding a new flag ok_alw for activating aerosol direct LW effect (for dust only). This change is bit comparable for flag_aerosol=6, flag_rrtm=1, NSW=6.

File:

• LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

@@ -1 +1 @@
 ! $Id$
 !
-SUBROUTINE readaerosol_optic_rrtm(debut, aerosol_couple,  &
+SUBROUTINE readaerosol_optic_rrtm(debut, aerosol_couple, ok_alw,  &
      new_aod, flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, &
      pdtphys, pplay, paprs, t_seri, rhcl, presnivs, &
@@ -30 +30 @@
   LOGICAL, INTENT(IN)                      :: debut
   LOGICAL, INTENT(IN)                      :: aerosol_couple
+  LOGICAL, INTENT(IN)                      :: ok_alw
   LOGICAL, INTENT(IN)                      :: new_aod
   INTEGER, INTENT(IN)                      :: flag_aerosol
@@ -281 +282 @@
   ! Calculate the total mass of all soluble aersosols
   ! to be revisited for AR6
-  mass_solu_aero(:,:)    = sulfacc(:,:)    + bcsol(:,:)    + pomsol(:,:)   + nitracc(:,:) + ssacu(:,:)
-  mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:) 
+  mass_solu_aero(:,:)    = sulfacc(:,:)    + bcsol(:,:)    + pomsol(:,:)    + nitracc(:,:)    + ssacu(:,:)
+  mass_solu_aero_pi(:,:) = sulfacc_pi(:,:) + bcsol_pi(:,:) + pomsol_pi(:,:) + nitracc_pi(:,:) + ssacu_pi(:,:)
 
   !****************************************************************************************
@@ -310 +311 @@
        pplay, t_seri,                &
        tausum_aero, tau3d_aero )
+
+  !--call LW optical properties for tropospheric aerosols 
+  CALL aeropt_lw_rrtm(ok_alw, pdel, zrho, flag_aerosol, m_allaer, m_allaer_pi)
 
   ! Diagnostics calculation for CMIP5 protocol