% Historical shoulder plot

hist_flow = xs_hist_cell{1,4};
time = 1:577;
time = time/6;

figure(1)
set(gcf,'color','white');

for i = [4 28 40 52 72]
    
    plot(time,hist_flow(:,i));
    hold on
    ylim([20 180])
    xlim([0 40])
    
end

set(gca,'fontsize',14);
xlabel('Time (hrs)');
ylabel('Discharge (m^3/s)');
legend('1 km','7 km','10 km','13 km','18 km');