clear all
close all
%script to plot the TSG data from Perseverance
B=dir('SBE3845999_2026012*.csv');
SST=[];
SSS=[];
Lon=[];
Lat=[];
for i=1:length(B)
    A=Importfile_TSG(B(i).name);
    SST=[SST;A.sst];
    SSS=[SSS;A.sss];
    Lon=[Lon;A.lon];
    Lat=[Lat;A.lat];
end

figure(1)
load coastlines
plot(coastlon,coastlat,'-','LineWidth',2)
axis([150 180 -90 -40])
hold on
scatter(Lon,Lat,20,SST)
hold off
colorbar
xlabel('Longitude')
ylabel('Latitude')
title('Temperature de la surface')

figure(2)
load coastlines
plot(coastlon,coastlat,'-','LineWidth',2)
axis([150 180 -90 -40])
hold on
I=find(SSS>30);
scatter(Lon(I),Lat(I),20,SSS(I))
hold off
colorbar
xlabel('Longitude')
ylabel('Latitude')
title('Salinite de la surface')



