%plotPT

load p_salinity_Kehua_Liu_NaCl.mat
T_salinity=(0:0.2:40);
salinity=(0:0.5:40)/100;
T=(0:0.2:40);
N=length(T);
M=length(S);

for j=1:M

for i=1:N
   P7(i)=interpolation_pressure(p_salinity,T_salinity,salinity,S(j)/10,T(i));
end


T(31)
S_press(j)=P7(31);

end





% figure(1)
% plot(T,P7,'k.')
% ylabel('Pressure (MPa)','FontSize',15)
% xlabel('Temperature (\circC)','FontSize',15)
% ylim([0 16])
% xlim([0 10])
% set(gca,'XTick',(0:1:10),'YTick',(0:2:16))