function [] = plot_smtbl_sct_diff(lat2d, lon2d, bb_diff, drysmc_diff, f11_diff, maxsmc_diff,...
    refsmc_diff, satpsi_diff, satdk_diff, satdw_diff, wltsmc_diff, qtz_diff)

close all
make_it_tight = true;
subplot = @(m,n,p) subtightplot (m, n, p, [0.04 0.01], [0.06 0.025], [0.03 0.03]);
if ~make_it_tight,  clear subplot;  end

set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [0 0 1000 600 ])
% figure default 
figure(1)
%---------------------------------------------------------------
subplot(3,3,1)
surf(lon2d, lat2d,bb_diff); view(2); shading interp;
caxis([ -1 1]);colormap (bluewhitered(20)); 
h=colorbar; title(h, '[-]')
title_left('a) Beta ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'xticklabel',{[]})
box on;
%---------------------------------------------------------------
subplot(3,3,2)
surf(lon2d, lat2d,drysmc_diff); view(2); shading interp;
max_lim = max(abs(drysmc_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m^3/m^{3}]')
title_left('b) Dry soil moisture thershold  ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'xticklabel',{[]})
set(gca,'yticklabel',{[]})
box on;

%---------------------------------------------------------------
subplot(3,3,3)
surf(lon2d, lat2d,maxsmc_diff); view(2); shading interp;
max_lim = max(abs(maxsmc_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m^3/m^{3}]')
title_left('c) Porosity  ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'xticklabel',{[]})
set(gca,'yticklabel',{[]})
box on;

%---------------------------------------------------------------
subplot(3,3,4)
surf(lon2d, lat2d,refsmc_diff); view(2); shading interp;
max_lim = max(abs(refsmc_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m^3/m^{3}]')
title_left('d) Field capacity ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'xticklabel',{[]})
box on;

%---------------------------------------------------------------
subplot(3,3,5)
surf(lon2d, lat2d,satpsi_diff); view(2); shading interp;
max_lim = max(abs(satpsi_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m]' )
title_left('e) Soil potential')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'xticklabel',{[]})
set(gca,'yticklabel',{[]})
box on;

%---------------------------------------------------------------
subplot(3,3,6)
surf(lon2d, lat2d,satdk_diff); view(2); shading interp;
max_lim = max(abs(satdk_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m/s]')
title_left('f) Saturated hydraulic conductivity  ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'yticklabel',{[]})
box on;

%---------------------------------------------------------------
subplot(3,3,7)
surf(lon2d, lat2d,satdw_diff); view(2); shading interp;
max_lim = max(abs(satdw_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m^2/s]')
title_left('g) Saturated hydraulic diffusivity  ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
box on;

%---------------------------------------------------------------
subplot(3,3,8)
surf(lon2d, lat2d,wltsmc_diff); view(2); shading interp;
max_lim = max(abs(wltsmc_diff(:)));
caxis([-1*max_lim max_lim]);
colormap (bluewhitered);
h=colorbar; title(h, '[m^3/m^{3}]')
title_left('h) Wilting point soil moisture  ')
ylim([35.1, 38.1])
xlim([-99.3,-95.6])
set(gca,'yticklabel',{[]})
box on;

%---------------------------------------------------------------
end