function analyze_this(info, dirs)

% read dimensions
d       = xb_read_dims;

% read measurements
x_comp  = mean(load(fullfile(dirs.data, info.run, 'x_comp.grd')),1)';
x_meas  = mean(load(fullfile(dirs.data, info.run, 'x_meas.grd')),1)';
z_comp  = mean(load(fullfile(dirs.data, info.run, 'z_comp.dat')),2);
z_meas  = mean(load(fullfile(dirs.data, info.run, 'z_meas.dep')),1)';

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

% plot output
xb_plot_profile(xbo, ...
    'measured',     [-x_meas+20 z_meas], ...
    'xbeach',       [-x_comp+20 z_comp]    );

set(gca,'XLim',[-1100 -600]);

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