Index: /issm/trunk-jpl/src/m/consistency/checkfield.m
===================================================================
--- /issm/trunk-jpl/src/m/consistency/checkfield.m	(revision 25635)
+++ /issm/trunk-jpl/src/m/consistency/checkfield.m	(revision 25636)
@@ -139,5 +139,4 @@
 end
 
-
 %check cell
 if getfieldvalue(options,'cell',0);
@@ -187,5 +186,5 @@
 	lowerbound=getfieldvalue(options,'>=');
 	field2=reshape(field,prod(size(field)),1);
-	if getfieldvalue(options,'timeseries',0), field2=reshape(field(1:end-1,:),prod(size(field(1:end-1,:))),1); end
+	if getfieldvalue(options,'timeseries',0),       field2=reshape(field(1:end-1,:),prod(size(field(1:end-1,:))),1); end
 	if getfieldvalue(options,'singletimeseries',0), field2=reshape(field(1,:),prod(size(field(1,:))),1); end
 	if any(field2<lowerbound),
@@ -197,5 +196,5 @@
 	lowerbound=getfieldvalue(options,'>');
 	field2=reshape(field,prod(size(field)),1);
-	if getfieldvalue(options,'timeseries',0), field2=reshape(field(1:end-1,:),prod(size(field(1:end-1,:))),1); end
+	if getfieldvalue(options,'timeseries',0),       field2=reshape(field(1:end-1,:),prod(size(field(1:end-1,:))),1); end
 	if getfieldvalue(options,'singletimeseries',0), field2=reshape(field(1,:),prod(size(field(1,:))),1); end
 	if any(field2<=lowerbound),
Index: /issm/trunk-jpl/src/m/exp/contourlevelzero.m
===================================================================
--- /issm/trunk-jpl/src/m/exp/contourlevelzero.m	(revision 25635)
+++ /issm/trunk-jpl/src/m/exp/contourlevelzero.m	(revision 25636)
@@ -9,7 +9,10 @@
 %   See also: PLOT_CONTOUR
 
+if isempty(mask)
+	error('mask provided is empty');
+end
+
 %process data 
 if dimension(md.mesh)==3,
-% 	error('contourlevelzero error message: routine not supported for 3d meshes, project on a layer');
 	x = md.mesh.x2d;
 	y = md.mesh.y2d;
@@ -27,7 +30,6 @@
 end
 
-if isempty(mask), error('mask provided is empty'); end
 if dimension(md.mesh)==3,
-	if length(mask)~=md.mesh.numberofvertices2d, error('mask provided should be specified at the vertices of the mesh'); end
+	if length(mask)~=md.mesh.numberofvertices2d, error('Since this is a 3d model, you need to project the mask on a layer (e.g. project2d(md,mask,0))'); end
 else
 	if length(mask)~=md.mesh.numberofvertices, error('mask provided should be specified at the vertices of the mesh'); end