Index: /trunk/LMDZ.MARS/libf/phymars/improvedclouds.F
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/improvedclouds.F	(revision 530)
+++ /trunk/LMDZ.MARS/libf/phymars/improvedclouds.F	(revision 531)
@@ -170,10 +170,10 @@
 
         do i=1,nbin_cld-1
-          rad_cld(i+1)  = rad_cld(i) * CBRT(vrat_cld) !!**(1./3.)
+          rad_cld(i+1)  = rad_cld(i) * vrat_cld**(1./3.)
           vol_cld(i+1)  = vol_cld(i) * vrat_cld
         enddo
         
         do i=1,nbin_cld
-          rb_cld(i+1)= CBRT( (2.*vrat_cld) / (vrat_cld+1.) ) *
+          rb_cld(i+1)= ( (2.*vrat_cld) / (vrat_cld+1.) )**(1./3.) *
      &      rad_cld(i)
           dr_cld(i)  = rb_cld(i+1) - rb_cld(i)
@@ -364,5 +364,5 @@
           rhocloud(ig,l) = min(max(rhocloud(ig,l),rho_ice),rho_dust)
           rice(ig,l) =
-     &      CBRT( Mo / No * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
+     &      CBRT( real(Mo)/real(No) * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
 c          nuice(ig,l)=nuice_ref ! used for rad. transfer calculations
           if ((Mo.lt.1.e-20) .or. (No.le.1)) rice(ig,l) = 1.e-8
@@ -484,5 +484,5 @@
           rice_out(ig,l)=rice(ig,l)
           rice(ig,l) =
-     &      CBRT( Mo / No * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
+     &      CBRT( real(Mo)/real(No) * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
           if ((Mo.lt.1.e-20) .or. (No.le.1)) rice(ig,l) = 1.e-8
           rice_out(ig,l)=rice(ig,l)-rice_out(ig,l)