Index: /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90
===================================================================
--- /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3686)
+++ /trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90	(revision 3687)
@@ -450,4 +450,7 @@
 !--------------------------------------------------
 !     Effective radius and variance of the aerosols
+!     Madeleine's PhD (eq. 2.3-2.4):
+!     --> r_eff = <r^3> / <r^2>
+!     --> nu_eff = <r^4>*<r^2> / <r^3>^2 - 1
 !--------------------------------------------------
       ! Radiative Hazes
@@ -456,5 +459,5 @@
             ! Spherical aerosols
             sig = 0.2
-            where (mp2m_rc_sph(:,:) > 1e-10)
+            where (mp2m_rc_sph(:,:) > 1e-9)
                reffrad(:,:,1) = mp2m_rc_sph(:,:) * exp(5.*sig**2 / 2.)
             elsewhere
@@ -468,5 +471,5 @@
             ! Fractal aerosols
             sig = 0.35
-            where (mp2m_rc_fra(:,:) > 1e-10)
+            where (mp2m_rc_fra(:,:) > 1e-8)
                reffrad(:,:,2) = mp2m_rc_fra(:,:) * exp(5.*sig**2 / 2.)
                elsewhere
@@ -543,5 +546,5 @@
 !     Net atmospheric radiative cooling rate from C2H2 (K.s-1):
 !     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-    !   dtlw_hcn_c2h2=0.
+      dtlw_hcn_c2h2=0.
       if (cooling) then
          CALL cooling_hcn_c2h2(ngrid,nlayer,pplay,&
@@ -1103,5 +1106,7 @@
         end do
         ! adding c2h2 if cooling active
-        dtlw(ig,l)=dtlw(ig,l)+dtlw_hcn_c2h2(ig,l)
+        if (cooling) then
+         dtlw(ig,l)=dtlw(ig,l)+dtlw_hcn_c2h2(ig,l)
+        endif
 
         !SW