Index: trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F	(revision 2322)
+++ trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F	(revision 2323)
@@ -338,4 +338,6 @@
 c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       zqi(1:ngrid,1:nlay,1:nq)=pq(1:ngrid,1:nlay,1:nq)
+     &                         +pdqfi(1:ngrid,1:nlay,1:nq)*ptimestep
+      zq0(1:ngrid,1:nlay,1:nq)=pq(1:ngrid,1:nlay,1:nq) !MVals: keep the input value
      &                         +pdqfi(1:ngrid,1:nlay,1:nq)*ptimestep
       zt(1:ngrid,1:nlay)=pt(1:ngrid,1:nlay)