% --------------------------------
% this script plots figure 9a
% --------------------------------


% tau=0.01 ---------------------------------------------------------------

% load data
DBI = load('solubilityDBI_tau0.01.dat');
DBN = load('solubilityDBN_tau0.01.dat');

% extract packing fractions and temp for DBI
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);
ind = find(tempDBI > 1.7);
phiDBI = phiDBI(ind);
tempDBI = tempDBI(ind);

% spline interpolation for smoothing
X = (min(tempDBI):0.001:max(tempDBI));
Vq = interp1(tempDBI,phiDBI,X,'spline');

% extract packing fractions and temp for DBN
phiDBN = DBN(:,1);
tempDBN = DBN(:,2);

% plot data
plot(phiDBN, tempDBN,'color',[.0,0.,0.],'lineWidth',2);
hold on;
plot(Vq, X,'color',[51/255,34/255,136/255],'lineWidth',2);

% tau=0.05 ---------------------------------------------------------------

% load data
DBI = load('solubilityDBI_c4_f1.dat');
DBN = load('solubilityDBN_c4_f1.dat');

% extract packing fractions and temp for DBI
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);
ind = find(tempDBI > 1.7);
phiDBI = phiDBI(ind);
tempDBI = tempDBI(ind);

% spline interpolation for smoothing
X = (min(tempDBI):0.001:max(tempDBI));
Vq = interp1(tempDBI,phiDBI,X,'spline');

% extract packing fractions and temp for DBN
phiDBN = DBN(:,1);
tempDBN = DBN(:,2);

% plot data
plot(Vq, X,'color',[68/255,170/255,153/255],'lineWidth',2);
hold on;

% tau=0.1 ----------------------------------------------------------------

% load data
DBI = load('solubilityDBI_tau0.25.dat');
DBN = load('solubilityDBN_tau0.25.dat');

% extract packing fractions and temp for DBI
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);
ind = find(tempDBI > 1.7);
phiDBI = phiDBI(ind);
tempDBI = tempDBI(ind);

% spline interpolation for smoothing
X = (min(tempDBI):0.001:max(tempDBI));
Vq = interp1(tempDBI,phiDBI,X,'spline');

% extract packing fractions and temp for DBN
phiDBN = DBN(:,1);
tempDBN = DBN(:,2);

% plot data
plot(Vq, X,'color',[204/255,102/255,119/255],'lineWidth',2);
hold on;

% tau=0.25 ---------------------------------------------------------------

% load data
DBI = load('solubilityDBI_tau0.35.dat');
DBN = load('solubilityDBN_tau0.35.dat');

% extract packing fractions and temp for DBI
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);
ind = find(tempDBI > 1.7);
phiDBI = phiDBI(ind);
tempDBI = tempDBI(ind);

% spline interpolation for smoothing
X = (min(tempDBI):0.001:max(tempDBI));
Vq = interp1(tempDBI,phiDBI,X,'spline');

% extract packing fractions and temp for DBN
phiDBN = DBN(:,1);
tempDBN = DBN(:,2);

% plot data
plot(Vq, X,'color',[136/255,34/255,85/255],'lineWidth',2);
hold on;
xlim([0,3e-3])
ylim([1.4,2.3])


% figure properties ----------------------
set(gca, 'FontSize',20)
xlabel('\phi');
ylabel('T');
legend({'DBN','\tau=0.01','\tau=0.05','\tau=0.25','\tau=0.35'},'Location','southeast')
saveas(gcf, 'fig9a','epsc');