% -----------------------------
% This script plots figure 7
% -----------------------------


% sw=0 -----------------------------------

% load data
DBI = load('solubilityDBI_sw0.dat');
DBN = load('solubilityDBN_sw0.dat');

% extract DBI solubility info
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 solubility info
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);


% sw=-10 -----------------------------------

% load data
DBI = load('solubilityDBI_sw-10.dat');
DBN = load('solubilityDBN_sw-10.dat');

% extract DBI solubility info
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 data
plot(Vq, X,'color',[0.55,0.8,0.92], 'lineWidth',2);

% sw=-15 -----------------------------------

% load data
DBI = load('solubilityDBI_sw-15.dat');
DBN = load('solubilityDBN_sw-15.dat');

% extract DBI solubility info
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 data
plot(Vq, X,'color',[0.79,0.4,0.47], 'lineWidth',2);
hold on;

% sw=-30 -----------------------------------

% load data
DBI = load('solubilityDBI_sw-30.dat');
DBN = load('solubilityDBN_sw-30.dat');

% extract DBI solubility info
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 data
plot(Vq, X,'color',[0.11,0.47,0.24], 'lineWidth',2);
hold on;

% sw=-50 -----------------------------------

% load data
DBI = load('solubilityDBI_sw-50.dat');
DBN = load('solubilityDBN_sw-50.dat');

% extract DBI solubility info
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 data
plot(Vq, X,'color',[221/255,204/255,119/255], 'lineWidth',2);
hold on;

% figure properties ----------------------------------
xlim([0,2.e-3])
ylim([1.5,2.3])
xlabel('\phi');
ylabel('T');
set(gca, 'FontSize',20)
legend({'DBN','\Delta s_w=0','\Delta s_w=-10','\Delta s_w=-15','\Delta s_w=-30','\Delta s_w=-50'},'location','southeast','Fontsize',14)
saveas(gcf, 'fig7','epsc');