Changeset 26792 for issm/trunk-jpl/src/m/exp/contourlevelzero.m

Timestamp:

01/19/22 06:58:30 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added handy output for contourlevelzero

File:

• issm/trunk-jpl/src/m/exp/contourlevelzero.m (modified) (3 diffs)

issm/trunk-jpl/src/m/exp/contourlevelzero.m

@@ -1 @@
-function contours=contourlevelzero(md,mask,level)
+function contours=contourlevelzero(md,mask,level,varargin)
 %CONTOURLEVELZERO - figure out the zero level (or offset thereof, specified by the level value) 
 %                   of a vectorial mask, and vectorialize it into an exp or shp compatible 
@@ -7 +7 @@
 %      contours=contourlevelzero(md,mask,level)
 %
+%   Example:
+%      contours=contourlevelzero(md, md.results.TransientSolution(end).MaskOceanLevelset,0);
+%      contours=contourlevelzero(md, md.results.TransientSolution(end).MaskOceanLevelset,0,'output','vector');
+%
+%   Supported options:
+%      'output': either 'matrix' (with NaN delimiters) or 'struct' (default)
+%
 %   See also: PLOT_CONTOUR
+
+%Process options
+options = pairoptions(varargin{:});
 
 %process data 
@@ -228 +238 @@
 
 end
+
+%process output
+if strcmp(getfieldvalue(options,'output',1),'matrix')
+        x = cell2mat(cellfun(@(x) [x;NaN],{contours(:).x},'UniformOutput',0)');
+        y = cell2mat(cellfun(@(x) [x;NaN],{contours(:).y},'UniformOutput',0)');
+
+        contours = [x y];
+end