function hm=plottetra(node,elem,varargin)
%
% hm=plottetra(node,elem,opt)
%
% plot 3D surface meshes
% 
% author: Qianqian Fang <fangq at nmr.mgh.harvard.edu>
%
% input: 
%      node: a node coordinate list, 3 columns for x/y/z; if node has a 
%            4th column, it will be used to set the color at each node.
%      elem: a tetrahedral element list; if elem has a 5th column,
%            it will be used to separate the mesh into 
%            sub-domains and display them in different colors.
%      opt:  additional options for a patch object, see plotmesh
%
% output:
%   hm: handle or handles (vector) to the plotted surfaces
%
% example:
%
%   h=plottetra(node,elem);
%   h=plottetra(node,elem,'facealpha',0.5);
% 
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%

rngstate = rand ('state');

randseed=hex2dec('623F9A9E'); % "U+623F U+9A9E"
if(isoctavemesh) randseed=randseed+3; end

if(~isempty(getvarfrom('base','ISO2MESH_RANDSEED')))
        randseed=getvarfrom('base','ISO2MESH_RANDSEED');
end
rand('state',randseed);

if(~iscell(elem))
	if(size(elem,2)>4)
		tag=elem(:,5);
		types=unique(tag);
		hold on;
		h=[];
		for i=1:length(types)
			idx=find(tag==types(i));
			face=volface(elem(idx,1:4));
	                if(size(node,2)==3)
			    h=[h plotsurf(node,face,'facecolor',rand(3,1),varargin{:})];
			else
                            h=[h plotsurf(node,face,varargin{:})];
			end
		end
	else
		face=volface(elem(:,1:4));
		h=plotsurf(node,face,varargin{:});
	end
end

if(~isempty(h)) 
  axis equal;
end
if(~isempty(h) & nargout>=1)
  hm=h;
end

rand ('state',rngstate);