function analyze_this(info, dirs)

% read dimensions
d       = xb_read_dims;

% read measurements
data    = load(fullfile(dirs.data, ['data_' info.run '.mat']));

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

% plot output
xb_plot_profile(xbo, ...
    'measured',     [data.x data.zfin'], ...
    'xbeach',       [data.x data.z_finX], ...
    'other',        [data.(['x_petra' info.run(1)]) data.(['Pmod_' info.run(1)])]);

set(gca, 'XLim', [4500 5500]);

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