Index: ../trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp
===================================================================
--- ../trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18366)
+++ ../trunk-jpl/src/c/solutionsequences/solutionsequence_fct.cpp	(revision 18367)
@@ -247,10 +247,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*(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.);
 		d = 1.;
