Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18360)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18361)
@@ -232,8 +232,8 @@
 
 		/*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.);
@@ -389,4 +389,7 @@
 	VecToMPISerial(&ml_serial  ,Ml_petsc,IssmComm::GetComm());
 
+	/*Push uL to elements*/
+	InputUpdateFromSolutionx(femmodel,u_serial);
+
 	/*Anti diffusive fluxes*/
 	Vec         Fbar            = NULL;