%-------------------------------
% this script plots figure 6
%-------------------------------

% f=0.5, nw=10-------------------------------------------------------

% load data
DBI = load('solubilityDBI_nw10_f0.5.dat');
DBN = load('solubilityDBN_nw10_f0.5.dat');

% extract DBI solubility data
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);

% sort data
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);

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

% extract DBN data
phiDBN = DBN(:,1);
tempDBN = DBN(:,2);

% plot DBN solubility line
plot(phiDBN, tempDBN,'k','lineWidth',2);
hold on

% plot DBI solubility line
plot(Vq, X,'color',[248/255,181/255,113/255],'lineWidth',2);

% f=0.5, nw=30--------------------------------------------------------

% load data
DBI = load('solubilityDBI_nw30_f0.5.dat');

% extract DBI solubility data
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);

% sort data
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);

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

% plot 
plot(Vq, X,'color',[106/255,175/255,226/255],'lineWidth',2);

% f=0.5, nw=50--------------------------------------------------------

% load data
DBI = load('solubilityDBI_nw35_f0.5.dat');

% extract DBI solubility data
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);

% sort data
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);

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

% plot
plot(Vq, X,'color',[42/255,125/255,194/255],'lineWidth',2);

% f=0.5, nw=50---------------------------------------------------------

% load data
DBI = load('solubilityDBI_nw50_f0.5.dat');

% extract DBI solubility data
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);

% sort data
[tempDBI, I] = sort(tempDBI);
phiDBI = phiDBI(I);

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


% plot
plot(Vq, X,'color',[45/255,67/255,151/255],'lineWidth',2);


% plot properties -----------------------------------------------------
set(gca, 'FontSize',20)
xlim([0,2e-3]);
ylim([0.7,1.2])
xlabel('\phi');
ylabel('T');


legend('DBN','n_w=10','n_w=30','n_w=35','n_w=50','Location','southeast')
saveas(gcf, 'fig6','epsc');