Index: /trunk/LMDZ.GENERIC/README
===================================================================
--- /trunk/LMDZ.GENERIC/README	(revision 746)
+++ /trunk/LMDZ.GENERIC/README	(revision 747)
@@ -746,2 +746,6 @@
    => This solves a water conservation problem arising when the code tries to
    evaporate over dry land.
+
+== 26/07/2012 == AS
+- Bug fix in physiq : the size of OLR_nu is L_NSPECTI and not L_NSPECTV
+- A more robust aerosol_mod + iniaerosol : problems with ifort+parallel solved, while still OK with other compilers and seq runs.
Index: /trunk/LMDZ.GENERIC/libf/phystd/aerosol_mod.F90
===================================================================
--- /trunk/LMDZ.GENERIC/libf/phystd/aerosol_mod.F90	(revision 746)
+++ /trunk/LMDZ.GENERIC/libf/phystd/aerosol_mod.F90	(revision 747)
@@ -1,13 +1,16 @@
 !==================================================================
 module aerosol_mod
+implicit none
+save
 !==================================================================
 
-!  aerosol indexes: these are initialized in iniaerosol.F and should be 0 if the
+!  aerosol indexes: these are initialized to be 0 if the
 !                 corresponding aerosol was not activated in callphys.def
-      integer, save :: iaero_co2 
-      integer, save :: iaero_h2o
-      integer, save :: iaero_dust
-      integer, save :: iaero_h2so4
-      logical, save :: noaero
+!                 -- otherwise a value is given in iniaerosol
+      integer :: iaero_co2 = 0 
+      integer :: iaero_h2o = 0
+      integer :: iaero_dust = 0
+      integer :: iaero_h2so4 = 0
+      logical :: noaero = .false.
       
 !==================================================================
Index: /trunk/LMDZ.GENERIC/libf/phystd/iniaerosol.F
===================================================================
--- /trunk/LMDZ.GENERIC/libf/phystd/iniaerosol.F	(revision 746)
+++ /trunk/LMDZ.GENERIC/libf/phystd/iniaerosol.F	(revision 747)
@@ -21,39 +21,30 @@
       integer ia
 
-      ! Initialize aerosol indexes to zero:
-      ! NB: iaero_* indexes are in module aerosol_mod
-
-      iaero_co2   =0
-      iaero_h2o   =0
-      iaero_dust  =0
-      iaero_h2so4 =0
-      noaero = .false.
- 
       ia=0
       if (aeroco2) then
          ia=ia+1
          iaero_co2=ia
-         write(*,*) 'CO2 aerosol = ', iaero_co2
       endif
+      write(*,*) '--- CO2 aerosol = ', iaero_co2
  
       if (aeroh2o) then
          ia=ia+1
          iaero_h2o=ia
-         write(*,*) 'H2O aerosol = ', iaero_h2o
       endif
+      write(*,*) '--- H2O aerosol = ', iaero_h2o
 
       if (dusttau.gt.0) then
          ia=ia+1
          iaero_dust=ia
-         write(*,*) 'Dust aerosol = ', iaero_dust
       endif
+      write(*,*) '--- Dust aerosol = ', iaero_dust
 
       if (aeroh2so4) then
          ia=ia+1
          iaero_h2so4=ia
-         write(*,*) 'H2SO4 aerosol = ', iaero_h2so4
       endif
+      write(*,*) '--- H2SO4 aerosol = ', iaero_h2so4
       
-      write(*,*) 'Number of aerosols= ', ia
+      write(*,*) '=== Number of aerosols= ', ia
       
 ! For the zero aerosol case, we currently make a dummy co2 aerosol which is zero everywhere.