Index: /issm/trunk-jpl/src/m/plot/plot_BC.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_BC.py	(revision 21655)
+++ /issm/trunk-jpl/src/m/plot/plot_BC.py	(revision 21656)
@@ -44,11 +44,13 @@
 							 y[np.where(~np.isnan(md.stressbalance.spcvz))],
 							 marker='o',c='y',s=80,label='vz Dirichlet',linewidth=0)
-		ax.scatter(x[np.where(~np.isnan(md.hydrology.spcepl_head))],
-							 y[np.where(~np.isnan(md.hydrology.spcepl_head))],
-							 marker='v',c='r',s=240,label='EPL Head',linewidth=0)
-		ax.scatter(x[np.where(~np.isnan(md.hydrology.spcsediment_head))],
-							 y[np.where(~np.isnan(md.hydrology.spcsediment_head))],
-							 marker='^',c='b',s=240,label='IDS head',linewidth=0)
-		
+		try:
+			ax.scatter(x[np.where(~np.isnan(md.hydrology.spcepl_head))],
+								 y[np.where(~np.isnan(md.hydrology.spcepl_head))],
+								 marker='v',c='r',s=240,label='EPL Head',linewidth=0)
+			ax.scatter(x[np.where(~np.isnan(md.hydrology.spcsediment_head))],
+								 y[np.where(~np.isnan(md.hydrology.spcsediment_head))],
+								 marker='^',c='b',s=240,label='IDS head',linewidth=0)
+		except AttributeError:
+			print ('Not treating Hydrologydc, skipping these boundaries')
 		ax.set_xlim(XLims)
 		ax.set_ylim(YLims)