% 3NN model correlation length

close all;

linewidth = 1;
smoothflag = true;
usebufflag = true; % use stored distribution bins or not
% set up figure
% labels = {'(a)', '(b)', '(c)', '(d)'};
widthpix = 320;
heightpix = 200;

leftmargin = 80;
bottommargin = 65;
topmargin = 20;
rightmargin = 30;
centerwmargin = 60;
centerhmargin = 60;
widthtotal = widthpix*2+leftmargin+centerwmargin+rightmargin;
heighttotal = heightpix*3+topmargin+centerhmargin*2+bottommargin;
left = leftmargin/widthtotal;
bottom = bottommargin/heighttotal;
width = widthpix/widthtotal;
height = heightpix/heighttotal;

centerW = centerwmargin/widthtotal;
centerH = centerhmargin/widthtotal;

psize = 26; % font size
msize = 8; % marker size

fig = figure('rend','painters','pos',[100 100 widthtotal heighttotal]);

%% panel (a): 0.15_0.3

panelA = axes('Position', [left, bottom+(height+centerH)*2, width, height]);

kappa = 0.15;
kappa2 = 0.3;
hs = 0;
kappastr = sprintf('%g_%g', kappa, kappa2);

sizesA = [6, 8, 10, 12, 14, 16];

n=numel(sizesA);
rpp=(1:n)';
corders = [(rpp-1)/(n-1), zeros(n,1), 1 - (rpp-1)/(n-1)];

plot(nan, nan);hold on;
for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/2DYBT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    if mod(asize,4)==0
        linespec = '-';
    else
        linespec = '--';
    end
    plot(1./data(:,3), log10( abs(data(:, 7)) /asize), linespec, 'color', corders(rp,:));
end
plot([0.48,0.48], [-2,2], '--k');
xlim([0.2, 1]);
ylim([log10(7e-2), log10(1e2) ]);

xticks([0.2,0.6,1]);
yticks([ -1, 0, 1, 2]);
set(gca, 'yticklabels', {' '; '10^0'; ' '; '10^2'});

xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];

set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (b): 0.3_0.6

panelB = axes('Position', [left+width+centerW, bottom+(height+centerH)*2, width, height]);

kappa = 0.3;
kappa2 = 0.6;
kappastr = sprintf('%g_%g', kappa, kappa2);

plot(nan, nan);hold on;
for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/2DYBT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    if mod(asize,4)==0
        linespec = '-';
    else
        linespec = '--';
    end
    plot(1./data(:,3), log10( abs(data(:, 7)) /asize), linespec, 'color', corders(rp,:));
end
plot([0.555,0.555], [-2,2], '--k');
xlim([0.2, 1]);
ylim([log10(7e-2), log10(1e2) ]);

xticks([0.2,0.6,1]);
yticks([ -1, 0, 1, 2]);
set(gca, 'yticklabels', {' '; '10^0'; ' '; '10^2'});

xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];

set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (c): 0.3_1

panelC = axes('Position', [left, bottom+(height+centerH), width, height]);

kappa = 0.3;
kappa2 = 1;
kappastr = sprintf('%g_%g', kappa, kappa2);

plot(nan, nan);hold on;
for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/2DYBT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    if mod(asize,4)==0
        linespec = '-';
    else
        linespec = '--';
    end
    plot(1./data(:,3), log10( abs(data(:, 7)) /asize), linespec, 'color', corders(rp,:));
end
plot([1.03,1.03], [-2,2], '--k');
xlim([0.5, 1.5]);
ylim([log10(7e-2), log10(1e2) ]);

xticks([0.5,1,1.5 ]);
yticks([ -1, 0, 1, 2]);
set(gca, 'yticklabels', {' '; '10^0'; ' '; '10^2'});

ylabel('$\xi_1/L$', 'interpreter', 'latex');
xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];

set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (d) 0.8_1.6, Y

panelD = axes('Position', [left+width+centerW, bottom+(height+centerH), width, height]);

kappa = 0.8;
kappa2 = 1.6;
kappastr = sprintf('%g_%g', kappa, kappa2);
sizesD = [8,12,16,20,24,28];
plot(nan, nan);hold on;
for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/2DYBT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    if mod(asize,4)==0
        linespec = '-';
    else
        linespec = '--';
    end
    plot(1./data(:,3), log10( abs(data(:, 7)) /asize), linespec, 'color', corders(rp,:));
end

xlim([0.2, 1.2]);
ylim([log10(7e-2), log10(1e2) ]);

xticks([0.2,0.7,1.2 ]);
yticks([ -1, 0, 1, 2]);
set(gca, 'yticklabels', {' '; '10^0'; ' '; '10^2'});

xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];

set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (e) 0.8_0.3, Z
panelE = axes('Position', [left, bottom, width, height]);

kappa = 0.8;
kappa2 = 0.3;
kappastr = sprintf('%g_%g', kappa, kappa2);

n=numel(sizesD);

plot(nan, nan);hold on;
for rp=1:n
    asize = sizesD(rp);
    fname = strcat('../data/2DZBT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    plot(1./data(:,3), log10( abs(data(:, 7)) /asize), linespec, 'color', corders(rp,:));
end
plot([0.817,0.817], [-2,2], '--k');
xlim([0.5, 1.5]);
ylim([log10(8e-2), log10(1e2) ]);

xticks([0.5,1,1.5 ]);
yticks([ -1, 0, 1, 2]);
set(gca, 'yticklabels', {' '; '10^0'; ' '; '10^2'});

xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];

xlabh=xlabel('$T$', 'interpreter', 'latex');
set(xlabh,'position',get(xlabh,'position') + [0.57, -0.3 0]);

set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (f): 0.8_1.6, Z

panelF = axes('Position', [left+width+centerW, bottom, width, height]);


kappa = 0.8;
kappa2 = 1.6;
kappastr = sprintf('%g_%g', kappa, kappa2);

n=numel(sizesD);

plot(nan, nan);hold on;
for rp=1:n
    asize = sizesD(rp);
    fname = strcat('../data/2DZBT/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    if ~exist(fname, 'file'), continue, end
    data = dlmread(fname);
    plot(1./data(:,3), log10( abs(data(:, 7)) /asize), linespec, 'color', corders(rp,:));
end

xlim([0.2, 1.2]);
ylim([log10(7e-2), log10(1e2) ]);

xticks([0.2,0.7,1.2 ]);
yticks([ -1, 0, 1, 2]);
set(gca, 'yticklabels', {' '; '10^0'; ' '; '10^2'});

xx = gca;
xx.XRuler.TickLength = [0.02, 0];
xx.YRuler.TickLength = [0.03, 0];


set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% labels

annotation('textbox',[0.445, 0.88, 0.1, 0.1], 'String', '(a)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.91, 0.88, 0.1, 0.1], 'String', '(d)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.445, 0.557, 0.1, 0.1], 'String', '(b)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.91, 0.557, 0.1, 0.1], 'String', '(e)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.445, 0.235, 0.1, 0.1], 'String', '(c)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.91, 0.235, 0.1, 0.1], 'String', '(f)', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);

annotation('textbox',[0.33, 0.85, 0.1, 0.1], 'String', sprintf("$_\\parallel\\mathrm{TM}$\n$\\kappa=0.15$\n$\\kappa'=0.3$"), 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize-2);
annotation('textbox',[0.82, 0.85, 0.1, 0.1], 'String', sprintf("$_\\parallel\\mathrm{TM}$\n$\\kappa=0.3$\n$\\kappa'=0.6$"), 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize-2);
annotation('textbox',[0.33, 0.52, 0.1, 0.1], 'String', sprintf("$_\\parallel\\mathrm{TM}$\n$\\kappa=0.3$\n$\\kappa'=1$"), 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize-2);
annotation('textbox',[0.82, 0.52, 0.1, 0.1], 'String', sprintf("$_\\parallel\\mathrm{TM}$\n$\\kappa=0.8$\n$\\kappa'=1.6$"), 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize-2);
annotation('textbox',[0.33, 0.21, 0.1, 0.1], 'String', sprintf("$_/\\mathrm{TM}$\n$\\kappa=0.8$\n$\\kappa'=0.3$"), 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize-2);
annotation('textbox',[0.82, 0.21, 0.1, 0.1], 'String', sprintf("$_/\\mathrm{TM}$\n$\\kappa=0.8$\n$\\kappa'=1.6$"), 'interpreter', 'latex', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize-2);


print('../figures/artfig_3NNclength.eps', '-depsc');