Index: /issm/trunk-jpl/src/m/plot/plot_elementnumbering.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_elementnumbering.py	(revision 21474)
+++ /issm/trunk-jpl/src/m/plot/plot_elementnumbering.py	(revision 21475)
@@ -9,5 +9,5 @@
 from plot_icefront import plot_icefront
 
-def plot_elementnumbering(md,options,fig,ax):
+def plot_elementnumbering(md,options,fig,axgrid,gridindex):
 	'''
 	plot_elementnumbering - plot element numberign (starting at 1 matlab and c convention)
@@ -18,6 +18,14 @@
 		See also: PLOTMODEL
 	'''
+	x,y,z,elements,is2d,isplanet=processmesh(md,[],options)
+
+	ax=axgrid[gridindex]
+	fig.delaxes(axgrid.cbar_axes[gridindex])
 	
-	x,y,z,elements,is2d,isplanet=processmesh(md,[],options)
+	if is2d:
+		ax.triplot(x,y,elements)
+	else:
+		ax=inset_locator.inset_axes(axgrid[gridindex],width='100%',height='100%',loc=3,borderpad=0,axes_class=Axes3D)
+
 	XLims=[np.min(x),np.max(x)]
 	YLims=[np.min(y),np.max(y)]
@@ -37,8 +45,7 @@
 		ax.text(np.mean(x[element]),np.mean(y[element]),str(i+1),ha='center',va='center')
 		
-
 	#apply options
-	options.addfielddefault('title','Element numbers')
+	options.addfielddefault('title','Element numbers (matlab indexation)')
 	options.addfielddefault('colorbar','off')
-	applyoptions(md,[],options,fig,ax)
+	applyoptions(md,[],options,fig,axgrid,gridindex)
 	
Index: /issm/trunk-jpl/src/m/plot/plot_manager.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_manager.py	(revision 21474)
+++ /issm/trunk-jpl/src/m/plot/plot_manager.py	(revision 21475)
@@ -67,15 +67,15 @@
 			return
 		elif data=='BC': 
-			plot_BC(md,options,fig,ax)
+			plot_BC(md,options,fig,axgrid,gridindex)
 			return
 		elif data=='elementnumbering': 
-			plot_elementnumbering(md,options,fig,ax)
+			plot_elementnumbering(md,options,fig,axgrid,gridindex)
 			return
 		elif data=='vertexnumbering': 
-			plot_vertexnumbering(md,options,fig,ax)
+			plot_vertexnumbering(md,options,fig,axgrid,gridindex)
 			return
 		elif data=='none':
 			print 'no data provided to plot (TODO: write plot_none.py)'
-			applyoptions(md,[],options,fig,ax)
+			applyoptions(md,[],options,fig,axgrid,gridindex)
 			return
 		else:
Index: /issm/trunk-jpl/src/m/plot/plot_mesh.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_mesh.py	(revision 21474)
+++ /issm/trunk-jpl/src/m/plot/plot_mesh.py	(revision 21475)
@@ -27,9 +27,5 @@
 	if is2d:
 		ax.triplot(x,y,elements)
-
 	else:
-		# fig.delaxes(ax)
-		# geom=int(str(axgrid.get_geometry()[0])+str(axgrid.get_geometry()[1])+str(gridindex+1))
-		# ax = fig.add_subplot(geom,projection='3d')
 		ax=inset_locator.inset_axes(axgrid[gridindex],width='100%',height='100%',loc=3,borderpad=0,axes_class=Axes3D)
 		
Index: /issm/trunk-jpl/src/m/plot/plot_vertexnumbering.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_vertexnumbering.py	(revision 21474)
+++ /issm/trunk-jpl/src/m/plot/plot_vertexnumbering.py	(revision 21475)
@@ -1,5 +1,2 @@
-
-
-
 try:
 	import pylab as p
@@ -12,37 +9,45 @@
 from plot_icefront import plot_icefront
 
-def plot_vertexnumbering(md,options,fig,ax):
-    '''
+def plot_vertexnumbering(md,options,fig,axgrid,gridindex):
+	'''
 	PLOT_VERTEXNUMBERING - plot vertex numbering
-
-   	Usage:
-      	plot_vertexnumbering(md,options,fig,axes);
-
+	
+	Usage:
+	plot_vertexnumbering(md,options,fig,axes);
+	
  	See also: PLOTMODEL
-
+	
 	'''
 	#process data and model
-    x,y,z,elements,is2d,isplanet=processmesh(md,[],options)
-    XLims=[np.min(x),np.max(x)]
-    YLims=[np.min(y),np.max(y)]
-    #plot mesh
-    ax.triplot(x,y,elements)
-    highlightpos=options.getfieldvalue('highlight','none')
-    if highlightpos!='none':
-        #if just one element duplicate it to avoid coloring issues
-        if type(highlightpos)==int:
-            highlightpos=[highlightpos,highlightpos]
-        #convert from to matlab numbering
-        highlightpos=[pos-1 for pos in highlightpos]
-        colors=np.asarray([0.5 for element in elements[highlightpos]])
-        ax.tripcolor(x,y,vertices[highlightpos],facecolors=colors,alpha=0.5)
-    # and numbers
+	x,y,z,elements,is2d,isplanet=processmesh(md,[],options)
+
+	ax=axgrid[gridindex]
+	fig.delaxes(axgrid.cbar_axes[gridindex])
 	
-    for i,element in enumerate(elements):
-        for j in [0,1,2]:
-            ax.text(x[element][j],y[element][j],str(element[j]))	
+	if is2d:
+		ax.triplot(x,y,elements)
+	else:
+		ax=inset_locator.inset_axes(axgrid[gridindex],width='100%',height='100%',loc=3,borderpad=0,axes_class=Axes3D)
+
+
+	XLims=[np.min(x),np.max(x)]
+	YLims=[np.min(y),np.max(y)]
+	#plot mesh
+	ax.triplot(x,y,elements)
+	highlightpos=options.getfieldvalue('highlight','none')
+	if highlightpos!='none':
+		#if just one element duplicate it to avoid coloring issues
+		if type(highlightpos)==int:
+			highlightpos=[highlightpos,highlightpos]
+		#convert from to matlab numbering
+		highlightpos=[pos-1 for pos in highlightpos]
+		colors=np.asarray([0.5 for element in elements[highlightpos]])
+		ax.tripcolor(x,y,vertices[highlightpos],facecolors=colors,alpha=0.5)
+	# and numbers
+	for i,Xcoord in enumerate(x):
+		ax.text(x[i],y[i],str(i+1))	
 	
-    #apply options
-    options.addfielddefault('title','Vertex numbers')
-    options.addfielddefault('colorbar','off')
-    applyoptions(md,[],options,fig,ax)
+	#apply options
+	options.addfielddefault('title','Vertex numbers (matlab indexation)')
+	options.addfielddefault('colorbar','off')
+	applyoptions(md,[],options,fig,axgrid,gridindex)