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

season = 5;

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

% sdf1 = sdf1(1:200,100:192,:,:);
% sdf2 = sdf2(1:200,100:192,:,:);
% slope= slope(1:200,100:192);

figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 800 200])
x = [-0.25:0.0001:0.25];

%-----------------------------------------------------
subplot(1,2,1)
    
    sdata = sdf1(:,:,6,season);
    [xx, grd_1] = func_calc_slope(slope, sdata); hold on; 
    plot(xx, grd_1, 'r','LineWidth', 1.5);
    sdata = sdf1(:,:,7,season);
    [xx, grd_1] = func_calc_slope(slope, sdata); hold on; 
    plot(xx, grd_1, 'g','LineWidth', 1.5); 
    sdata = sdf1(:,:,8,season);
    [xx, grd_1] = func_calc_slope(slope, sdata); hold on; 
    plot(xx, grd_1, 'b','LineWidth', 1.5); 
    grid on; 
    plot(0:0.01:1.2, 0*(0:0.01:1.2), 'k.');
    ylim([-8 8])
    legend('routing 100m - routing 250m','routing 100m - routing 500m', 'routing 100m - routing 1000m');
    ylabel('d(SM) / d(topo. slope)')

    title_left('a) LH')

%-----------------------------------------------------
subplot(1,2,2)
    
    sdata = sdf2(:,:,6,season);
    [xx, grd_1] = func_calc_slope(slope, sdata); hold on; 
    plot(xx, grd_1, 'r','LineWidth', 1.5);
    sdata = sdf2(:,:,7,season);
    [xx, grd_1] = func_calc_slope(slope, sdata); hold on; 
    plot(xx, grd_1, 'g','LineWidth', 1.5); 
    sdata = sdf2(:,:,8,season);
    [xx, grd_1] = func_calc_slope(slope, sdata); hold on; 
    plot(xx, grd_1, 'b','LineWidth', 1.5); 
    grid on; 
    plot(0:0.01:1.2, 0*(0:0.01:1.2), 'k.');
    ylim([-8 8])
    title_left('b) SH')
%-----------------------------------------------------


end