function f = plotPressureWaterfall(cycle, numCycles)
figure('name', "Waterfall", 'NumberTitle', 'off');
hold on;

colormap(jet);
cmap = colormap;
 
for i=1:numCycles
 %   x = datenum(cycle(i).dt)';
    x = datetime(cycle(i).dt)';
    y = cycle(i).pressureDBar';
    z = cycle(i).vrsVolts';

    c = (mean(cycle(i).thermTemp) * 10 * (25.6/25.0));
    c = round(c);

    %lineColor = colorValue * ones(1, length(x));
    plot3(x, y, z, 'Color', cmap(c,:));
    % h = patch(x, y, z, c);
    % h.FaceColor='none';
    % h.EdgeColor='interp';
    % patch([x nan], [y nan], [z nan], [c nan],'EdgeColor','interp','FaceColor','none')
    % if(mod(i,2) == 0)
    %     colorValue = colorValue + 5;
    % end
end
set(gca, 'ydir', 'reverse');
zLimits = zlim;
set(gca, 'ztick', [zLimits(1):0.001:zLimits(2)+.001]);
view(3);
datetick('x', 'yyyymmmdd HH:MM')
grid on
set(gca, 'CLim', [0, 25.6]); 
colorbar;
xStart = datenum(cycle(1).dt(1));
xEnd = datenum(cycle(length(cycle)).dt(end));
xlim([xStart xEnd]);
end