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