Changeset 24762

Timestamp:

04/30/20 08:01:46 (5 years ago)

Author:

Mathieu Morlighem

Message:

NEW: better handling of plotmodel with levelset: use 'gridded',1,'levelset',md.results.TransientSolution(end).MaskIceLevelset

File:

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

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

@@ -10 +10 @@
 [x y z elements is2d isplanet]=processmesh(md,[],options);
 [data datatype]=processdata(md,data,options);
+
+islevelset = exist(options,'levelset');
+if islevelset
+        levelset = getfieldvalue(options,'levelset');
+        [levelset datatype]=processdata(md,levelset,options);
+end
 
 %check is2d
@@ -29 +35 @@
 if size(data_grid,1)<3 | size(data_grid,2)<3,
         error('data_grid size too small in plot_gridded, check posting and units');
+end
+
+%Mask values if levelset>0
+if islevelset
+        ls_grid=InterpFromMeshToGrid(elements,x,y,levelset,x_m,y_m,NaN);
+        data_grid(ls_grid>0) = NaN;
 end