Index: trunk/LMDZ.MARS/libf/phymars/improvedCO2clouds_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/improvedCO2clouds_mod.F	(revision 2151)
+++ trunk/LMDZ.MARS/libf/phymars/improvedCO2clouds_mod.F	(revision 2152)
@@ -198,6 +198,6 @@
 c       dr_cld is the width of each rad_cldco2 bin.
 
-        vrat_cld = dlog(rmax_cld/rmin_cld) / float(nbinco2_cld-1) *3.
-        vrat_cld = dexp(vrat_cld)
+        vrat_cld = log(rmax_cld/rmin_cld) / float(nbinco2_cld-1) *3.
+        vrat_cld = exp(vrat_cld)
         rb_cldco2(1)  = rbmin_cld
         rad_cldco2(1) = rmin_cld
@@ -226,5 +226,5 @@
         print*,'-----------------------------------'
         do i=1,nbinco2_cld+1
-            rb_cldco2(i) = dlog(rb_cldco2(i))  !! we save that so that it is not computed
+            rb_cldco2(i) = log(rb_cldco2(i))  !! we save that so that it is not computed
                                          !! at each timestep and gridpoint
         enddo
@@ -402,5 +402,5 @@
 
               Mo =4.*pi*rho_dust*No*rdust(ig,l)**(3.)
-     &           *dexp(9.*nuiceco2_ref/2.)/3. ! in Madeleine et al 2011
+     &           *exp(9.*nuiceco2_ref/2.)/3. ! in Madeleine et al 2011
      
               No = zq(ig,l,igcm_dust_number)* tauscaling(ig)+1.e-30
@@ -410,5 +410,5 @@
 
               Rn = rdust(ig,l)
-              Rn = -dlog(Rn) 
+              Rn = -log(Rn) 
               Rm = Rn - 3. * sigma_iceco2*sigma_iceco2  
               n_derf = derf( (rb_cldco2(1)+Rn) *dev2)
@@ -448,5 +448,5 @@
                 No = zq(ig,l,igcm_ccn_number)* tauscaling(ig) + 1.e-30
                 Rn = rice(ig,l)
-                Rn = -dlog(Rn) 
+                Rn = -log(Rn) 
                 Rm = Rn - 3. * sigma_ice*sigma_ice  
                 n_derf = derf( (rb_cldco2(1)+Rn) *dev3)
@@ -473,5 +473,5 @@
               dMh2o = 0.
               do i = 1, nbinco2_cld
-                 Proba    =1.0-dexp(-1.*microtimestep*rate(i))
+                 Proba    =1.0-exp(-1.*microtimestep*rate(i))
                  dN       = dN + n_aer(i) * Proba
                  dM       = dM + m_aer(i) * Proba             
@@ -479,5 +479,5 @@
               if (co2useh2o) then
                  do i = 1, nbinco2_cld
-                    Probah2o = 1.0-dexp(-1.*microtimestep*rateh2o(i))
+                    Probah2o = 1.0-exp(-1.*microtimestep*rateh2o(i))
                     dNh2o    = dNh2o + n_aer_h2oice(i) * Probah2o
                     dMh2o    = dMh2o + m_aer_h2oice(i) * Probah2o
@@ -535,6 +535,7 @@
      &         * tauscaling(ig)+1.e-30.ge. 1)THEN   ! we trigger crystal growth
               
-              call updaterice_microco2(zq(ig,l,igcm_co2_ice),
-     &            zq(ig,l,igcm_ccnco2_mass),zq(ig,l,igcm_ccnco2_number),
+              call updaterice_microco2(dble(zq(ig,l,igcm_co2_ice)),
+     &            dble(zq(ig,l,igcm_ccnco2_mass)),
+     &            dble(zq(ig,l,igcm_ccnco2_number)),
      &            tauscaling(ig),riceco2(ig,l),rhocloudco2(ig,l))