function [] = func_plot_if_lateral(lon2d, lat2d, hgt,soil_1, sm_label,season)
% Summary of this function goes here
%   
make_it_tight = true;
subplot = @(m,n,p) subtightplot (m, n, p, [0.045 0.075], [0.1 0.025], [0.05 0.05]);
if ~make_it_tight,  clear subplot;  end

figure
set(gcf,'PaperPositionMode', 'auto')
set(gcf,'Position', [100 100 900 800])
%-----------------------------------------------------
subplot(3,2,1)
surf(lon2d, lat2d, hgt,soil_1(:,:,1,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]);

%-----------------------------------------------------
subplot(3,2,3)
surf(lon2d, lat2d, hgt,soil_1(:,:,2,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]);

%-----------------------------------------------------
subplot(3,2,4)
surf(lon2d, lat2d, hgt,soil_1(:,:,3,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]);

%-----------------------------------------------------
subplot(3,2,5)
surf(lon2d, lat2d, hgt,soil_1(:,:,4,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(3,2,5),'Position')
h = colorbar('southoutside')
set(h, 'Position', [hp3(1) + 0.1  hp3(2)-0.05  hp3(3)*1.45 0.01 ]);
h.Label.String = [sm_label,' (m^3 m^{-3})'];

zlim([0 800]);

%-----------------------------------------------------
subplot(3,2,6)
surf(lon2d, lat2d, hgt,soil_1(:,:,5,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]);

end

