Index: /trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90
===================================================================
--- /trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90	(revision 3377)
+++ /trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90	(revision 3378)
@@ -629,5 +629,5 @@
     ! use Pluto.old routines
      if (is_master) write(*,*) trim(rname)//": use vdifc from old Pluto ?"
-     oldplutovdifc=.true. ! default value
+     oldplutovdifc=.false. ! default value
      call getin_p("oldplutovdifc",oldplutovdifc)
      if (is_master) write(*,*) trim(rname)//": oldplutovdifc = ",oldplutovdifc
@@ -635,5 +635,5 @@
      if (is_master) write(*,*) trim(rname)//&
        ": call pluto.old correlated-k radiative transfer ?"
-     oldplutocorrk=.true. ! default value
+     oldplutocorrk=.false. ! default value
      call getin_p("oldplutocorrk",oldplutocorrk)
      if (is_master) write(*,*) trim(rname)//": oldplutocorrk = ",oldplutocorrk
@@ -641,5 +641,5 @@
      if (is_master) write(*,*) trim(rname)//&
        ": call pluto.old sedimentation ?"
-     oldplutosedim=.true. ! default value
+     oldplutosedim=.false. ! default value
      call getin_p("oldplutosedim",oldplutosedim)
      if (is_master) write(*,*) trim(rname)//": oldplutosedim = ",oldplutosedim
Index: /trunk/LMDZ.PLUTO/libf/phypluto/radinc_h.F90
===================================================================
--- /trunk/LMDZ.PLUTO/libf/phypluto/radinc_h.F90	(revision 3377)
+++ /trunk/LMDZ.PLUTO/libf/phypluto/radinc_h.F90	(revision 3378)
@@ -11,9 +11,9 @@
 !======================================================================
 !
-!     RADINC.H 
+!     RADINC.H
 !
 !     Includes for the radiation code; RADIATION LAYERS, LEVELS,
 !     number of spectral intervals. . .
-! 
+!
 !======================================================================
 
@@ -27,10 +27,10 @@
 
 !     L_NLEVRAD corresponds to the surface - i.e., the GCM Level that
-!     is at the surface.  PLEV(L_NLEVRAD) = P(J,I)+PTROP, 
+!     is at the surface.  PLEV(L_NLEVRAD) = P(J,I)+PTROP,
 !     PLEV(2) = PTROP, PLEV(1) = ptop
 
 !     L_NLAYRAD is the number of radiation code layers
 !     L_NLEVRAD is the number of radiation code levels.  Level N is the
-!               top of layer N. 
+!               top of layer N.
 !
 !     L_NSPECTI is the number of IR spectral intervals
@@ -39,5 +39,5 @@
 !               GAUSS POINT 17 (aka the last one) is the special case
 !
-!     L_NPREF   is the number of reference pressures that the 
+!     L_NPREF   is the number of reference pressures that the
 !               k-coefficients are calculated on
 !     L_PINT    is the number of Lagrange interpolated reference
@@ -74,5 +74,5 @@
       real,    parameter :: L_TAUMAX  = 35
 
-      ! For Planck function integration: 
+      ! For Planck function integration:
       ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       ! Integration boundary temperatures are NTstart/NTfac and Ntstop/NTfac
@@ -86,7 +86,7 @@
 
       ! Maximum number of grain size classes for aerosol convolution:
-      ! This must correspond to size of largest dataset used for aerosol 
+      ! This must correspond to size of largest dataset used for aerosol
       ! optical properties in datagcm folder.
-      integer, parameter :: nsizemax = 60
+      integer, parameter :: nsizemax = 200
 
       character(len=100),save :: corrkdir
@@ -105,5 +105,5 @@
   real*8, intent(in) :: tplanckmax
   real*8, intent(in) :: dtplanck
-  
+
   L_NLAYRAD = nbp_lev
   L_LEVELS  = 2*(nbp_lev-1)+3
@@ -113,5 +113,5 @@
   NTstart = int(tplanckmin * NTfac)
   NTstop  = int(tplanckmax * NTfac)
- 
+
   end subroutine