Index: /trunk/LMDZ.GENERIC/libf/phystd/interpolateH2Ocont.F90
===================================================================
--- /trunk/LMDZ.GENERIC/libf/phystd/interpolateH2Ocont.F90	(revision 536)
+++ /trunk/LMDZ.GENERIC/libf/phystd/interpolateH2Ocont.F90	(revision 537)
@@ -51,5 +51,7 @@
       amagatF=(273.15/temp)*(presF/101325.0)
 
-      if(firstcall)then
+      if(firstcall)then ! called by sugas_corrk only
+         print*,'----------------------------------------------------'
+         print*,'Initialising H2O continuum from MT_CKD data...'
 
 !     1.1 Open the ASCII files
Index: /trunk/LMDZ.GENERIC/libf/phystd/radinc_h.F90
===================================================================
--- /trunk/LMDZ.GENERIC/libf/phystd/radinc_h.F90	(revision 536)
+++ /trunk/LMDZ.GENERIC/libf/phystd/radinc_h.F90	(revision 537)
@@ -8,5 +8,5 @@
 !======================================================================
 !
-!     RADINC.H    RADiation INCludes
+!     RADINC.H 
 !
 !     Includes for the radiation code; RADIATION LAYERS, LEVELS,
@@ -54,5 +54,5 @@
 !     NSIZEMAX  The maximum number of aerosol particle sizes
 !
-!----------------------------------------------------------------------C
+!----------------------------------------------------------------------
 
       integer, parameter :: L_NLAYRAD  = llm
@@ -76,19 +76,9 @@
       integer, parameter :: NTstar = 500
       integer, parameter :: NTstop = 9000 ! new default for all non hot Jupiter runs
-      !integer, parameter :: NTstop = 6000 ! for GJ581d / earlymars runs
 
-! Maximum number of grain size classes
-! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-!   This parameter has to be set to the maximum number of particle
-!   sizes contained in the optical parameter database; For example,
-!   if only one grain size is used to describe dust, and 30 are used
-!   to describe water-ice crystals in the visible and 15 in the IR,
-!   nsizemax has to be set to 30.
-!   If only one grain size is considered for all the aerosols, set
-!   this parameter to 1 and convolution will be turned off during
-!   the radiative calculations.
-
+      ! Maximum number of grain size classes for aerosol convolution:
+      ! This must correspond to size of largest dataset used for aerosol 
+      ! optical properties in datagcm folder.
       integer, parameter :: nsizemax = 60
-!      integer, parameter :: nsizemax = 1
 
       character (len=20) :: corrkdir