Index: /issm/trunk/src/m/solutions/cielo/control_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/control_core.m	(revision 1224)
+++ /issm/trunk/src/m/solutions/cielo/control_core.m	(revision 1225)
@@ -27,5 +27,5 @@
 %Take car of Stokes : compute slope and get spc once for all
 if isstokes,
-	[inputs m_ds]=ControlProcessStokes(models,inputs);
+	[inputs m_ds]=ControlPrepareStokes(models,inputs);
 end
 
Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m	(revision 1224)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core_nonlinear.m	(revision 1225)
@@ -49,4 +49,7 @@
 		[soln(count).u_f]=Solver(K_ff,p_f,[],m.parameters);
 
+		%Get residue
+		residu=norm(K_ff*soln(count).u_f-p_f,2);
+
 		%Merge back to g set
 		[soln(count).u_g]= Mergesolutionfromftog( soln(count).u_f, m.Gmn, m.ys, m.nodesets ); 
@@ -67,4 +70,7 @@
 		ndu=norm(dug,2); 
 		nu=norm(soln(count-1).u_g,2); 
+
+		%Residue criterion
+		displaystring(m.parameters.debug,'%s%g','      convergence criterion: norm(residue)=',residu);
 
 		%Relative criterion