% letter plot for alpha_2

close all;
linewidth = 1.25;
smoothflag = true;
% set up figure
widthpix = 400;
heightpix = 300;

leftmargin = 75;
bottommargin = 72;
topmargin = 20;
rightmargin = 10;
centermargin = 0;
widthtotal = widthpix+leftmargin+rightmargin;
heighttotal = heightpix+topmargin+bottommargin;
left = leftmargin/widthtotal;
bottom = bottommargin/heighttotal;
center = centermargin/widthtotal;
width = widthpix/widthtotal;
height = heightpix/heighttotal;

psize = 26; % font size
msize = 9; % marker size
fig = figure('rend','painters','pos',[100 100 widthtotal heighttotal]);

% time evolution 8D
panelB = axes('Position', [left, bottom, width, height]);

dim = '8D';
ndim = str2double(dim(1:end-1));

phis = [12, 14.8, 17.6 19.2271, 20.8, 24];

n=length(phis);
rpp=(1:n)';
colors = get(gca, 'colororder');
markers = {'--', '--', '--', '-', '-.', '-.'};
linewidths = {1.5, 1.5, 1.5, 1.25, 1.75, 1.75};
lgds = {};
for rp=1:length(phis)
    phi = phis(rp);
    phistr = num2str(phi);
    fname = strcat('../data/rlgdynamics/', dim, '/', phistr, '/msd.dat');
    data = dlmread(fname, '', 1, 0);
    xs = data(:,1);
    ys = ndim /(ndim+2) * data(:,3) ./ (data(:,2).*data(:,2)) - 1;
    if rp==length(phis)
        xs = [xs; 1e7];
        ys = [ys; 1.33];
    end
    ys = ndim * smooth(ys);
    plot(xs*sqrt(ndim), ys, markers{rp}, 'color', colors(rp,:), 'linewidth', linewidths{rp});
    lgds{end+1} = num2str(phi/ndim, '%.3g');
    if 1==rp, hold on, end
end
lgds{end-2} = '$\hat\varphi_\mathrm{d}$';
legend(lgds, 'interpreter', 'latex', 'location', 'northwest', 'box', 'off');
set(gca, 'xscale', 'log');
xlabel('$\hat{t}$', 'interpreter', 'latex');
ylabel('$\hat\alpha_2$', 'interpreter', 'latex');
xlim([1e-1,1e7]);
ylim([-2, 25]);
xticks([1e0,1e3, 1e6]);
yticks([0, 10, 20]);
set(gca, 'fontname', 'times new roman', 'fontsize', psize);

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