Index: LMDZ5/trunk/libf/phylmd/conf_phys_m.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/conf_phys_m.F90	(revision 2737)
+++ LMDZ5/trunk/libf/phylmd/conf_phys_m.F90	(revision 2738)
@@ -17,5 +17,5 @@
        iflag_cld_th, &
        iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
-       ok_ade, ok_aie, ok_cdnc, aerosol_couple, &
+       ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple, &
        flag_aerosol, flag_aerosol_strat, new_aod, &
        flag_bc_internal_mixture, bl95_b0, bl95_b1,&
@@ -61,4 +61,5 @@
     ! ok_instan:  sorties instantanees
     ! ok_ade, ok_aie: apply or not aerosol direct and indirect effects
+    ! ok_alw: activate aerosol LW effect
     ! ok_cdnc, ok cloud droplet number concentration
     ! flag_aerosol_strat : flag pour les aerosols stratos
@@ -74,5 +75,5 @@
     LOGICAL              :: ok_LES
     LOGICAL              :: callstats
-    LOGICAL              :: ok_ade, ok_aie, ok_cdnc, aerosol_couple
+    LOGICAL              :: ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple
     INTEGER              :: flag_aerosol
     INTEGER              :: flag_aerosol_strat
@@ -93,5 +94,5 @@
     LOGICAL, SAVE       :: ok_LES_omp   
     LOGICAL, SAVE       :: callstats_omp
-    LOGICAL, SAVE       :: ok_ade_omp, ok_aie_omp, ok_cdnc_omp, aerosol_couple_omp
+    LOGICAL, SAVE       :: ok_ade_omp, ok_aie_omp, ok_alw_omp, ok_cdnc_omp, aerosol_couple_omp
     INTEGER, SAVE       :: flag_aerosol_omp
     INTEGER, SAVE       :: flag_aerosol_strat_omp
@@ -347,4 +348,12 @@
     ok_ade_omp = .FALSE.
     CALL getin('ok_ade', ok_ade_omp)
+
+    !Config Key  = ok_alw
+    !Config Desc = Aerosol longwave effect or not?
+    !Config Def  = .FALSE.
+    !Config Help = Used in radlwsw.F
+    !
+    ok_alw_omp = .FALSE.
+    CALL getin('ok_alw', ok_alw_omp)
 
     !
@@ -2161,4 +2170,5 @@
     ok_ade = ok_ade_omp
     ok_aie = ok_aie_omp
+    ok_alw = ok_alw_omp
     ok_cdnc = ok_cdnc_omp
     aerosol_couple = aerosol_couple_omp
@@ -2468,4 +2478,5 @@
     write(lunout,*)' ok_ade = ',ok_ade
     write(lunout,*)' ok_aie = ',ok_aie
+    write(lunout,*)' ok_alw = ',ok_alw
     write(lunout,*)' aerosol_couple = ', aerosol_couple
     write(lunout,*)' flag_aerosol = ', flag_aerosol
Index: LMDZ5/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2737)
+++ LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2738)
@@ -980,8 +980,9 @@
     ! Parameters
     LOGICAL ok_ade, ok_aie    ! Apply aerosol (in)direct effects or not
+    LOGICAL ok_alw            ! Apply aerosol LW effect or not
     LOGICAL ok_cdnc ! ok cloud droplet number concentration (O. Boucher 01-2013)
     REAL bl95_b0, bl95_b1   ! Parameter in Boucher and Lohmann (1995)
-    SAVE ok_ade, ok_aie, ok_cdnc, bl95_b0, bl95_b1
-    !$OMP THREADPRIVATE(ok_ade, ok_aie, ok_cdnc, bl95_b0, bl95_b1)
+    SAVE ok_ade, ok_aie, ok_alw, ok_cdnc, bl95_b0, bl95_b1
+    !$OMP THREADPRIVATE(ok_ade, ok_aie, ok_alw, ok_cdnc, bl95_b0, bl95_b1)
     LOGICAL, SAVE :: aerosol_couple ! true  : calcul des aerosols dans INCA
     ! false : lecture des aerosol dans un fichier
@@ -1142,5 +1143,5 @@
             fact_cldcon, facttemps,ok_newmicro,iflag_radia, &
             iflag_cld_th,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
-            ok_ade, ok_aie, ok_cdnc, aerosol_couple,  &
+            ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple,  &
             flag_aerosol, flag_aerosol_strat, new_aod, &
             flag_bc_internal_mixture, bl95_b0, bl95_b1, &
@@ -1587,5 +1588,5 @@
             jjmp1,nlevSTD,clevSTD,rlevSTD, dtime,ok_veget, &
             type_ocean,iflag_pbl,iflag_pbl_split,ok_mensuel,ok_journe, &
-            ok_hf,ok_instan,ok_LES,ok_ade,ok_aie,  &
+            ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, &
             read_climoz, phys_out_filestations, &
             new_aod, aerosol_couple, &
@@ -3293,5 +3294,5 @@
        !jq - introduce the aerosol direct and first indirect radiative forcings
        !jq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
-       IF (flag_aerosol .gt. 0) THEN
+       IF (flag_aerosol .GT. 0) THEN
           IF (iflag_rrtm .EQ. 0) THEN !--old radiation
              IF (.NOT. aerosol_couple) THEN
@@ -3312,7 +3313,7 @@
 #ifdef CPP_RRTM
                 IF (NSW.EQ.6) THEN
-                   !--new aerosol properties
+                   !--new aerosol properties SW and LW
                    !
-                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, &
+                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, ok_alw, &
                         new_aod, flag_aerosol, flag_bc_internal_mixture, itap, jD_cur-jD_ref, &
                         pdtphys, pplay, paprs, t_seri, rhcl, presnivs,  &
@@ -3320,5 +3321,5 @@
                         tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,  &
                         tausum_aero, tau3d_aero)
-
+                   !
                 ELSE IF (NSW.EQ.2) THEN 
                    !--for now we use the old aerosol properties
@@ -3339,4 +3340,8 @@
                    piz_aero_sw_rrtm(:,:,2,:)=piz_aero(:,:,2,:)
                    cg_aero_sw_rrtm (:,:,2,:)=cg_aero (:,:,2,:)
+                   !
+                   !--no LW optics
+                   tau_aero_lw_rrtm(:,:,:,:) = 1.e-15
+                   !
                 ELSE
                    abort_message='Only NSW=2 or 6 are possible with ' &
@@ -3344,9 +3349,4 @@
                    CALL abort_physic(modname,abort_message,1)
                 ENDIF
-
-                !--call LW optical properties for tropospheric aerosols 
-                !--only works for INCA aerosol (aerosol_couple = TRUE)
-                CALL aeropt_lw_rrtm(aerosol_couple,paprs,tr_seri)
-                !
 #else
                 abort_message='You should compile with -rrtm if running ' &
@@ -3357,5 +3357,5 @@
              ENDIF
           ENDIF
-       ELSE
+       ELSE   !--flag_aerosol = 0 
           tausum_aero(:,:,:) = 0.
           mass_solu_aero(:,:) = 0.
Index: LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90	(revision 2737)
+++ LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90	(revision 2738)
@@ -83,7 +83,8 @@
   INTEGER, PARAMETER :: la865 = 5
   INTEGER, PARAMETER :: nbre_RH=12
-  INTEGER, PARAMETER :: naero_soluble=7   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.
-                                          !  4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
-  INTEGER, PARAMETER :: naero_insoluble=3 !  1- Dust; 2- BC insoluble; 3- POM insoluble
+  INTEGER, PARAMETER :: naero_soluble=9   !  1- BC soluble; 2- POM soluble; 3- SO4 coarse
+                                          !  4- SO4 acc; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
+                                          !  8- NO3 coarse; 9- NO3 accumulation
+  INTEGER, PARAMETER :: naero_insoluble=4 !  1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
 
   REAL :: zrho
@@ -101,7 +102,5 @@
   REAL :: zdh(klon,klev)
   
-   ! Soluble components 1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-coarse; 6 seasalt coarse; 7 seasalt acc.
   REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! Ext. coeff. ** m2/g 
-   ! Insoluble components 1- Dust: 2- BC; 3- POM
   REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! Ext. coeff. ** m2/g 
 
@@ -110,6 +109,6 @@
   !
   INTEGER, PARAMETER ::  nbclassbc = 5  ! Added by Rong Wang/OB for the 5 fractions
-                                       ! of BC in the soluble mode:
-                                       ! bc_content/0.001, 0.01, 0.02, 0.05, ! 0.1/
+                                        ! of BC in the soluble mode:
+                                        ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
   ! for Maxwell-Garnet internal mixture
   ! Detailed theory can be found in R. Wang Estimation of global black carbon ! direct
@@ -124,8 +123,6 @@
   INTEGER :: n, classbc
 
-! From here on we look at the optical parameters at 5 wavelengths:  
-! 443nm, 550, 670, 765 and 865 nm 
-!                                   le 12 AVRIL 2006 
-!  
+! From here on we look at the optical parameters at 5 wavelengths 443, 550, 670, 765 and 865 nm 
+
  DATA alpha_aers_5wv/ & 
                                 ! bc soluble 
@@ -170,13 +167,27 @@
         2.340, 3.358, 3.803, 4.303, 4.928, 4.696, 5.629, 7.198,10.308,13.342,19.120,34.296, &
         1.789, 2.626, 2.999, 3.422, 3.955, 3.730, 4.519, 5.864, 8.593,11.319,16.653,31.331, &
-        1.359, 2.037, 2.343, 2.693, 3.139, 2.940, 3.596, 4.729, 7.076, 9.469,14.266,28.043 /
+        1.359, 2.037, 2.343, 2.693, 3.139, 2.940, 3.596, 4.729, 7.076, 9.469,14.266,28.043, &
+                           ! NO3 coarse (CS)
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+                           ! NO3 accumulation (AS)
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
 
   DATA alpha_aeri_5wv/ &
                                  ! dust insoluble 
-        0.759, 0.770, 0.775, 0.775, 0.772, & 
+        0.759, 0.770, 0.775, 0.775, 0.772,  & 
                                  !!jb bc insoluble 
         11.536,10.033, 8.422, 7.234, 6.270, & 
                                  ! pom insoluble 
-        5.042, 3.101, 1.890, 1.294, 0.934/ 
+        5.042, 3.101, 1.890, 1.294, 0.934,  &
+                                 ! NO3 insoluble 
+        0.0, 0.0, 0.0, 0.0, 0.0 /
 
   ! internal mixture
@@ -252,5 +263,5 @@
      aerosol_name(1) = id_CIDUSTM_phy
   ELSEIF (flag_aerosol .EQ. 6) THEN 
-     nb_aer = 10
+     nb_aer = 13
      ALLOCATE (aerosol_name(nb_aer)) 
      aerosol_name(1) = id_ASSO4M_phy      
@@ -264,4 +275,7 @@
      aerosol_name(9) = id_CIDUSTM_phy
      aerosol_name(10)= id_CSSO4M_phy
+     aerosol_name(11)= id_CSNO3M_phy
+     aerosol_name(12)= id_ASNO3M_phy
+     aerosol_name(13)= id_CINO3M_phy
   ENDIF
 
@@ -315,4 +329,12 @@
         soluble=.TRUE.
         spsol=7
+    ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
+        soluble=.TRUE.
+        spsol=8
+        fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
+    ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
+        soluble=.TRUE.
+        spsol=9
+        fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
     ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN 
         soluble=.FALSE.
@@ -324,4 +346,8 @@
         soluble=.FALSE.
         spinsol=3
+    ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN 
+        soluble=.FALSE.
+        spinsol=4
+        fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
     ELSE 
         CYCLE
Index: LMDZ5/trunk/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90	(revision 2737)
+++ LMDZ5/trunk/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90	(revision 2738)
@@ -53,7 +53,9 @@
   INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
 
-  INTEGER, PARAMETER :: naero_soluble=7    ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
-                                           ! 5- seasalt super coarse  6- seasalt coarse   7- seasalt acc.
-  INTEGER, PARAMETER :: naero_insoluble=3  ! 1- Dust; 2- BC insoluble; 3- POM insoluble
+  INTEGER, PARAMETER :: naero_soluble=9    ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
+                                           ! 5- seasalt super coarse  6- seasalt
+                                           ! coarse  7- seasalt acc.
+                                           ! 8- NO3 coarse 9- NO3 acc.
+  INTEGER, PARAMETER :: naero_insoluble=4  ! 1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
 
   INTEGER, PARAMETER :: nbre_RH=12
@@ -80,16 +82,16 @@
   ! Proprietes optiques
   !
-  REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)   !--unit m2/g SO4
-  REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble)
-  REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)      !--unit 
-  REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble)
-  REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)     !-- unit
-  REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble)        !-- unit
+  REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)   !--unit m2/g
+  REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble)         !--unit m2/g
+  REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)      !--unitless
+  REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble)            !--unitless
+  REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)     !--unitless
+  REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble)           !--unitless
   ! 
   ! BC internal mixture
   !
   INTEGER, PARAMETER ::  nbclassbc = 5  ! Added by Rong Wang/OB for the 5 fractions
-                                       ! of BC in the soluble mode:
-                                       ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
+                                        ! of BC in the soluble mode:
+                                        ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
   ! for Maxwell-Garnet internal mixture
   ! Detailed theory can be found in R. Wang Estimation of global black carbon direct
@@ -100,5 +102,5 @@
   REAL :: piz_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
   !
-  INTEGER :: id, classbc, classbc_pi
+  INTEGER :: aerindex, classbc, classbc_pi
   REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi
 
@@ -126,4 +128,11 @@
   0.377, 0.377, 0.377, 0.377, 0.377, 0.415, 0.456, 0.542, 0.627, 0.636, 0.646, 0.738, &
   0.052, 0.052, 0.052, 0.052, 0.052, 0.057, 0.063, 0.075, 0.087, 0.088, 0.089, 0.102, &
+       ! sulfate coarse CS
+  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
+  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
+  0.715, 0.715, 0.715, 0.788, 0.865, 0.945, 1.029, 1.171, 1.557, 1.936, 1.936, 3.459, &
+  0.736, 0.736, 0.736, 0.811, 0.891, 0.973, 1.060, 1.206, 1.603, 1.994, 1.994, 3.563, &
+  0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, &
+  0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, &
        ! sulfate AS   
   6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
@@ -133,11 +142,4 @@
   0.312, 0.312, 0.312, 0.344, 0.378, 0.413, 0.450, 0.511, 0.680, 0.846, 0.846, 1.511, &
   0.121, 0.121, 0.121, 0.134, 0.147, 0.161, 0.175, 0.199, 0.264, 0.329, 0.329, 0.588, &
-       ! sulfate coarse CS
-  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
-  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
-  0.715, 0.715, 0.715, 0.788, 0.865, 0.945, 1.029, 1.171, 1.557, 1.936, 1.936, 3.459, &
-  0.736, 0.736, 0.736, 0.811, 0.891, 0.973, 1.060, 1.206, 1.603, 1.994, 1.994, 3.563, &
-  0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, &
-  0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, &
         ! seasalt seasalt Super Coarse Soluble (SS)
   0.214, 0.267, 0.287, 0.310, 0.337, 0.373, 0.421, 0.494, 0.625, 0.742, 0.956, 1.480, &
@@ -160,5 +162,19 @@
   1.389, 2.066, 2.371, 2.719, 3.161, 2.970, 3.623, 4.747, 7.064, 9.419,14.120,27.565, &
   0.309, 0.497, 0.587, 0.693, 0.832, 0.777, 0.986, 1.364, 2.207, 3.136, 5.174,12.065, &
-  0.021, 0.037, 0.045, 0.054, 0.067, 0.065, 0.087, 0.129, 0.230, 0.353, 0.654, 1.885  /
+  0.021, 0.037, 0.045, 0.054, 0.067, 0.065, 0.087, 0.129, 0.230, 0.353, 0.654, 1.885, & 
+        ! NO3 coarse (CS)
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        ! NO3 accumulation (AS)
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
 
   DATA alpha_aeri_6bands/  & 
@@ -168,5 +184,7 @@
   6.497, 6.497, 5.900, 4.284, 2.163, 0.966, &
        ! pom insoluble AI
-  6.443, 6.443, 4.381, 1.846, 0.377, 0.052 /
+  6.443, 6.443, 4.381, 1.846, 0.377, 0.052, &
+       ! NO3 insoluble CI
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
 
   DATA cg_aers_6bands/ &
@@ -185,4 +203,11 @@
   0.363, 0.363, 0.363, 0.363, 0.363, 0.360, 0.357, 0.352, 0.350, 0.349, 0.349, 0.347, &
   0.142, 0.142, 0.142, 0.142, 0.142, 0.139, 0.137, 0.133, 0.131, 0.131, 0.130, 0.129, &
+       ! sulfate coarse CS
+  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
+  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
+  0.799, 0.799, 0.799, 0.787, 0.777, 0.768, 0.760, 0.747, 0.736, 0.729, 0.725, 0.716, &
+  0.797, 0.797, 0.797, 0.782, 0.770, 0.760, 0.750, 0.735, 0.722, 0.714, 0.709, 0.698, &
+  0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, &
+  0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, &
        ! sulfate AS 
   0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
@@ -192,11 +217,4 @@
   0.362, 0.362, 0.362, 0.359, 0.356, 0.354, 0.352, 0.348, 0.345, 0.343, 0.342, 0.340, &
   0.137, 0.137, 0.137, 0.135, 0.133, 0.132, 0.130, 0.128, 0.126, 0.125, 0.124, 0.122, &
-       ! sulfate coarse CS
-  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
-  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
-  0.799, 0.799, 0.799, 0.787, 0.777, 0.768, 0.760, 0.747, 0.736, 0.729, 0.725, 0.716, &
-  0.797, 0.797, 0.797, 0.782, 0.770, 0.760, 0.750, 0.735, 0.722, 0.714, 0.709, 0.698, &
-  0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, &
-  0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, &
         ! seasalt seasalt Super Coarse Soluble (SS)
   0.797, 0.800, 0.801, 0.802, 0.804, 0.822, 0.825, 0.828, 0.832, 0.835, 0.838, 0.843, &
@@ -219,5 +237,19 @@
   0.561, 0.585, 0.593, 0.601, 0.609, 0.637, 0.651, 0.669, 0.691, 0.705, 0.723, 0.745, &
   0.392, 0.427, 0.439, 0.451, 0.464, 0.480, 0.500, 0.526, 0.563, 0.588, 0.621, 0.671, &
-  0.144, 0.170, 0.179, 0.189, 0.201, 0.207, 0.224, 0.248, 0.285, 0.315, 0.359, 0.439 /
+  0.144, 0.170, 0.179, 0.189, 0.201, 0.207, 0.224, 0.248, 0.285, 0.315, 0.359, 0.439, &
+        ! NO3 coarse (CS)
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+        ! NO3 accumulation (AS)
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0   /
 
   DATA cg_aeri_6bands/ &
@@ -227,5 +259,7 @@
   0.721, 0.721, 0.643, 0.513, 0.321, 0.153, &
        ! pom insoluble AI
-  0.687, 0.687, 0.658, 0.564, 0.363, 0.142  /
+  0.687, 0.687, 0.658, 0.564, 0.363, 0.142, &
+       ! NO3 insoluble CI
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0  /
 
   DATA piz_aers_6bands/&
@@ -244,4 +278,11 @@
   0.816, 0.816, 0.816, 0.816, 0.816, 0.839, 0.859, 0.888, 0.908, 0.910, 0.911, 0.925, &
   0.663, 0.663, 0.663, 0.663, 0.663, 0.607, 0.562, 0.492, 0.446, 0.441, 0.437, 0.404, &
+       ! sulfate coarse CS
+  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
+  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
+  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
+  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
+  0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
+  0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, &
        ! sulfate AS
   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
@@ -251,11 +292,4 @@
   0.988, 0.988, 0.988, 0.989, 0.990, 0.990, 0.991, 0.992, 0.993, 0.993, 0.994, 0.994, &
   0.256, 0.256, 0.256, 0.263, 0.268, 0.273, 0.277, 0.284, 0.290, 0.294, 0.296, 0.301, &
-       ! sulfate coarse CS
-  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
-  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
-  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
-  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
-  0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
-  0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, &
         ! seasalt seasalt Super Coarse Soluble (SS)
   0.983, 0.982, 0.982, 0.982, 0.982, 0.992, 0.993, 0.994, 0.995, 0.996, 0.997, 0.998, &
@@ -278,5 +312,19 @@
   0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
   0.997, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, &
-  0.985, 0.989, 0.990, 0.990, 0.991, 0.996, 0.997, 0.998, 0.999, 0.999, 0.999, 1.000 /
+  0.985, 0.989, 0.990, 0.990, 0.991, 0.996, 0.997, 0.998, 0.999, 0.999, 0.999, 1.000, &
+        ! NO3 coarse (CS)
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+        ! NO3 accumulation (AS)
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0   /
 
   DATA piz_aeri_6bands/ &
@@ -286,5 +334,7 @@
   0.460, 0.460, 0.445, 0.394, 0.267, 0.121, &
        ! pom insoluble AI
-  0.973, 0.973, 0.972, 0.940, 0.816, 0.663 /
+  0.973, 0.973, 0.972, 0.940, 0.816, 0.663, &
+       ! NO3 insoluble CI
+  1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /
 
 ! Added by R. Wang (July 31 2016)
@@ -440,5 +490,5 @@
      aerosol_name(1) = id_CIDUSTM_phy
   ELSEIF (flag_aerosol .EQ. 6) THEN 
-     nb_aer = 10
+     nb_aer = 13
      ALLOCATE (aerosol_name(nb_aer)) 
      aerosol_name(1) = id_ASSO4M_phy      
@@ -452,4 +502,7 @@
      aerosol_name(9) = id_CIDUSTM_phy
      aerosol_name(10)= id_CSSO4M_phy
+     aerosol_name(11)= id_CSNO3M_phy
+     aerosol_name(12)= id_ASNO3M_phy
+     aerosol_name(13)= id_CINO3M_phy
   ENDIF
 
@@ -491,9 +544,9 @@
         soluble=.TRUE.
         spsol=2 
-     ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN 
+     ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN 
         soluble=.TRUE.
         spsol=3
         fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
-     ELSEIF  (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
+     ELSEIF  (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
         soluble=.TRUE.
         spsol=4
@@ -508,4 +561,12 @@
          soluble=.TRUE.
          spsol=7
+     ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
+         soluble=.TRUE.
+         spsol=8
+         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
+     ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
+         soluble=.TRUE.
+         spsol=9
+         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
      ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN 
          soluble=.FALSE.
@@ -517,4 +578,8 @@
          soluble=.FALSE.
          spinsol=3
+     ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN
+         soluble=.FALSE.
+         spinsol=4
+         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
      ELSE 
          CYCLE
@@ -522,5 +587,5 @@
 
     !--shortname for aerosol index
-    id=aerosol_name(m)
+    aerindex=aerosol_name(m)
 
     IF (soluble) THEN
@@ -532,6 +597,6 @@
            DO i=1, klon
 
-             tmp_var=m_allaer(i,k,spsol)/1.e6*zdh(i,k)*fac
-             tmp_var_pi=m_allaer_pi(i,k,spsol)/1.e6*zdh(i,k)*fac
+             tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
+             tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
 
              ! Calculate the dry BC/(BC+SUL) mass ratio for all (natural+anthropogenic) aerosols
@@ -580,7 +645,7 @@
                              cg_MG_6bands(RH_num(i,k),inu,classbc))
 
-               tau_ae(i,k,id,inu)    = tmp_var*tau_ae2b_int
-               piz_ae(i,k,id,inu)    = piz_ae2b_int
-               cg_ae(i,k,id,inu)     = cg_ae2b_int
+               tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
+               piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
+               cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
 
                !--natural aerosols
@@ -597,7 +662,7 @@
                              cg_MG_6bands(RH_num(i,k),inu,classbc_pi))
 
-               tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
-               piz_ae_pi(i,k,id,inu) = piz_ae2b_int
-               cg_ae_pi(i,k,id,inu)  = cg_ae2b_int
+               tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
+               piz_ae_pi(i,k,aerindex,inu) = piz_ae2b_int
+               cg_ae_pi(i,k,aerindex,inu)  = cg_ae2b_int
                         
              ENDDO
@@ -610,6 +675,6 @@
          DO k=1, klev
            DO i=1, klon
-             tmp_var=m_allaer(i,k,spsol)/1.e6*zdh(i,k)*fac
-             tmp_var_pi=m_allaer_pi(i,k,spsol)/1.e6*zdh(i,k)*fac
+             tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
+             tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
 
              DO inu=1,NSW
@@ -627,8 +692,8 @@
                              cg_aers_6bands(RH_num(i,k),inu,spsol))
 
-               tau_ae(i,k,id,inu)    = tmp_var*tau_ae2b_int
-               tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
-               piz_ae(i,k,id,inu)    = piz_ae2b_int
-               cg_ae(i,k,id,inu)     = cg_ae2b_int
+               tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
+               tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
+               piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
+               cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
                         
              ENDDO
@@ -638,6 +703,6 @@
          !--external mixture case for soluble BC
          IF (spsol.EQ.1) THEN
-           piz_ae_pi(:,:,id,:) = piz_ae(:,:,id,:)
-           cg_ae_pi(:,:,id,:)  = cg_ae(:,:,id,:)
+           piz_ae_pi(:,:,aerindex,:) = piz_ae(:,:,aerindex,:)
+           cg_ae_pi(:,:,aerindex,:)  = cg_ae(:,:,aerindex,:)
          ENDIF
 
@@ -648,6 +713,6 @@
        DO k=1, klev
          DO i=1, klon
-           tmp_var=m_allaer(i,k,naero_soluble+spinsol)/1.e6*zdh(i,k)*fac
-           tmp_var_pi=m_allaer_pi(i,k,naero_soluble+spinsol)/1.e6*zdh(i,k)*fac
+           tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
+           tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
 
            DO inu=1,NSW
@@ -656,8 +721,8 @@
              cg_ae2b_int = cg_aeri_6bands(inu,spinsol) 
 
-             tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
-             tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int
-             piz_ae(i,k,id,inu) = piz_ae2b_int
-             cg_ae(i,k,id,inu)= cg_ae2b_int
+             tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
+             tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi*tau_ae2b_int
+             piz_ae(i,k,aerindex,inu) = piz_ae2b_int
+             cg_ae(i,k,aerindex,inu)= cg_ae2b_int
            ENDDO
          ENDDO
Index: LMDZ5/trunk/libf/phylmd/rrtm/aeropt_lw_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/aeropt_lw_rrtm.F90	(revision 2737)
+++ LMDZ5/trunk/libf/phylmd/rrtm/aeropt_lw_rrtm.F90	(revision 2738)
@@ -2,105 +2,72 @@
 ! aeropt_lw_rrtm.F90 2014-05-13 C. Kleinschmitt
 !                    2016-05-03 O. Boucher 
+!                    2016-12-17 O. Boucher
 !
 ! This routine feeds aerosol LW properties to RRTM
 ! we only consider absorption (not scattering)
+! we only consider dust for now
 
-SUBROUTINE AEROPT_LW_RRTM(aerosol_couple,paprs,tr_seri)
+SUBROUTINE AEROPT_LW_RRTM(ok_alw, pdel, zrho, flag_aerosol, m_allaer, m_allaer_pi)
 
   USE dimphy
   USE aero_mod
-  USE infotrac_phy
   USE phys_state_var_mod, ONLY : tau_aero_lw_rrtm
   USE YOERAD, ONLY : NLW
+  USE YOMCST
 
   IMPLICIT NONE
 
-  INCLUDE "YOMCST.h"
   INCLUDE "clesphys.h"
   !
   ! Input arguments:
   !
-  LOGICAL, INTENT(IN)                         :: aerosol_couple
-  REAL, DIMENSION(klon,klev+1), INTENT(in)    :: paprs
-  REAL, DIMENSION(klon,klev,nbtr), INTENT(in) :: tr_seri
-  !
-  REAL, DIMENSION(klon,klev) :: zdp, mass_temp
-  ! 
+  LOGICAL, INTENT(IN)                                 :: ok_alw
+  INTEGER, INTENT(IN)                                 :: flag_aerosol
+  REAL, DIMENSION(klon,klev),     INTENT(in)          :: pdel, zrho
+  REAL, DIMENSION(klon,klev,naero_tot),   INTENT(in)  :: m_allaer, m_allaer_pi
   !
   INTEGER inu, i, k
-  INTEGER :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
-  INTEGER :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
-
+  REAL :: zdh(klon,klev)
+  REAL :: tmp_var, tmp_var_pi
   !
   !--absorption coefficient for CIDUST
   REAL:: alpha_abs_CIDUST_16bands(nbands_lw_rrtm)   !--unit m2/g 
   DATA alpha_abs_CIDUST_16bands /              & 
-       0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & 
-       0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  /
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
   !
   IF (NLW.NE.nbands_lw_rrtm) THEN
-     print *,'Erreur NLW doit etre egal a 16 pour cette routine'
-     stop
+     PRINT *,'Erreur NLW doit etre egal a 16 pour cette routine'
+     STOP
   ENDIF
   ! 
-  IF (aerosol_couple) THEN 
-  ! 
-      do i=1,nbtr
-      select case(trim(solsym(i)))
-         case ("ASBCM")
-            id_ASBCM = i
-         case ("ASPOMM")
-            id_ASPOMM = i
-         case ("ASSO4M")
-            id_ASSO4M = i
-         case ("ASMSAM")
-            id_ASMSAM = i
-         case ("CSSO4M")
-            id_CSSO4M = i
-         case ("CSMSAM")
-            id_CSMSAM = i
-         case ("SSSSM")
-            id_SSSSM = i
-         case ("CSSSM")
-            id_CSSSM = i
-         case ("ASSSM")
-            id_ASSSM = i
-         case ("CIDUSTM")
-            id_CIDUSTM = i
-         case ("AIBCM")
-            id_AIBCM = i
-         case ("AIPOMM")
-            id_AIPOMM = i
-         case ("ASNO3M")
-            id_ASNO3M = i
-         case ("CSNO3M")
-            id_CSNO3M = i
-         case ("CINO3M")
-            id_CINO3M = i
-         end select
-    enddo
+  IF (ok_alw) THEN                                   !--aerosol LW effects
+   !
+   IF (flag_aerosol.EQ.5.OR.flag_aerosol.EQ.6) THEN  !--Dust
+    !
+    zdh(:,:)=pdel(:,:)/(RG*zrho(:,:))      !-- m
     !
     DO k=1, klev
-       zdp(:,k) = (paprs(:,k)-paprs(:,k+1))/RG !--kg/m2
-    ENDDO
-    !
-    !--for now only consider dust
-    mass_temp(:,:)=tr_seri(:,:,id_CIDUSTM)  !--kg/kg
-    !
-    DO inu=1,NLW
-      !
-      !--total aerosol
-      tau_aero_lw_rrtm(:,:,2,inu) = mass_temp(:,:)*zdp(:,:)*1000.*alpha_abs_CIDUST_16bands(inu)
-      !--no aerosol at all
-      tau_aero_lw_rrtm(:,:,1,inu) = 0.0
+      DO i=1, klon
+         !
+         tmp_var   =m_allaer(i,k,id_CIDUSTM_phy)   /1.e6*zdh(i,k) !-- g/m2
+         tmp_var_pi=m_allaer_pi(i,k,id_CIDUSTM_phy)/1.e6*zdh(i,k) !-- g/m2
+         !
+         DO inu=1, NLW
+           !
+           !--total aerosol
+           tau_aero_lw_rrtm(i,k,2,inu) = MAX(1.e-15,tmp_var*alpha_abs_CIDUST_16bands(inu))
+           !--natural aerosol 
+           tau_aero_lw_rrtm(i,k,1,inu) = MAX(1.e-15,tmp_var_pi*alpha_abs_CIDUST_16bands(inu))
+           !
+         ENDDO
+      ENDDO
       !
     ENDDO
     !
-    !--avoid very small values
-    tau_aero_lw_rrtm = MAX(tau_aero_lw_rrtm,1.e-15)
-    ! 
-  ELSE !--not aerosol_couple
+   ENDIF
+   !
+  ELSE !--no aerosol LW effects
     !
-    !--no LW effects if not coupled to INCA
     tau_aero_lw_rrtm = 1.e-15
   ENDIF
Index: LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90	(revision 2737)
+++ LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90	(revision 2738)
@@ -1,5 +1,5 @@
 ! $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,4 +30,5 @@
   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,6 +282,6 @@
   ! 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,4 +311,7 @@
        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
