clear all:
close all;

%nRuns=100;
numOfBoxes=5;
nParticles=16*(numOfBoxes^3)
s=0.02;

if s==0.01;
    sStr='0.01';
    phiStr={'0.710','0.720','0.730','0.740','0.745'};
end
if s==0.02;
    sStr='0.02';
    phiStr={'0.680','0.690','0.700','0.710','0.720','0.725','0.730'};
end
if s==0.03;
    sStr='0.03';
    phiStr={'0.680','0.690','0.700','0.710','0.715'};
end

%cM={'sr','xb','og'};
%cM={'sr','sb','sg','sm','sc','sk'};
cL={'-r','-b','-g','-m','-c','-k','-r','-b','-g','-m','-c'};
%cLM={'s-r','s-b','s-g','s-m','s-c','s-k'};

cM={'sr','sb','sg','sm','sc','sk','or','ob','og','om','oc','ok'};
cLM={'s-r','s-b','s-g','s-m','s-c','s-k','o-r','o-b','o-g','o-m','o-c','o-k'};

dim1=500;
dim2=400;
fig1=figure('position',[0,0,dim1,dim2]);

topDir1 = ['./gardnerHS1/deltaHS1/'];
topDir2 = ['./gardnerHS1/deltaHS1AB/'];
systemSize = ['N' num2str(nParticles) '_s' sStr ];
midDir1 = [ topDir1 systemSize filesep ];
midDir2 = [ topDir2 systemSize filesep ];

for i =1:numel(phiStr)
    i
    phiTemp = phiStr{i};
    dataFile1 = [ midDir1 'delta' phiTemp '.dat' ];
    dataFile2 = [ midDir2 'dAB' phiTemp '.dat' ];

    D1 = load(dataFile1);
    D2 = load(dataFile2);

    loglog( D1(:,1), D1(:,2), strjoin(cM(i)),'DisplayName',['d, \phi =' strjoin(phiStr(i))]);
    hold on;
    D2 = load(dataFile2);
    loglog( D2(:,1), D2(:,2), strjoin(cLM(i)),'DisplayName',['dAB, \phi =' strjoin(phiStr(i))]);
end

%axis equal
%xlim([1e+0,2e+7])
%ylim([5e-5,2e-2])

title([ 'HS1   MSD   nParticles = ' num2str(nParticles) '   s = ' num2str(s)]);
ylabel('$$\Delta/\bar{\sigma}^2$$','Interpreter','Latex')
xlabel('time (MC steps)');
legend('-DynamicLegend','Location', 'Best');