function [loop] = get_Solid_El_Loop (nodes)
%   ------------------------------------------------------
%  Get node loop to draw an element by closing to first node.
%          Draw any interior nodes at the first node.
%   ------------------------------------------------------
% loop   = corners for nodes line polygon

if ( nodes == 8 )
  loop = [ 1:4, 1, 5:8, 5, 8, 4, 3, 7, 1]; % default for H8
%           4   5    9  10 11 12 13 14
elseif (nodes == 6)
  loop = [ 1, 2, 3, 1, 4, 2, 3, 4, 1 ]
end % if
% end function get_Solid_El_Loop