Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 2722)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 2723)
@@ -327,4 +327,18 @@
 	if any(md.spcthickness(find(md.gridonboundary))~=1),
 		error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcthickness']);
+	end
+end
+
+%BALANCEDVELOCITIES
+if md.analysis_type==BalancedvelocitiesAnalysisEnum
+
+	%VELOCITIES MELTING AND ACCUMULATION
+	fields={'vx','vy','accumulation','melting'};
+	checksize(md,fields,[md.numberofgrids 1]);
+	checknan(md,fields);
+
+	%SPC
+	if any(md.spcvelocity(find(md.gridonboundary),[1:2])~=1),
+		error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcvelocity']);
 	end
 end