% script to test vaklodingen
urls={'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/vaklodingen/vaklodingenKB138_1514.nc',
    'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/vaklodingen/vaklodingenKB134_1110.nc',
    'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/vaklodingen/vaklodingenKB135_1110.nc',
    'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/vaklodingen/vaklodingenKB135_1312.nc',
    'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/vaklodingen/vaklodingenKB134_1312.nc',
    'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/vaklodingen/vaklodingenKB136_0504.nc'}
for i=1:length(urls)
    url = urls{i}
time = nc_varget(url,'time');
if ~find(time == 12784)
    error(['couldn''t find time = 12784 in ' url])
end
info = nc_info(url);
z = nc_varget(url, 'z', [0 250 250], [-1 100 100]);
if min(z(:)) < -100
    error(['found a value which is too low' num2str(min(z(:)))])
end
if max(z(:)) > 100
    error(['found a value which is too high' num2str(min(z(:)))])
end
end