function tetraeder_explosiv(Tets,KK, shiftScale)
% tetraeder_explosiv(Tets,KK, shiftScale)
%
% The tetrahedrons are shown in an exploded view
% INPUT:    Test[nelem,4] - element connectivity of the tets
%           KK[3,nnode]   - vertex coordinates
%           shiftScale    - expansion factor (0 == do nothing)

% Tetraeder in Explosionsdarstellung
nelem = size(Tets,1);        % Anzahl der Tetraeder
nnode = size(KK,2);          % Anzahl der Eckpunkte
MM  = mean(KK);              % Mittelpunkt aller Eckpunkte

for k=1:nelem
    id = Tets(k,:);          % Eckpunkte aktuelles Tetraeder
    XYZ = KK(id,:);          % Koordinaten aktuelles Tetraeder
    mid = mean(XYZ);         % Mittelpunkt aktuelles Tetraeder
    shift = (mid-MM)*shiftScale;    % Explosionsverschiebung 
    XYZ = XYZ+shift;
    %XYZ = XYZ+repmat(shift,4,1);
    elemk = [1 2 3; 1 2 4; 1 3 4; 2 3 4 ]; % Oberflaechen des Tets in lokaler Numerierung
    
    axis manual
    trisurf(elemk,XYZ(:,1),XYZ(:,2),XYZ(:,3),k)
    hold on
end
% xlabel('x'),ylabel('y'),zlabel('z')
% title(['Explosionsdarstellung: ',num2str(nnode),' Punkte fuer ',num2str(nelem),' Tetraeder'])
% view([-14,42])
% alpha(0.7)
% axis equal

end