clear all; close all; clc;
%
% jfm 2018
% figure 14  
% data structure [RT KL_normalised]
%
load data_fig14_KLvsRT.mat
%
% free-slip data from hw18 and hw14
% no-slip data from hw16
loglog(data.hw18(:,1), data.hw18(:,2),...
    'ok','MarkerFaceColor','none','MarkerSize',4); 
hold on
loglog(data.hw14(:,1), data.hw14(:,2),...
    'sk','MarkerFace','none','MarkerSize',4);
loglog(data.hw16(:,1), data.hw16(:,2),...
     '^k','MarkerFace','none','MarkerSize',4);
%
% fitted lines
pslope05fs=loglog([30 500], 1.6*([30 500].^-0.5),'--k','LineWidth',2);
pslope025fs=loglog([300 8000], 0.3425*[300 8000].^-0.25,'-k','LineWidth',2);
pslope05=loglog([30 500], 0.945*[30 500].^-0.5,':k','LineWidth',1);
pslope025=loglog([300 8000], 0.1975*[300 8000].^-0.25,'-.k','LineWidth',1);
%
legend({'present DNS, n=1/2', ...
    'free-slip (HW14), n=1/2',...
    'no-slip (HW16), n=2/3',...
    'eq(4.8), c=1.6, n=1/2, q=1/2',...
    'eq(4.8), c=0.3425, n=1/2, q=1/4 ',...
    'eq(4.8), c=0.945, n=2/3, q=1/2',...
    'eq(4.8), c=0.1975, n=2/3, q=1/4'})
%
xlabel('R_T');
ylabel('$$\mathrm {\frac{K_L}{  u_\infty Sc^{-n}}}$$',...
    'interpreter','latex');
xlim([50 3000]);
ylim([0.01 0.5]);
grid on;
set(legend,'Box','off')