% percofig dynamics collapse
close all;
pccl = {1.128, 3.510, 6.247, 9.170, 12.22, 15.4, 18.6};
muminus = [3, 2.528, 1.3377, 0.73, 0.31];
%mus = [1, 1.31, 2.88, 4.37,5.717]; % mu
mus = [1, 1.31, 2.88, 4.37,5.717];
dwprimes = [2, 3.036, 6.25, 14, 39];
dimstrs = {'4D', '5D'};

pickA =  {{'5', '5.8', '6', '6.15', '6.4', '6.5', '6.8', '7.2'}, ... '6.15', 'phip', '8'
    {'8', '8.4', '8.6', '8.9', '9.4', '9.8', '10.3', '12'}, ...
    };
% set up figure
% labels = {'(a)', '(b)', '(c)', '(d)'};

%% figure setup
widthpix = 375;
heightpix = 300;

leftmargin = 85;
bottommargin = 80;
topmargin = 20;
rightmargin = 20;
centermargin = 0;
widthtotal = 2*widthpix+leftmargin+rightmargin;
heighttotal = heightpix+topmargin+bottommargin;
left = leftmargin/widthtotal;
bottom = bottommargin/heighttotal;
width = widthpix/widthtotal;
height = heightpix/heighttotal;
center = centermargin/widthtotal;
psize = 28; % font size
msize = 9; % marker size

fig = figure('rend','painters','pos',[100 100 widthtotal heighttotal]);
% set color order
n=8;
rpp=(1:n)';
CollorOrders = [(rpp-1)/(n-1), zeros(n,1), 1 - (rpp-1)/(n-1)];
colororder(CollorOrders);
%% panels

for rp=1:length(dimstrs)
 
    dimstr = dimstrs{rp};
    ndim = str2double(dimstr(1:end-1));
    picks = pickA{rp};
    nfnames = length(picks);
    pc = pccl{ndim-1};
    
    if 1==rp
        nowpanel = axes('Position', [left, bottom, width, height]);
        intc = 0.75;
    elseif 2==rp
        nowpanel = axes('Position', [left+width, bottom, width, height]);
        intc = 0.9;
    end
    
    for rq=1:nfnames
        strpf = picks{rq};
        if strcmp(strpf, 'phip')
            dbpf = pc;
            strpf = num2str(dbpf);
        else
            dbpf = str2double(strpf);
        end
        fname = strcat('../data/rlgdynamics/', dimstr, '/msd_', strpf, '.dat');

        data=dlmread(fname,'');
        epsilons = (dbpf-pc)/pc;
        xs = abs(epsilons)^(mus(ndim)+muminus(ndim))*data(:,1);
        ys = abs(epsilons)^(muminus(ndim))*data(:,2);
        plot(xs, ys, 'linewidth', 1);
        if 1==rq, hold on, end
    end
    % subdiffusion line
    xs = linspace(log(1e-8), log(1e8));
    ys = 2/dwprimes(ndim)*xs+intc;
    plot(exp(xs), exp(ys), '--k', 'linewidth', 1.25);
    xlim([1e-8, 1e4]);
    ylim([1e-4, 1e4]);
    set(gca, 'xscale', 'log');
    set(gca, 'yscale', 'log');
    if 1==rp
        %set(nowpanel, 'xticklabels', []);
        xticks([1e-8, 1e-4, 1e0]);
        yticks([1e-4, 1e0, 1e4]);
        ylabh = ylabel('$\widetilde{\Delta}$', 'interpreter', 'latex');
        set(ylabh, 'Units', 'Normalized', 'Position', [-0.11, 0.5, 0]);
        xlabel('$\tilde{t}$', 'interpreter', 'latex');
    else
        %yticks([1e-4, 1e0]);
        xticks([1e-8, 1e-4, 1e0, 1e4]);
        set(nowpanel, 'yticklabels', []);
        xlabel('$\tilde{t}$', 'interpreter', 'latex');
    end
    legend(picks, 'location', 'southeast', 'NumColumns', 2);
    legend boxoff;
    set(nowpanel,  'xscale', 'log', 'yscale', 'log');
    set(nowpanel, 'fontname', 'times new roman', 'fontsize', psize);
end

annotation('textbox',[0.1, 0.86, 0.1, 0.1], 'String', '(a) 4D', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);
annotation('textbox',[0.54, 0.86, 0.1, 0.1], 'String', '(b) 5D', 'EdgeColor','none', 'Fontname', 'Times New Roman','fontsize',psize);

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