Index: trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F	(revision 2984)
+++ trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F	(revision 2988)
@@ -298,8 +298,10 @@
         spenttime = 0.
         ending_ts=.false.
-        print*, "ig,l :", ig, l
         DO while (.not.ending_ts)
           microtimestep=ptimestep/real(zimicro(ig,l))
-          zq0(1:ngrid,1:nlay,1:nq) = zq
+c         Initialize tracers for scavenging + hdo computations (JN)
+          DO i=1,nq
+             zq0(ig,l,i)=zq(ig,l,i)
+          ENDDO !i=1,nq
 
           ! Check if we are integrating over ptimestep
@@ -577,7 +579,8 @@
      &          subpdtcloud(ig,l))*microtimestep
 
-c         Prevent negative tracers
-          WHERE( zq(1:ngrid,1:nlay,1:nq) < 1.e-30 )
-     &       zq(1:ngrid,1:nlay,1:nq) = 1.e-30
+c         Prevent negative tracers ! JN
+          DO i=1,nq
+            IF(zq(ig,l,i).lt.1.e-30) zq(ig,l,i)=1.e-30
+          ENDDO !i=1,nq
 
 c         Increment time spent in here
@@ -758,9 +761,9 @@
 c       zimicro = 30
 c      Coefficients good enough for present-day Mars :
-       alpha = 1.87485684e+09
-       beta = 1.45655856e+00
+c       alpha = 1.87485684e+09
+c       beta = 1.45655856e+00
 c      Coefficients covering high obliquity scenarios :
-c       alpha=3.36711332e+15
-c       beta=1.98636414e+00
+       alpha=3.36711332e+15
+       beta=1.98636414e+00
 c       nimicro=min(max(alpha*abs(potcond)**beta,5.),7000.)
 c       zimicro=ceiling((ptimestep/timeleft)*nimicro)