load("24403401_QC.mat");
s.profile = 1:length(s.lat);

%

station = [s.profile; s.profile];
depths = [zeros(size(s.profile)); max(s.pres)];

figure("Name","Depth vs Profile")
plot(station,depths)
set(gca,'Ydir','reverse')
title('Mission Depth Profile')

%

figure("Name","vrs std")
scatter(max(s.pres),max(s.vrs_std,[],1,"omitmissing"),'k.')
title('max vrs std vs max pressure')

%

figure("Name","vk std")
scatter(max(s.pres),max(s.vk_std,[],1,"omitmissing"),'k.')
title('max vk std vs max pressure')

%

figure("Name","Ik")
scatter(max(s.pres),max(s.Ik,[],1,"omitmissing"),'k.')
title('max Ik vs max pressure')

%

figure("Name","Ib")
scatter(max(s.pres),max(s.Ib,[],1,"omitmissing"),'k.')
title('max Ib vs max pressure')

%

figure("Name","Ik and Ib")
scatter(max(s.pres),max(s.Ik,[],1,"omitmissing"),5,max(s.Ib,[],1,"omitmissing"),'filled')
title('max Ik vs max pressure, colored by max Ib')
colorbar