function analyze_this(info, dirs)

re      = regexp(info.run,'_','split');

% read dimensions
d       = xb_read_dims;

% read measurements
comp    = load(fullfile(dirs.data, re{1}, re{2}, 'results.mat'));
meas    = load(fullfile(dirs.data, re{1}, re{2}, 'results_meas.mat'));

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

% plot output
xb_plot_profile(xbo, ...
    'measured',     [-meas.xend_meas./1.0401 meas.zend_meas], ...
    'xbeach',       [-comp.x./1.0401 comp.zend]                      );

set(gca, 'XLim', [-4677 -4300]);

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