    %% Plot VRS, VK, VRS STD, and VK STD profiles for each float
    f2 = figure ('unit', 'inch', 'position',[1 1 16 9]);

    % VRS vs Depth
    nexttile;
    hold on; box on;
    ax = gca;
    ax.FontSize = 16;
    colormap("jet");
    colorValue = -5;
    for k = 1:max(dataTable{WMO_float}.Station)
        indexes =find(dataTable{WMO_float}.Station == k);
        x = datenum(dataTable{WMO_float}.dt(indexes))';
        y =dataTable{WMO_float}.Vrs_V(indexes)';
        z = dataTable{WMO_float}.Pressure_dbar(indexes)';

        plot3(x, y, z);
        %lineColor = colorValue * ones(1, length(x));
        lineColor = dataTable{WMO_float}.Temperature_degC(indexes)';
        patch([x nan], [y nan], [z nan], [lineColor nan],'EdgeColor','interp','FaceColor','none')
        %colorValue = colorValue + 5;
    end
    set(gca, 'ydir', 'reverse');
    set(gca, 'zdir', 'reverse');
    datetick('x', 'yyyymmmdd')
    grid on
    %clim([min(dataTable{WMO_float}.Temperature_degC),max(dataTable{WMO_float}.Temperature_degC)])
    %colormap(jet(max(dataTable{WMO_float}.Temperature_degC)))
    %colororder(jet)
    cb = colorbar;
    % cb.Layout.Tile ='east';
    % cb.Label.String = 'Profile';
    ylabel("VRS (Volts)",'FontSize',18);
    view(3);