close all

yLimits = [10*floor(m/10) 10];

subplot(2,2,1), plot(profileDateTimeNum, profileDepth);
set(gcf, 'units', 'normalized', 'Position', [.1, .1, .5, .5]);
set(gca, 'Ylim', yLimits);
grid;
datetick('x');
title('Profile Depth vs. Time');

subplot(2,2,2);
line(cpT, -cpP, 'Color', 'k')
set(gca, 'YLim', yLimits);
ax1 = gca; % current axes
ax1.XColor = 'k';
ax1.YColor = 'k';

ax1_pos = ax1.Position; % position of first axes
ax2 = axes('Position', ax1_pos,...
    'XAxisLocation', 'top', ...
    'YAxisLocation', 'right', ...
    'XColor', 'r', 'YColor', 'r', ...
    'Color', 'none');
set(gca, 'Ylim', yLimits);
line(cpS, -cpP, 'Parent', ax2, 'Color', 'r')

grid;
set(ax1.XLabel, 'String', 'CP Temperature');
set(ax2.XLabel, 'String', 'CP Salinity');

subplot(2,2,3), plot(optodeO2Conc(2:end), optodeDepth(2: end));
set(gca, 'Ylim', yLimits);
grid;
title('O2 Concentration vs Depth');

subplot(2,2,4)
line(flbbFL(3:end), flbbDepth(3:end), 'Color', 'g')
set(gca, 'Ylim', yLimits);
ax1 = gca; % current axes
ax1.XColor = 'g';
ax1.YColor = 'g';

ax1_pos = ax1.Position; % position of first axes
ax2 = axes('Position', ax1_pos,...
    'XAxisLocation', 'top', ...
    'YAxisLocation', 'right', ...
    'XColor', 'r', 'YColor', 'r', ...
    'Color', 'none');
set(gca, 'Ylim', yLimits);
line(flbbBB(3:end), flbbDepth(3:end), 'Parent', ax2, 'Color', 'r')

grid;
set(ax1.XLabel, 'String', 'FLBB FL');
set(ax2.XLabel, 'String', 'FLBB BB');



