function analyze_this(info, dirs)

% read measurements
pre     = load('pre.obs');
post    = load('post.obs');
x       = sort(unique([pre(:,1);post(:,1)]));
M       = [-x interp1(pre(:,1),pre(:,2),x) interp1(post(:,1),post(:,2),x)];
t       = [0 42];

% read dimensions
d       = xb_read_dims;
[m i]   = closest(t*3600, d.t);

% read output
xbo1    = xb_read_output(pwd, 'vars', {'zb'}, 'index', i-1);

% plot output
xb_write_skilltable(xbo1, {M(:,[1 end])}, 'vars', 'zb', 'file', fullfile(dirs.output, 'skills.tex'));

% remove all but the first and last outputs
for iv = 2:length(xbo1.data)
    if length(size(xbo1.data(iv).value)) == 3
        xbo1.data(iv).value = xbo1.data(iv).value([1 end],:,:);
    end
end

xb_plot_profile(xbo1,'measured',M)
title(['profiles ' info.test ' ' info.run], 'interpreter', 'none')
xb_write_plot(gcf, dirs.output, 'final_profiles')