Changeset 25797

Timestamp:

12/01/20 13:30:52 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added levelset option to plot_overlay.m

File:

• issm/trunk-jpl/src/m/plot/plot_overlay.m (modified) (2 diffs)

issm/trunk-jpl/src/m/plot/plot_overlay.m

@@ -16 +16 @@
         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 +74 @@
 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