Changeset 2518

Timestamp:

10/23/09 10:33:15 (15 years ago)

Author:

Mathieu Morlighem

Message:

minor

File:

• issm/trunk/src/m/classes/public/contourenvelope.m (modified) (2 diffs)

issm/trunk/src/m/classes/public/contourenvelope.m

@@ -1 @@
-function segments=contourenvelope(md,file)
+function segments=contourenvelope(md,varargin)
 %CONTOURENVELOPE - build a set of segments enveloping a contour .exp
 %
 %   Usage:
-%      segments=contourenvelope(md,file)
+%      segments=contourenvelope(md,varargin)
+%
+%   Example:
+%      segments=contourenvelope(md,'Stream.exp');
+%      segments=contourenvelope(md);
 
 %some checks
-if ~exist(file),
-        error(['thicknessevolution error message: file ' file ' not found']);
+if nargin>2,
+        help contourenvelope
+        error('contourenvelope error message: bad usage');
+end
+if nargin==2,
+        file=varargin{1};
+        if ~exist(file),
+                error(['thicknessevolution error message: file ' file ' not found']);
+        end
 end
 
@@ -20 +31 @@
 
 %get nodes inside profile
-%get flag list of elements and nodes inside the contour
-nodein=ContourToMesh(md.elements,md.x,md.y,expread(file,1),'node',1);
-elemin=(sum(nodein(md.elements),2)==size(md.elements,2));
-%modify element connectivity
-elemout=find(~elemin);
 elementconnectivity=md.elementconnectivity;
-elementconnectivity(elemout,:)=0;
-elementconnectivity(find(ismember(elementconnectivity,elemout)))=0;
+if nargin==2,
+        %get flag list of elements and nodes inside the contour
+        nodein=ContourToMesh(md.elements,md.x,md.y,expread(file,1),'node',1);
+        elemin=(sum(nodein(md.elements),2)==size(md.elements,2));
+        %modify element connectivity
+        elemout=find(~elemin);
+        elementconnectivity(elemout,:)=0;
+        elementconnectivity(find(ismember(elementconnectivity,elemout)))=0;
+end
 
 %Find element on boundary
 %First: find elements on the boundary of the domain
 flag=elementconnectivity;
-flag(find(flag))=elemin(flag(find(flag)));
+if nargin==2,
+        flag(find(flag))=elemin(flag(find(flag)));
+end
 elementonboundary=double(sum(flag,2)<3 & sum(flag,2)>0);