function [] = func_plot_soil_res(lon2d, lat2d, hgt,lh,sh, season)
% Summary of this function goes here
%   
make_it_tight = true;
subplot = @(m,n,p) subtightplot (m, n, p, [0.045 0.05], [0.12 0.025], [0.05 0.05]);
if ~make_it_tight,  clear subplot;  end

ann_list = char(97:122);
season_list = {'DJF', 'MAM', 'JJA', 'SON', 'ANN'};

iplot = 1;

% lon2d = lon2d(1:200,100:192);
% lat2d = lat2d(1:200,100:192);
% lh    = lh(1:200,100:192,:,:);
% sh    = sh(1:200,100:192,:,:);
% hgt   = hgt(1:200,100:192);

figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 1000 600])

for ics = 9:11
%-----------------------------------------------------
subplot(2,3,iplot)
surf(lon2d, lat2d, hgt,lh(:,:,ics,season)); 
shading interp; view(20,50); 
if iplot ==1; zlabel('LH'); end
caxis([-8 8])
colormap (gca, bluewhitered  )
% h = colorbar;
zlim([0 800]);
if iplot ==1; title('routing\_100m - routing\_250m'); end
if iplot ==2; title('routing\_250m - routing\_500m'); end
if iplot ==3; title('routing\_500m - routing\_1000m'); end

ann = [ann_list(iplot), ')'];
text(-100, 35, 1200, ann);

%-----------------------------------------------------
subplot(2,3,iplot+3)
surf(lon2d, lat2d, hgt,sh(:,:,ics,season)); 
shading interp; view(20,50); 
% xlabel('Longitude')
if iplot ==1; zlabel('SH'); end
caxis([-8 8])
colormap (gca, bluewhitered  )
% h = colorbar;
zlim([0 800]);
ann = [ann_list(iplot+3), ')'];
text(-100, 35, 1200, ann);
%-----------------------------------------------------

iplot = iplot + 1;


end

hp3 = get(subplot(2,3,4),'Position')
h = colorbar('southoutside')
set(h, 'Position', [hp3(1) + 0.27  hp3(2)-0.07  hp3(3)*1.45 0.01 ]);
h.Label.String = ['SH or LH diff. (W m^{-2})'];

print('-depsc', ['1km_LHSH_res_cumulative']);

%% 
iplot = 1;

figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 1000 600])

for ics = 6:8
%-----------------------------------------------------
subplot(2,3,iplot)
surf(lon2d, lat2d, hgt,lh(:,:,ics,season)); 
shading interp; view(20,50); 
if iplot ==1; zlabel('LH'); end
caxis([-8 8])
colormap (gca, bluewhitered  )
zlim([0 800]);
if iplot ==1; title('routing\_100m - routing\_250m'); end
if iplot ==2; title('routing\_100m - routing\_500m'); end
if iplot ==3; title('routing\_100m - routing\_1000m'); end
ann = [ann_list(iplot), ')'];
text(-100, 35, 1200, ann);

%-----------------------------------------------------
subplot(2,3,iplot+3)
surf(lon2d, lat2d, hgt,sh(:,:,ics,season)); 
shading interp; view(20,50); 
if iplot ==1; zlabel('SH'); end
caxis([-8 8])
colormap (gca, bluewhitered  )
% h = colorbar;
zlim([0 800]);
% title('SM2')
ann = [ann_list(iplot+3), ')'];
text(-100, 35, 1200, ann);
%-----------------------------------------------------

iplot = iplot + 1;

end

hp3 = get(subplot(2,3,4),'Position')
h = colorbar('southoutside')
set(h, 'Position', [hp3(1) + 0.27  hp3(2)-0.07  hp3(3)*1.45 0.01 ]);
h.Label.String = ['SH or LH diff. (W m^{-2})'];

print('-depsc', ['1km_LHSH_res']);


end

