mission = 7;
filename = char(gdata.Filename(mission));
load(filename, 's'); % .mat file
%%

idiel = gdata.StartProf(mission):gdata.EndProf(mission); % QC from spreadsheet
iupper = 1:5; % Defining surface layer 0 to approx 2.5 m
Temp = s.tc(:,idiel);
depth = s.depth(:,idiel);
sdn = s.sdn(idiel);
prof = 1;

x = sdn;
y = nanmean(depth(1:100,:),2);
z = Temp(1:100,:);
%% 
contourf(x,y,z)
datetick
colorbar