function [] = func_plot_soil_slope_bar_p1(sdf1, sdf2, slope)

season = 5;

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

figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 300 300])
x = [-0.25:0.0001:0.25];
ymax = 3;
%-----------------------------------------------------
subplot(2,1,1)
    
    sdata = sdf1(:,:,2,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-sr','LineWidth', 1.5);
    sdata = sdf1(:,:,3,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-gx','LineWidth', 1.5); 
    sdata = sdf1(:,:,4,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-bo','LineWidth', 1.5); 
    grid on; 
    plot(0:0.01:0.6, 0*(0:0.01:0.6), 'k.');
    ylim([-1*ymax ymax])
    legend('routing 100m - routing 250m','routing 100m - routing 500m', 'routing 100m - routing 1000m', 'Location', 'Southeast');
    ylabel('LH diff (W m^{-2})')
    text(0.05,ymax*-0.85, 'a) LH')
    title('LSM 4km')
%-----------------------------------------------------
subplot(2,1,2)
    
    sdata = sdf2(:,:,2,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-sr','LineWidth', 1.5);
    sdata = sdf2(:,:,3,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-gx','LineWidth', 1.5); 
    sdata = sdf2(:,:,4,season);
    [xx, grd_1, err] = func_calc_slope_barplot(slope, sdata); hold on; 
    plot(xx, grd_1, '-bo','LineWidth', 1.5); 
    grid on; 
    plot(0:0.01:0.6, 0*(0:0.01:0.6), 'k.');
    ylim([-1*ymax ymax])
    ylabel('SH diff (W m^{-2})')
    text(0.05,ymax*-0.85,'b) SH')
%-----------------------------------------------------

end