Index: /issm/trunk-jpl/src/m/plot/plot_icefront.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_icefront.py	(revision 21610)
+++ /issm/trunk-jpl/src/m/plot/plot_icefront.py	(revision 21611)
@@ -22,6 +22,4 @@
 	noice=np.where(np.sum(md.mask.ice_levelset[elements],1)==3)
 
-	#hydro neumann
-	hydro_neumann=np.where(md.hydrology.neumannflux!=0)
 	#plot mesh
 	ax.triplot(x,y,elements)
@@ -31,7 +29,4 @@
 		colors=np.asarray([0.5 for element in elements[icefront]])
 		ax.tripcolor(x,y,elements[icefront],facecolors=colors,alpha=0.5,label='elements on ice front')
-	if len(hydro_neumann[0])>0:
-		colors=np.asarray([0.5 for element in elements[hydro_neumann]])
-		ax.tripcolor(x,y,elements[hydro_neumann],facecolors=colors,alpha=0.5,label='non zero neumann flux for the hydrology')
 
 	#apply options
Index: /issm/trunk-jpl/src/m/plot/plot_unit.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_unit.py	(revision 21610)
+++ /issm/trunk-jpl/src/m/plot/plot_unit.py	(revision 21611)
@@ -147,5 +147,5 @@
 	elif datatype==2:
 		if is2d:
-			if options.exist('mask'):
+			if np.ma.is_masked(data):
 				EltMask=np.asarray([np.any(np.in1d(index,np.where(data.mask))) for index in elements])
 				triangles=mpl.tri.Triangulation(x,y,elements,EltMask)
@@ -193,5 +193,4 @@
 			ax.set_ylim([min(y),max(y)])
 			ax.set_zlim([min(z),max(z)])
-
 			#raise ValueError('plot_unit error: 3D element plot not supported yet')
 		return