function plotLongitudinalModulusComparison(ctrl,determinedModulus)

load(ctrl.indentationResultFile)
micromechanicalTesting = determinedModulus;
coxInverted = 15.2;


figure('color','w','units','centimeters','OuterPosition',[10 10 16 8]);

% boxchart([4.*ones(size(micromechanicalTesting))],[micromechanicalTesting],'Orientation','horizontal','handlevisibility','off',...
%          'BoxFaceColor',[0 0 0],'MarkerColor',[0 0 0])
boxchart([4.*ones(size(micromechanicalTesting))],[micromechanicalTesting],'Orientation','horizontal','handlevisibility','off',...
        'BoxFaceColor',ctrl.colors(1,:),'MarkerColor',ctrl.colors(1,:),'WhiskerLineColor',ctrl.colors(1,:))
hold on
% plot(mean(micromechanicalTesting,'omitnan'), 4,'k+','handlevisibility','off','linewidth',1.5)
plot(mean(micromechanicalTesting,'omitnan'), 4,'+','color',ctrl.colors(1,:),'handlevisibility','off','linewidth',1.5)

% plot([finalState.pyrMean(1) finalState.pyrMedian(1)]  ,1.*[1 1],'k-','MarkerFaceColor','k','MarkerSize',8,'handlevisibility','off')
% plot([finalState.hemiMean(1) finalState.hemiMedian(1)],2.*[1 1],'k-','MarkerFaceColor','k','MarkerSize',8,'handlevisibility','off')
% plot([finalState.nanoMean(1) finalState.nanoMedian(1)],3.*[1 1],'k-','MarkerFaceColor','k','MarkerSize',8,'handlevisibility','off')
plot([finalState.pyrMean(1) finalState.pyrMedian(1)]  ,1.*[1 1],'-','MarkerFaceColor',ctrl.colors(1,:),'color',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')
plot([finalState.hemiMean(1) finalState.hemiMedian(1)],2.*[1 1],'-','MarkerFaceColor',ctrl.colors(1,:),'color',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')
plot([finalState.nanoMean(1) finalState.nanoMedian(1)],3.*[1 1],'-','MarkerFaceColor',ctrl.colors(1,:),'color',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')

% plot(finalState.pyrMean(1),1,'wo','MarkerFaceColor','k','MarkerSize',8,'displayname','Mean M_L used in Eq. (16)')
% plot(finalState.hemiMean(1),2,'wo','MarkerFaceColor','k','MarkerSize',8,'handlevisibility','off')
% plot(finalState.nanoMean(1),3,'wo','MarkerFaceColor','k','MarkerSize',8,'handlevisibility','off')  
plot(finalState.pyrMean(1),1,'wo','MarkerFaceColor',ctrl.colors(1,:),'MarkerSize',8,'displayname','Mean M_L used in Eq. (16)')
plot(finalState.hemiMean(1),2,'wo','MarkerFaceColor',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')
plot(finalState.nanoMean(1),3,'wo','MarkerFaceColor',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')  

% plot(finalState.pyrMedian(1),1,'ws','MarkerFaceColor',0.*[0.5 0.5 0.5],'MarkerSize',8,'displayname','Median M_L used in Eq. (16)')
% plot(finalState.hemiMedian(1),2,'ws','MarkerFaceColor',0.*[0.5 0.5 0.5],'MarkerSize',8,'handlevisibility','off')
% plot(finalState.nanoMedian(1),3,'ws','MarkerFaceColor',0.*[0.5 0.5 0.5],'MarkerSize',8,'handlevisibility','off')  
plot(finalState.pyrMedian(1),1,'ws','MarkerFaceColor',ctrl.colors(1,:),'MarkerSize',8,'displayname','Median M_L used in Eq. (16)')
plot(finalState.hemiMedian(1),2,'ws','MarkerFaceColor',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')
plot(finalState.nanoMedian(1),3,'ws','MarkerFaceColor',ctrl.colors(1,:),'MarkerSize',8,'handlevisibility','off')  
    
% plot(coxInverted.*[1],5,'dk','linewidth',1,'markerfacecolor','k','handlevisibility','off')
plot(coxInverted.*[1],5,'d','color',ctrl.colors(1,:),'linewidth',1,'markerfacecolor',ctrl.colors(1,:),'handlevisibility','off')


yticks([1:5])
yticklabels({'AFM-NI, Pyramid','AFM-NI, Hemisphere','NI','Single fiber tensile test',"Cox' theory"})
ylim([0 6])

legend('location','southeast')
xlabel('E_L [GPa]')
xlim([0 35])
set(gca,ctrl.axisInstructions{:})
box on

print([ctrl.plotDir 'figure8_ELComparison'],'-dpng','-r1200')
print([ctrl.plotDir 'figure8_ELComparison'],'-depsc')