Index: ../trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp
===================================================================
--- ../trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18360)
+++ ../trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18361)
@@ -231,10 +231,10 @@
 		}
 
 		/*Compute Qis and Ris*/
-		double Qi_plus  = ml_serial[row]/deltat*(3. - u[row]);
-		double Qi_minus = ml_serial[row]/deltat*(2. - u[row]);
-		//double Qi_plus  = ml_serial[row]/deltat*(ulmax[row] - u[row]);
-		//double Qi_minus = ml_serial[row]/deltat*(ulmin[row] - u[row]);
+		//double Qi_plus  = ml_serial[row]/deltat*(3. - u[row]);
+		//double Qi_minus = ml_serial[row]/deltat*(2. - u[row]);
+		double Qi_plus  = ml_serial[row]/deltat*(ulmax[row] - u[row]);
+		double Qi_minus = ml_serial[row]/deltat*(ulmin[row] - u[row]);
 		_assert_(Qi_plus  >= 0.);
 		_assert_(Qi_minus <= 0.);
 		d = 1.;
@@ -388,6 +388,9 @@
 	VecToMPISerial(&u_serial   ,u       ,IssmComm::GetComm());
 	VecToMPISerial(&ml_serial  ,Ml_petsc,IssmComm::GetComm());
 
+	/*Push uL to elements*/
+	InputUpdateFromSolutionx(femmodel,u_serial);
+
 	/*Anti diffusive fluxes*/
 	Vec         Fbar            = NULL;
 	IssmDouble *Ri_minus_serial = NULL;
