conditions = {'KCSC','KCSP','KPSC','KPSP'};

selectedCondition = 1;

c = conditions{selectedCondition};

cd(c);

filenames = dir('*.mat');

hF = figure('Position',[441,94,964,999],'Color',[1,1,1]);
hold on

axis equal
view([158,-24]);

toPlot = [1,4,8,9,11,12,13,16,17,18,19];

CS_colors = [0,1,0;...      % CS 1
             0.3,0.7,1;...  % CS 3
             0.4,1,0.4;...  % CS 5
             1,0.4,0.4;...  % CS 9
             1,0.2,0.2;...  % CS 8
             0,1,1;...      % CS 7
             0.6,1,0.6;...  % CS 6
             1,0.8,0.8;...  % CS 11
             0.1,0.7,1;...  % CS 2
             0.2,1,0.2;...  % CS 4
             1,0.4,0.4];    % CS 10 
j = 1;
      
for iFile = toPlot
    
    data = [];
    
    load(filenames(iFile).name);
    
    displacement = data(:,1:3)*50 + data(:,4:6);
    
    data = data(:,4:6);
    
    k = boundary(data);
    hSurf = trisurf(k,data(:,1),data(:,2),data(:,3),...
        'Facecolor',CS_colors(j,:),...
        'FaceAlpha',0.3,...
        'EdgeColor','none');
    
    set(hSurf,'FaceLighting','gouraud')
    set(hSurf,'EdgeLighting','gouraud')
    set(hSurf,'SpecularStrength',0.3)
    
    xlabel('x');
    ylabel('y');
    zlabel('z');
    
    hArrow = arrow3(data(1:200:end,:),displacement(1:200:end,:),'-k1',1);
    
    j = j + 1;
    
end

ax = gca;

xTicks = {'-8','','','','0','','','','8',''};
yTicks = {'-4','','0','','4'};
zTicks = {'-8','','','','0','','','','8',''};

set(ax,'XTickLabel',xTicks);
set(ax,'YTickLabel',yTicks);
set(ax,'ZTickLabel',zTicks);

xlabel('µm');
ylabel('µm');
zlabel('µm');

axis equal
view([158,-24]);

box off

camlight('headlight');

a = gca;

a.XGrid = 'off';
a.YGrid = 'off';
a.ZGrid = 'off';

sgtitle(c,'Interpreter','none');

set(hF,'Renderer','painters');