% ----------------------------
% this script plots figure 9b
% ----------------------------


% load data
DBI = load('solubilityDBI_tau0.28_f1.1.dat');
DBN = load('solubilityDBN_tau0.28_f1.1.dat');

% extract packing fraction and temperature for DBI
phiDBI = DBI(:,1);
tempDBI = DBI(:,2);
[tempDBI, I] = sort(tempDBI);

% organize data
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 fraction and temperature 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);

% figure properties
xlim([0,3e-3]);
ylim([1.6,2.4])
set(gca, 'FontSize',20)
xlabel('\phi');
ylabel('T');
legend({'DBN','\tau=0.28, f=1.1'},'Location','southeast')
saveas(gcf, 'fig9b','epsc');