function analyze_this(info, dirs)

datadir = fullfile(dirs.data, info.run);

% read dimensions
d1      = xb_read_dims(datadir);
d2      = xb_read_dims;

% read measurements
lidar   = load(fullfile(dirs.data, info.run, 'lidar.xz'));
%topo    = load(fullfile(dirs.data, info.run, 'topo.xz'));

% read output
xbo1    = xs_peel(xb_read_output(datadir, 'vars', {'zb'}, 'index', [0 d1.globaltime]),'nested',true);
xbo2    = xb_read_output(pwd, 'vars', {'zb'}, 'index', [0 d2.globaltime]);

% plot output
xb_plot_profile(xbo2, ...
    'measured',     lidar, ...
    'xbeach',       [mean(xbo1.DIMS.x,1)' squeeze(mean(xbo1.zb(end,:,:),2))]    );

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