Index: trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90	(revision 732)
+++ trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90	(revision 736)
@@ -580,5 +580,5 @@
                       !recompute q1 with new water flux from surface  
                       zq(ig,1,iq) = (zmass(ig,1)*(pq(ig,1,iq)+ptimestep*pdqfi(ig,1,iq))  &
-		                            +zfluxq(ig,2)*zcq(ig,1)-dqsdif_total(ig))     &
+		                            +zfluxq(ig,2)*zcq(ig,2)-dqsdif_total(ig))     &
                                  / (zmass(ig,1)+(1.-zdq(ig,2))*zfluxq(ig,2))		     
                   end if