% correlation length in \perp-TM and \pallel-TM

close all;

msize = 8;
psize = 20;

% set up figure
widthpix = 400;
heightpix = 250;

leftmargin = 65;
bottommargin = 50;
topmargin = 10;
rightmargin = 12;
centermargin = 0;
widthtotal = widthpix+leftmargin+rightmargin;
heighttotal = heightpix*2+topmargin+bottommargin;
left = leftmargin/widthtotal;
bottom = bottommargin/heighttotal;
width = widthpix/widthtotal;
height = heightpix/heighttotal;
center = centermargin/widthtotal;
psize = 18; % font size
msize = 9; % marker size
fig = figure('rend','painters','pos',[100 100 widthtotal heighttotal]);

kappa = 0.6;
hs = 0;
kappastr = num2str(kappa);

%% panel (a)
panelA = axes('Position', [left, bottom+height, width, height]);

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

sizesA = 8:4:24;
%sizesB = 6:4:22;

modeltype = '2DXAT';
lgds = {};

for rp=1:n
    asize = sizesA(rp);
    fname = strcat('../data/', modeltype, '/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    data = dlmread(fname);
    plot(1./data(:,2), abs(data(:, 6)./asize ), 'color', corders(rp,:));
    if 1==rp, hold on, end
end
plot([0.89,0.89], [1e-2,1e8], '--k');
set(gca, 'yscale', 'log');
yticks([1e0, 1e4, 1e8]);
xlim([0.5, 2]);
ylim([1e-1, 1e8]);
set(gca, 'xticklabel', []);
xlabel('$T$', 'interpreter', 'latex');
ylabel('$\xi_1/L$', 'interpreter', 'latex');
set(gca, 'fontname', 'times new roman', 'fontsize', psize);

%% panel (a) inset

insetA = axes('Position', [left+0.45*width, bottom+1.45*height, 0.55*width, 0.55*height]);

asize = 24;
fname = strcat('../data/',modeltype,'/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
data = dlmread(fname);
plot(1./data(:,2), abs(data(:, [7, 8, 9, 10])), '-');
hold on;
plot([0.89,0.89], [1e-2,1e3], '--k');
ylim([1e0, 5e2]);
xlim([0.5, 2]);

xticks([0.5, 1, 1.5]);
yticks([1, 1e1, 1e2]);

h=legend({'$\xi_1, \xi_2$', '$\xi_3$', '$\xi_4$', '$\xi_5$'}, 'interpreter', 'latex', 'fontsize', psize);
legend boxoff;
set(h,'Position', get(h, 'Position')+ [0.01  0.01 0 0]);
set(gca, 'yscale', 'log');
xlabel('$T$', 'interpreter', 'latex');
ylabel('$\xi$', 'interpreter', 'latex');

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

%% panel (b)

panelB = axes('Position', [left, bottom, width, height]);

sizesB = 8:2:16;
modeltype = '2DYAT';
lgds = {};
prp=1;

for rp=1:n
    asize = sizesB(rp);
    fname = strcat('../data/', modeltype, '/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
    data = dlmread(fname);
    plot(1./data(:,2), abs(data(:, 6)./asize ), 'color', corders(rp,:));
    if 1==rp, hold on, end
end
plot([0.89,0.89], [1e-1,1e3], '--k');
set(gca, 'yscale', 'log');
xlim([0.5, 2]);
ylim([3e-1, 2e1]);
xlabel('$T$', 'interpreter', 'latex');
ylabel('$\xi_1/L$', 'interpreter', 'latex');

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

%% panel (b) inset

insetB = axes('Position', [left+0.45*width, bottom+0.45*height, 0.55*width, 0.55*height]);

asize = 16;
fname = strcat('../data/',modeltype,'/clength_', kappastr, '_', num2str(hs) , '_', num2str(asize), '.dat');
data = dlmread(fname);
plot(1./data(:,2), abs(data(:, [7, 8, 9])), '-');
hold on;
plot([0.89,0.89], [1e-2,1e3], '--k');
ylim([3e-1, 5e2]);
xlim([0.5, 2]);

xticks([0.5, 1, 1.5]);
yticks([1e0, 1e1, 1e2]);

legend({'$\xi_1, \xi_2$', '$\xi_3$', '$\xi_4$'}, 'interpreter', 'latex', 'fontsize', psize);
legend boxoff;
set(gca, 'yscale', 'log');
xlabel('$T$', 'interpreter', 'latex');
ylabel('$\xi$', 'interpreter', 'latex');

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

%% l=panel labels
annotation('textbox',[0.14, 0.5, 0.1, 0.1], 'String', '(a)', 'EdgeColor','none','fontname', 'times new roman','fontsize',psize+2);
annotation('textbox',[0.14, 0.05, 0.1, 0.1], 'String', '(b)', 'EdgeColor','none','fontname', 'times new roman','fontsize',psize+2);

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