function [xi, yi, altitude] = plot_interpolation(year,month,plot)
%%% plot rectified interpolation of 1 survey
% year = 2009;
% month = 12;

[cross_shore,altitude] = combine_surveys(year,month);
dist_alongshore = nc_varget('vlugtenburg.nc','dist_alongshore');

[xi,yi] = meshgrid(dist_alongshore,cross_shore);

if exist('plot')==1 && plot~=0
    figure
    plot3(xi,yi,altitude);
    
    figure
    pcolor(xi,yi,altitude);
    shading interp
end