Index: /issm/trunk-jpl/src/m/plot/plot_overlay.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_overlay.m	(revision 25796)
+++ /issm/trunk-jpl/src/m/plot/plot_overlay.m	(revision 25797)
@@ -16,4 +16,10 @@
 	radaronly=0;
 	[data datatype]=processdata(md,data,options);
+end
+
+islevelset = exist(options,'levelset');
+if islevelset
+	levelset = getfieldvalue(options,'levelset');
+	[levelset datatype]=processdata(md,levelset,options);
 end
 
@@ -68,4 +74,9 @@
 else
 	data_grid=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),data,x_m,y_m,NaN);
+end
+%Mask values if levelset>0
+if islevelset
+	ls_grid=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),levelset,x_m,y_m,+1);
+	data_grid(ls_grid>0) = NaN;
 end