maxdepth = max(abs(s.depth),[],'omitmissing');
time = datetime(datevec(s.sdn),'TimeZone','UTC');
plot(time,maxdepth)
set(gca,'Ydir','reverse')