function [] = func_plot_soil_vert_profile(lon2d, lat2d, hgt,soil_1,soil_2,soil_3,soil_4, season)
% Summary of this function goes here
%   
make_it_tight = true;
subplot = @(m,n,p) subtightplot (m, n, p, [0.055 0.075], [0.12 0.025], [0.05 0.05]);
if ~make_it_tight,  clear subplot;  end

season_list = {'DJF', 'MAM', 'JJA', 'SON'};
case_list   = {'both routing - no routing',...
               'surface routing - no routing',...
               'subsurface routing - no routing', ...
               'both routing - subsurface routing', ...
               'both routing - surface routing'...
               };

for ics = 1:5
figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 1000 600])
%-----------------------------------------------------
subplot(2,2,1)
surf(lon2d, lat2d, hgt,soil_1(:,:,ics,season)); 
shading interp; view(20,50); 
xlabel('Longitude')
ylabel('Latitude')
zlabel('Surface Elevation (m)')
caxis([-0.025 0.025])
colormap (gca, bluewhitered  )
% h = colorbar;
zlim([0 800]);
title('SM1')

%-----------------------------------------------------
subplot(2,2,2)
surf(lon2d, lat2d, hgt,soil_2(:,:,ics,season)); 
shading interp; view(20,50); 
xlabel('Longitude')
ylabel('Latitude')
caxis([-0.025 0.025])
colormap (gca, bluewhitered  )
% h = colorbar;
zlim([0 800]);
title('SM2')

%-----------------------------------------------------
subplot(2,2,3)
surf(lon2d, lat2d, hgt,soil_3(:,:,ics,season)); 
shading interp; view(20,50); 
xlabel('Longitude')
ylabel('Latitude')
zlabel('Surface Elevation (m)')
caxis([-0.025 0.025])
colormap (gca, bluewhitered  )

hp3 = get(subplot(2,2,3),'Position')
h = colorbar('southoutside')
set(h, 'Position', [hp3(1) + 0.1  hp3(2)-0.07  hp3(3)*1.45 0.01 ]);
h.Label.String = ['SM diff. (m^3 m^{-3})'];
% h = colorbar;
zlim([0 800]);
title('SM3')

%-----------------------------------------------------
subplot(2,2,4)
surf(lon2d, lat2d, hgt,soil_4(:,:,ics,season)); 
shading interp; view(20,50); 
xlabel('Longitude')
ylabel('Latitude')
caxis([-0.025 0.025])
colormap (gca, bluewhitered  )

zlim([0 800]);
title('SM4')

title_name = [char(case_list(ics)) , ' in ', char(season_list(season))];
suptitle(title_name)


print('-depsc', title_name);

end


end

