close all;
clear all;


load NPEO2008_CTD1_noHeader.txt -ASCII

depth1=NPEO2008_CTD1_noHeader(:,1);
temp1=NPEO2008_CTD1_noHeader(:,3);
salinity1=NPEO2008_CTD1_noHeader(:,5);


load NPEO2008_CTD2_noHeader.txt -ASCII

depth2=NPEO2008_CTD2_noHeader(:,1);
temp2=NPEO2008_CTD2_noHeader(:,3);
salinity2=NPEO2008_CTD2_noHeader(:,6);


figure(1)
plot(temp1,depth1,'red')
hold on;
plot(temp2,depth2,'blue')
set(gca, 'YDir','reverse')
title('Temperature profile')
xlabel('Temperature (C)') % x-axis label
ylabel('Depth (m)') % y-axis label

figure(2)
plot(salinity1,depth1,'red')
hold on;
plot(salinity2,depth2,'blue')
set(gca, 'YDir','reverse')
title('Salinity profile')
xlabel('Salinity (psu)') % x-axis label
ylabel('Depth (m)') % y-axis label

figure(3)
plot(salinity1,temp1,'red')
hold on;
plot(salinity2,temp2,'blue')
title('T-S diagram')
xlabel('Salinity (psu)') % x-axis label
ylabel('Temperature (C)') % y-axis label
