function analyze_this(info, dirs)

nr      = str2num(info.run(2:4));

% read dimensions
d       = xb_read_dims;

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

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

% plot output
xb_plot_profile(xbo, ...
    'measured',     profA(nr).november, ...
    'nonerodible',  profA(nr).dijk, ...
    'xbeach',       profA(nr).xbeachT2v18, ...
    'durosta',      profA(nr).durosta                            );

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