tr = jarkus_transects('url', '..\data\jarkus.nc', 'output', {'id'});
ID = tr.id;
irun = 0;

for id = ID
    basepath = fullfile(fileparts(which(mfilename)), '..', 'simulations', sprintf('%08i', id), 'durosplus', sprintf('run%02i', irun));
    ncduros = fullfile(basepath, sprintf('%08i_pp.nc', id));
    EVduros(:,id==ID) = nc_varget(ncduros, 'erosion_volume');
    
    ncxbeach = strrep(ncduros, 'durosplus', 'xbeach');
    EVxbeach(:,id==ID) = nc_varget(ncxbeach, 'erosion_volume');
end