function OK = convert_test

nc        = 'F:\opendap\thredds\rijkswaterstaat\vaklodingen\vaklodingenKB120_1918.nc'
D.x       = nc_varget(nc,'x');
D.y       = nc_varget(nc,'y');
D.lon     = nc_varget(nc,'lon');
D.lat     = nc_varget(nc,'lat')
[D.X,D.Y] = meshgrid(D.x,D.y);
[X,Y]     = convertCoordinates(D.lon,D.lat,'CS1.code',4326,'CS2.code',28992);

figure;pcolorcorcen(D.X-X);colorbarwithtitle('\epsilon_x');
figure;pcolorcorcen(D.Y-Y);colorbarwithtitle('\epsilon_y');


OK = abs(max(D.X(:)-X(:))) < 1 & ...
     abs(min(D.X(:)-X(:))) < 1 & ...
     abs(max(D.Y(:)-Y(:))) < 1 & ...
     abs(min(D.Y(:)-Y(:))) < 1   ;  