function analyze_this(info, dirs)

runs    = {'event1' 'event2' 'event3' 'event8' 'event9'};
nr      = find(strcmpi(runs,info.run));

% read dimensions
d       = xb_read_dims;

% read measurements
data    = load(fullfile(dirs.data, 'testbed.mat'));

% read output
xbo     = xb_read_output(pwd, 'vars', {'zb'}, 'index', [0 d.globaltime-1]);

% plot output
xb_plot_profile(xbo, ...
    'measured',     [data.topo(nr).xf(:) data.topo(nr).zf(:)], ...
    'other',        [data.litpr(nr).x(:) data.litpr(nr).z(:)], ...
    'xbeach',       [data.topo(nr).x(:) data.xb(nr).zmf(:)]         );

set(gca,'XLim',[1200 1640]);

xb_write_plot(gcf, dirs.output, 'profiles');