Index: /issm/trunk-jpl/src/m/classes/mesh3dprisms.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/mesh3dprisms.m	(revision 18076)
+++ /issm/trunk-jpl/src/m/classes/mesh3dprisms.m	(revision 18077)
@@ -18,6 +18,4 @@
 		hemisphere                  = NaN
 
-		elementonbase                = NaN
-		elementonsurface            = NaN
 		vertexonbase                 = NaN
 		vertexonsurface             = NaN
@@ -163,7 +161,5 @@
 			fielddisplay(obj,'numberoflayers','number of extrusion layers');
 			fielddisplay(obj,'vertexonbase','lower vertices flags list');
-			fielddisplay(obj,'elementonbase','lower elements flags list');
 			fielddisplay(obj,'vertexonsurface','upper vertices flags list');
-			fielddisplay(obj,'elementonsurface','upper elements flags list');
 			fielddisplay(obj,'uppervertex','upper vertex list (NaN for vertex on the upper surface)');
 			fielddisplay(obj,'upperelements','upper element list (NaN for element on the upper layer)');
Index: /issm/trunk-jpl/src/m/classes/mesh3dprisms.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/mesh3dprisms.py	(revision 18076)
+++ /issm/trunk-jpl/src/m/classes/mesh3dprisms.py	(revision 18077)
@@ -27,6 +27,4 @@
 		self.hemisphere                  = float('NaN');
 
-		self.elementonbase                = float('NaN');
-		self.elementonsurface            = float('NaN');
 		self.vertexonbase                 = float('NaN');
 		self.vertexonsurface             = float('NaN');
@@ -75,7 +73,5 @@
 		string="%s\n%s"%(string,fielddisplay(self,"numberoflayers","number of extrusion layers"))
 		string="%s\n%s"%(string,fielddisplay(self,"vertexonbase","lower vertices flags list"))
-		string="%s\n%s"%(string,fielddisplay(self,"elementonbase","lower elements flags list"))
 		string="%s\n%s"%(string,fielddisplay(self,"vertexonsurface","upper vertices flags list"))
-		string="%s\n%s"%(string,fielddisplay(self,"elementonsurface","upper elements flags list"))
 		string="%s\n%s"%(string,fielddisplay(self,"uppervertex","upper vertex list (-1 for vertex on the upper surface)"))
 		string="%s\n%s"%(string,fielddisplay(self,"upperelements","upper element list (-1 for element on the upper layer)"))
Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 18076)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 18077)
@@ -711,6 +711,4 @@
 
 			%bedinfo and surface info
-			md.mesh.elementonbase=project3d(md,'vector',ones(md.mesh.numberofelements2d,1),'type','element','layer',1);
-			md.mesh.elementonsurface=project3d(md,'vector',ones(md.mesh.numberofelements2d,1),'type','element','layer',md.mesh.numberoflayers-1);
 			md.mesh.vertexonbase=project3d(md,'vector',ones(md.mesh.numberofvertices2d,1),'type','node','layer',1);
 			md.mesh.vertexonsurface=project3d(md,'vector',ones(md.mesh.numberofvertices2d,1),'type','node','layer',md.mesh.numberoflayers);
@@ -943,6 +941,4 @@
 			if isfield(structmd,'upperelements'), md.mesh.upperelements=structmd.upperelements; end
 			if isfield(structmd,'lowerelements'), md.mesh.lowerelements=structmd.lowerelements; end
-			if isfield(structmd,'elementonbase'), md.mesh.elementonbase=structmd.elementonbase; end
-			if isfield(structmd,'elementonsurface'), md.mesh.elementonsurface=structmd.elementonsurface; end
 			if isfield(structmd,'nodeonsurface'), md.mesh.vertexonsurface=structmd.nodeonsurface; end
 			if isfield(structmd,'nodeonbase'), md.mesh.vertexonbase=structmd.nodeonbase; end
Index: /issm/trunk-jpl/src/m/classes/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.py	(revision 18076)
+++ /issm/trunk-jpl/src/m/classes/model.py	(revision 18077)
@@ -618,6 +618,4 @@
 
 		#bedinfo and surface info
-		md.mesh.elementonbase=project3d(md,'vector',numpy.ones(md.mesh.numberofelements2d,bool),'type','element','layer',1)
-		md.mesh.elementonsurface=project3d(md,'vector',numpy.ones(md.mesh.numberofelements2d,bool),'type','element','layer',md.mesh.numberoflayers-1)
 		md.mesh.vertexonbase=project3d(md,'vector',numpy.ones(md.mesh.numberofvertices2d,bool),'type','node','layer',1)
 		md.mesh.vertexonsurface=project3d(md,'vector',numpy.ones(md.mesh.numberofvertices2d,bool),'type','node','layer',md.mesh.numberoflayers)
Index: /issm/trunk-jpl/src/m/consistency/checkfield.m
===================================================================
--- /issm/trunk-jpl/src/m/consistency/checkfield.m	(revision 18076)
+++ /issm/trunk-jpl/src/m/consistency/checkfield.m	(revision 18077)
@@ -23,7 +23,4 @@
 %   Usage:
 %      md = checkfield(md,fieldname,options);
-%
-%   Example:
-%      md = checkfield(md,'mesh.elementonbase','size',[md.mesh.numberofelements 1],'values',[0 1]);
 
 %get options
Index: /issm/trunk-jpl/src/m/consistency/checkfield.py
===================================================================
--- /issm/trunk-jpl/src/m/consistency/checkfield.py	(revision 18076)
+++ /issm/trunk-jpl/src/m/consistency/checkfield.py	(revision 18077)
@@ -30,7 +30,4 @@
 	   Usage:
 	      md = checkfield(md,fieldname,options);
-
-	   Example:
-	      md = checkfield(md,'mesh.elementonbase','size',[md.mesh.numberofelements 1],'values',[0 1]);
 	"""
 
Index: /issm/trunk-jpl/src/m/contrib/bamg/yams.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/bamg/yams.m	(revision 18076)
+++ /issm/trunk-jpl/src/m/contrib/bamg/yams.m	(revision 18077)
@@ -108,6 +108,4 @@
 md.mesh.vertexonbase=ones(md.mesh.numberofvertices,1);
 md.mesh.vertexonsurface=ones(md.mesh.numberofvertices,1);
-md.mesh.elementonbase=ones(md.mesh.numberofelements,1);
-md.mesh.elementonsurface=ones(md.mesh.numberofelements,1);
 if strcmpi(Names.interp,'node'),
 	md.inversion.vx_obs=InterpFromGridToMesh(Vel.(Names.xname),Vel.(Names.yname),Vel.(Names.vxname),md.mesh.x,md.mesh.y,0);