function [xi, yi, altitude] = plot_interpolation_nr(survey_ind,pair_nr)
%%% plot rectified interpolation of 1 survey
% survey_ind = list_surveys; 

%year = 2009;
% month = 12;

[cross_shore,altitude] = combine_surveys_nr(survey_ind,pair_nr);
dist_alongshore = nc_varget('vlugtenburg.nc','dist_alongshore');

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

plot=1 ;

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