Index: /issm/trunk/src/m/classes/public/display/displaybc.m
===================================================================
--- /issm/trunk/src/m/classes/public/display/displaybc.m	(revision 2207)
+++ /issm/trunk/src/m/classes/public/display/displaybc.m	(revision 2208)
@@ -24,5 +24,4 @@
 disp(sprintf('\n      prognostic:'));
 fielddisplay(md,'spcthickness','constraints flag list (first column) and values (second column)');
-fielddisplay(md,'neumannvalues_prog','values of the Neumann [m/a]');
 
 disp(sprintf('\n      thermal:'));
Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 2207)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 2208)
@@ -204,4 +204,12 @@
 %DIAGNOSTIC
 if md.analysis_type==DiagnosticAnalysisEnum,
+
+	%CHECK THAT WE ARE NOT FULLY CONSTRAINED
+	if strcmpi(md.type,'2d'),
+		if isempty(find(~md.spcvelocity(:,1:2))),
+			disp(['model ' md.name ' is totally constrained, no need to solve!']);
+			bool=0;return;
+		end
+	end
 
 	%HUTTER ON ICESHELF WARNING