Changeset 1239

Timestamp:

07/06/09 14:52:27 (16 years ago)

Author:

Mathieu Morlighem

Message:

extended BasinConstrain for ~domain

File:

• issm/trunk/src/m/classes/public/BasinConstrain.m (modified) (1 diff)

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

@@ -12 +12 @@
 %   Example:
 %      md=BasinConstrain(md,'DomainOutline.exp');
+%      md=BasinConstrain(md,'~Iceshelves.exp');
 
 %now, flag grids and elements outside the domain outline.
-[gridondomain elementondomain]=ContourToMesh(md.elements,md.x,md.y,expread(domain,1),'element and node',2);
+if ischar(domain),
+        if isempty(domain),
+                elementondomain=zeros(md.numberofelements,1);
+                gridondomain=zeros(md.numberofgrids,1);
+                invert=0;
+        elseif strcmpi(domain,'all')
+                elementondomain=ones(md.numberofelements,1);
+                gridondomain=ones(md.numberofgrids,1);
+                invert=0;
+        else
+                %make sure that we actually don't want the elements outside the domain outline!
+                if strcmpi(domain(1),'~'),
+                        domain=domain(2:length(domain));
+                        invert=1;
+                else
+                        invert=0;
+                end
+                %ok, flag elements and nodes
+                [gridondomain elementondomain]=ContourToMesh(md.elements,md.x,md.y,expread(domain,1),'element and node',2);
+        end
+        if invert,
+                gridondomain=~gridondomain;
+                elementondomain=~elementondomain;
+        end
+else
+        error('BasinConstrain error message: domain type not supported yet');
+end
+
+%list of elements and nodes not on domain
 gridnotondomain=find(~gridondomain);
 elementnotondomain=find(~elementondomain);