root/gliderproc/trunk/MATLAB/opnml/FEM/drawelems3d.m

% DRAWELEMS3D draw 2-D FEM element configuration in 3-D
%
% DRAWELEMS3D draws element boundries given a valid grid structure.  
%
%  INPUT : fem_grid_struct - (from LOADGRID, see FEM_GRID_STRUCT)       
%           
% OUTPUT : hel - handle to the element object.
%
%   CALL : hel=drawelems3d(fem_grid_struct);
%
% Written by: Brian O. Blanton
% Summer 1998
%                  
function hel=drawelems3d(fem_grid_struct) 

% DEFINE ERROR STRINGS
err1=['Not enough input arguments; need a fem_grid_struct'];

% check arguments
if nargin ==0 
   error(err1);
end  

if ~is_valid_struct(fem_grid_struct)
   error('    Argument to DRAWELEMS must be a valid fem_grid_struct.')
end

% Extract grid fields from fem_grid_struct
%
elems=fem_grid_struct.e;
% COPY FIRST COLUMN TO LAST TO CLOSE ELEMENTS
%
elems=elems(:,[1 2 3 1]);

x=fem_grid_struct.x;
y=fem_grid_struct.y;
z=fem_grid_struct.z;

elems=elems';
[m,n]=size(elems);
xt=x(elems);
yt=y(elems);
zt=z(elems);
if n~=1 
   if m>n
      xt=reshape(xt,n,m);
      yt=reshape(yt,n,m);
      zt=reshape(zt,n,m);
   else
      xt=reshape(xt,m,n);
      yt=reshape(yt,m,n);
      zt=reshape(zt,m,n);
   end
   xt=[xt
       NaN*ones(size(1:length(xt)))];
   yt=[yt
       NaN*ones(size(1:length(yt)))];
   zt=[zt
       NaN*ones(size(1:length(zt)))];
end
xt=xt(:);
yt=yt(:);
zt=zt(:);
%
% DRAW GRID
%
hel=line(xt,yt,zt,'LineWidth',1,'LineStyle','-','Color',[.8 .8 .8]);
set(hel,'Tag','elements');
 
%
%        Brian O. Blanton
%        Curr. in Marine Sciences
%        15-1A Venable Hall
%        CB# 3300
%        Uni. of North Carolina
%        Chapel Hill, NC
%                 27599-3300
%
%        919-962-4466
%        blanton@cuda.chem.unc.edu
%
%        Summer 1997
%