function [] = Load_Vel_Plot(d1,d2,d3,HorizontalLoad,Slip,Test_Name,Time,VerticalLoad)

red_fac = 10;

%Coefficient of friction [-]
mu = (VerticalLoad)./(HorizontalLoad);
%Time step between each recorded piece of data [sec]
Stamp = diff(Time);
%Velocity
vel=diff(Slip)./Stamp;
vel=abs(vel);

%Plot whole experiment
d2 = d3;

ind_use = d1:d2;
pp = find(vel(d1:d2)<0.1);
ind_use(rem(pp,red_fac) > 0) = [];

cc = [Slip(ind_use)' Slip(ind_use)'] - Slip(d1);
xx = [vel(ind_use)' vel(ind_use)'];
yy = [mu(ind_use)' mu(ind_use)'];
zz = zeros(size(xx));
cc(cc>35) = 34.9999;

colorsbrewer = flip(brewermap(100,'RdYlGn'));

fh = figure;
surf(flip(xx),flip(yy),zz,flip(cc),'EdgeColor','interp','MarkerSize',0.01) ;
ylab = ylabel('$$\frac{F_v}{F_n}$$ [-]');
xlab = xlabel('Velocity [$$mm \cdot sec^{-1}$$]');
set(xlab,'Interpreter','latex','fontsize',12)
set(ylab,'Interpreter','latex','fontsize',12)
set(gca, 'XScale', 'log')
set(gca,'fontsize',11) 
set(fh, 'color', 'white')
set(gca, 'Box', 'on' );
xlim([1e-4 200])
ylim([0.4 0.9])
colormap(colorsbrewer) ;
view(2) %// view(0,90)          %// set view in X-Y plane
c = colorbar;
caxis([0 35])
c.Label.String = 'Slip [mm]';
c.Label.Interpreter = 'latex';
% 
% saveas(gcf,[Test_Name,'_FvFn_Vel.png']);

end

