function [] = func_plot_soil_slope_bar_p3(sdf1, sdf2, sdf3, sdf4, slope)

season = 5;

make_it_tight = true;
subplot = @(m,n,p) subtightplot (m, n, p, [0.05 0.055], [0.1 0.045], [0.05 0.05]);
if ~make_it_tight,  clear subplot;  end

figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 300 600])
x = [-0.25:0.0001:0.25];
ymax = 3*10^-3;
%-----------------------------------------------------
subplot(4,1,1)
    
    sdata = sdf1(:,:,2,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-xr','LineWidth', 1.5);
    sdata = sdf1(:,:,3,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-^b','LineWidth', 1.5); 
    grid on; 
    
    plot(0:0.01:0.35, 0*(0:0.01:0.35), 'k.');
    ylim([-1*ymax ymax])
    legend('routing 250m - routing 500m', 'routing 250m - routing 1000m', 'Location', 'Southeast');
%     ylabel('mean SM diff (m^3 m^{-3})')
    
    text(0.0125, ymax*-0.85, 'i) SM1')
    title('LSM 10km')
%-----------------------------------------------------
subplot(4,1,2)
    
    sdata = sdf2(:,:,2,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-xr','LineWidth', 1.5);
    sdata = sdf2(:,:,3,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-^b','LineWidth', 1.5); 

    grid on; 
    plot(0:0.01:0.35, 0*(0:0.01:0.35), 'k.');
    ylim([-1*ymax ymax])
    text(0.0125, ymax*-0.85,'j) SM2')
%-----------------------------------------------------
subplot(4,1,3)
    sdata = sdf3(:,:,2,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-xr','LineWidth', 1.5);
    sdata = sdf3(:,:,3,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-^b','LineWidth', 1.5); 

    grid on; 
    plot(0:0.01:0.35, 0*(0:0.01:0.35), 'k.');
    ylim([-1*ymax ymax])
%     ylabel('mean SM diff (m^3 m^{-3})')
%     xlabel('topo. slope')
    text(0.0125, ymax*-0.85,'k) SM3')

%-----------------------------------------------------
subplot(4,1,4)

    sdata = sdf4(:,:,2,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-xr','LineWidth', 1.5);
    sdata = sdf4(:,:,3,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-^b','LineWidth', 1.5); 

    grid on; 
    plot(0:0.01:0.35, 0*(0:0.01:0.35), 'k.');
    ylim([-1*ymax ymax])
%     ylabel('d(SM)/d(topo. slope)')
    xlabel('topo. slope')
    text(0.0125, ymax*-0.85,'l) SM4')


end