function drawbars(Bars,Nodes,labflg,ax,rgb,style,nfontsz)
%
%  Name: drawbars
% 
%  Usage: drawbars(Bars,Nodes,labflg,ax,rgb,style,nfontsz)
% 
%   Description:
%      Draws lines to represent
%      the bar connectivity in 'Bars'
%      'Nodes' is the node list
%      'ax' is axis handle
%      'rgb' is color vector
%      'labflg' flag controls display of labels
%      'style' is linestyle string

%  Version SWD970710
%
%  Copyright (c) 1996-1997 by Scott W. Doebling
%	Los Alamos National Laboratory
%	ESA-EA, doebling@lanl.gov
%

if ~exist('nfontsz'), nfontsz=14; end

axes(ax)

for i = 1:size(Bars,1),

   n1 = Bars(i,2);
   n2 = Bars(i,3);

%  Find locations of end-nodes in Node list

   in1 = find(Nodes(:,1) == n1);
   in2 = find(Nodes(:,1) == n2);

   if isempty(in1) | isempty(in2), disp('ERROR: DRAWBARS -- Node Label is in Bars but not in Nodes'), end

   h = line([Nodes(in1,2) Nodes(in2,2)], [Nodes(in1,3) Nodes(in2,3)], [Nodes(in1,4) Nodes(in2,4)] );

   %h = plot3([Nodes(in1,2) Nodes(in2,2)], [Nodes(in1,3) Nodes(in2,3)], [Nodes(in1,4) Nodes(in2,4)] );

   %xvec((i-1)*2 + 1 : i*2) = [Nodes(in1,2) Nodes(in2,2)];
   %yvec((i-1)*2 + 1 : i*2) = [Nodes(in1,3) Nodes(in2,3)];
   %zvec((i-1)*2 + 1 : i*2) = [Nodes(in1,4) Nodes(in2,4)];

   set(h,'Color','y')
   set(h,'LineStyle',style)

   if labflg,
	tx = (Nodes(in2,2) + Nodes(in1,2))/2;
	ty = (Nodes(in2,3) + Nodes(in1,3))/2;
	tz = (Nodes(in2,4) + Nodes(in1,4))/2;
	ht = text(tx,ty,tz,num2str(Bars(i,1)));
	set(ht,'FontSize',nfontsz)
	set(ht,'FontWeight','bold')
	set(ht,'VerticalAlignment','bottom')
	set(ht,'Color',rgb)
   end
end

return
