%% define source
nc = 'http://dtvirt5.deltares.nl:8080/thredds/dodsC/opendap/rijkswaterstaat/kustlidar/09bn2.nc'

%% show contents
nc_dump(nc)

%% get time and coordinates
D.t   = nc_cf_time(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');

for it  = 1:length(D.t)

%% get data of per timestep

   D.z = nc_varget(nc,'z',[it-1 0 0],[1 -1 -1]);
   
%% plot per timestep (either in x/y or lat/lon)

   pcolorcorcen(D.lon,D.lat,D.z,[.5 .5 .5])
   tickmap('ll','format','%0.2f'); %tickmap('xy')
   axislat(52)
   grid on
   colorbarwithtitle('z [m]')
   title(datestr(D.t(it),'yyyy-mmm'))
   print2a4([filename(nc),'_',datestr(D.t(it),'yyyy-mmm')],'v','t')

end