source: issm/trunk/src/m/utils/Geometry/FlagElements.m@ 2967

function flag=FlagElements(md,area),
%FLAGELEMENTS - flag the elements in an area
%
%   The area can be given with an exp file, a list of elements.
%
%   Usage: 
%      flag=FlagElements(md,area);
%
%   Example:
%      flag=FlagElements(md,'all');
%      flag=FlagElements(md,'');
%      flag=FlagElements(md,'Domain.exp');
%      flag=FlagElements(md,'~Domain.exp');
%      flag=FlagElements(md,md.elementonicesheet);

        if ischar(area),
                if isempty(area),
                        flag=zeros(md.numberofelements,1);
                        invert=0;
                elseif strcmpi(area,'all')
                        flag=ones(md.numberofelements,1);
                        invert=0;
                else
                        %make sure that we actually don't want the elements outside the domain outline!
                        if strcmpi(area(1),'~'),
                                area=area(2:length(area));
                                invert=1;
                        else
                                invert=0;
                        end
                        %ok, flag elements
                        flag=ContourToMesh(md.elements(:,1:3),md.x,md.y,expread(area,1),'element',1);
                end
                if invert,
                        flag=~flag;
                end
        elseif isfloat(area),
                if size(area,1)~=md.numberofelements,
                        setelementstypeusage();
                        error('Flaglist for area must be of same size as number of elements in model');
                end
                flag=area;
        else
                error('Invalid area option option');
        end
end