Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18366)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18367)
@@ -248,8 +248,8 @@
 		//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*(uLmax_global - u[row]);
-		double Qi_minus = ml_serial[row]/deltat*(uLmin_global - 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*(uLmax_global - u[row]);
+		//double Qi_minus = ml_serial[row]/deltat*(uLmin_global - u[row]);
 		_assert_(Qi_plus  >= 0.);
 		_assert_(Qi_minus <= 0.);