 year = 2009;
 month = 7:12;

for i = 1:length(month)

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

[along_i cross_i] = meshgrid(along_shore,cross_shore);

pcolor(along_i,cross_i,altitude)
shading interp
axis equal
pause

end

year = 2010;
month = 1:3;

for i = 1:length(month)

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

[along_i cross_i] = meshgrid(along_shore,cross_shore);

pcolor(along_i,cross_i,altitude)
shading interp
axis equal
pause

end