Index: /issm/trunk/src/m/solutions/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/diagnostic_core.m	(revision 5926)
+++ /issm/trunk/src/m/solutions/diagnostic_core.m	(revision 5927)
@@ -61,5 +61,5 @@
 	end
 
-	if dim==3 & (ismacayealpattyn | ishutter),
+	if dim==3 & (ismacayealpattyn | ishutter) & ~isstokes,
 	
 		displaystring(verbose,'\n%s',['computing vertical velocities...']);
Index: /issm/trunk/src/m/solvers/solver_stokescoupling_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solvers/solver_stokescoupling_nonlinear.m	(revision 5926)
+++ /issm/trunk/src/m/solvers/solver_stokescoupling_nonlinear.m	(revision 5927)
@@ -44,18 +44,19 @@
 		%Figure out if convergence have been reached
 		converged=convergence(K_ff,p_f,uf_horiz,old_uf_horiz,femmodel.parameters);
+		converged=1;
 
-		%Then compute vertical velocity
-		femmodel=SetCurrentConfiguration(femmodel,DiagnosticVertAnalysisEnum);
-		[K_gg,K_ff,K_fs,p_g,p_f,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
-
-		[K_ff, K_fs] = Reducematrixfromgtof( K_gg,  femmodel.nodesets,femmodel.parameters); 
-		p_f = Reduceloadfromgtof( p_g,  K_fs, femmodel.ys, femmodel.nodesets,femmodel.parameters);
-
-		displaystring(femmodel.parameters.Verbose>1,'%s%g','      condition number of stiffness matrix: ',condest(K_ff));
-		uf_vert=Solver(K_ff,p_f,[],femmodel.parameters);
-
-		ug_vert= Mergesolutionfromftog( uf_vert, femmodel.ys, femmodel.nodesets,femmodel.parameters); 
-
-		[femmodel.elements,femmodel.materials]=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ug_vert);
+%		%Then compute vertical velocity
+%		femmodel=SetCurrentConfiguration(femmodel,DiagnosticVertAnalysisEnum);
+%		[K_gg,K_ff,K_fs,p_g,p_f,kmax]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+%
+%		[K_ff, K_fs] = Reducematrixfromgtof( K_gg,  femmodel.nodesets,femmodel.parameters); 
+%		p_f = Reduceloadfromgtof( p_g,  K_fs, femmodel.ys, femmodel.nodesets,femmodel.parameters);
+%
+%		displaystring(femmodel.parameters.Verbose>1,'%s%g','      condition number of stiffness matrix: ',condest(K_ff));
+%		uf_vert=Solver(K_ff,p_f,[],femmodel.parameters);
+%
+%		ug_vert= Mergesolutionfromftog( uf_vert, femmodel.ys, femmodel.nodesets,femmodel.parameters); 
+%
+%		[femmodel.elements,femmodel.materials]=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ug_vert);
 
 	end