figure;

LON = []; LAT = []; Z = [];

for i = 119:121
    url = ['P:\mcdata\opendap\rijkswaterstaat\vaklodingen\vaklodingenKB' num2str(i) '_2726.nc'];

    info = nc_info(url);

    z = nc_varget(url, 'z', [info.Dimension(1).Length-1 1 1], [1 -1 -1]);
    x = nc_varget(url, 'x');
    y = nc_varget(url, 'y');
    
    x = x(2:end)-10;
    y = y(2:end)-10;
    
    [X Y] = meshgrid(x, y);
    [lon lat] = convertCoordinates(X,Y,'CS1.code',28992,'CS2.code',4326);
    
    LON = [LON lon];
    LAT = [LAT lat];
    Z = [Z z];
end

pcolor(LON, LAT, Z);
shading flat; axis equal; grid on; colorbar;

% lon = 4.xx
% lat = 52.xx