%This is file crystalplot.m   
%plots and labels the 11 crystals we care about
%must be called after setup.m has been run. 

i=1; % Defines which time frame you want to plot

% Now plot the crystals with dots and lines between some crystals
plot3(...
        x(i,[1 4 2 3 1]),y(i,[1 4 2 3 1]),z(i,[1 4 2 3 1]),'.r-',... % M1 SW M2 AW, central
        x(i,[7 8 12 13 11 7]),y(i,[7 8 12 13 11 7]),z(i,[7 8 12 13 11 7]),'.b-',... % mitral annulus
        x(i,[14]),y(i,[14]),z(i,[14]),'c.',... % base crystal
        x(i,16),y(i,16),z(i,16),'r.',...% apex       
        x(i,[14 4 16]),y(i,[14 4 16]),z(i,[14 4 16]),'c-',... % base->SW->Apex
        x(i,[14 7]),y(i,[14 7]),z(i,[14 7]),'c-',... %base -> T1 
        x(i,[14 8]),y(i,[14 8]),z(i,[14 8]),'c-',... %base -> T2
        x(i,[13 3 16]),y(i,[13 3 16]),z(i,[13 3 16]),'c-',... % PM->AW->Apex
        'MarkerSize',10);
       
% Now add text labels and bounds
    text(x(i,1),y(i,1),z(i,1),'M1')
    text(x(i,2),y(i,2),z(i,2),'M2')
    text(x(i,3),y(i,3),z(i,3),'AW')
    text(x(i,4),y(i,4),z(i,4),'SW')
%   text(x(i,5),y(i,5),z(i,5),'S1')  % excluded
%   text(x(i,6),y(i,6),z(i,6),'S2')  % excluded
    text(x(i,7),y(i,7),z(i,7),'T1')
    text(x(i,8),y(i,8),z(i,8),'T2')
%   text(x(i,9),y(i,9),z(i,9),'AL')     % excluded
%   text(x(i,10),y(i,10),z(i,10),'PL')  % excluded
    text(x(i,11),y(i,11),z(i,11),'P1')
    text(x(i,12),y(i,12),z(i,12),'P2')
    text(x(i,13),y(i,13),z(i,13),'PM')
    text(x(i,14),y(i,14),z(i,14),'base')
%   text(x(i,15),y(i,15),z(i,15),'com')  % excluded
    text(x(i,16),y(i,16),z(i,16),'apex')
    grid on
    axis equal
    xlabel('x')
    ylabel('y')
    zlabel('z')
    title(['11 crystals'])

% After you run this, you can click on the circle tab at the top of the graphics window
%    and rotate the graph with the mouse.