function analyze_this(info, dirs)

t       = [1 2 4 8];
d       = xb_read_dims;
[mn i]  = closest(1+t*3600, d.t);

% read measurements
load(fullfile(dirs.data, '1C.mat'));

% read output
xbo1    = xb_read_output(pwd, 'vars', {'zb'}, 'index', i-1);
xbo2    = xb_read_output(pwd, 'vars', {'zb' 'zs' 'H' 'u*' 'hh_var'});

% analyze output
xbw     = xb_get_hydro(xbo2, 'Trep', 8);
xbb     = xb_get_morpho(xbo2, 'level', 0);

% plot output
xb_plot_hydro(xbw, ...
    's',        [s.xps;s.setup]',    ...
    'Hrms_lf',  [s.xps;s.Hrms_hf]',  ...
    'Hrms_hf',  [s.xps;s.Hrms_lf]',  ...
    'Hrms_t',   [s.xps;s.Hrms]',     ...
    'urms_lf',  [s.xps;s.urms_lf1]', ...
    'urms_t',   [s.xps;s.urms1]'                );

set(gca,'XLim',[0 196]);

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

xb_plot_morpho(xbb);

set(subplot(3,1,1),'XLim',[0 200]);
set(subplot(3,1,2),'XLim',[0 8]);

pause(1); % this should be here because of a bug in Matlab (see service request 1-F227KF)

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

% write skill table
xb_write_skilltable(                                               ...
    xs_join(xbw, xbb),                                             ...
    {[s.xps;s.Hrms_hf]' [s.xps;s.Hrms_lf]' [s.xps;s.urms1]'},      ...
    'vars', {'Hrms_hf' 'Hrms_lf' 'urms_t'},                        ...
    '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', [s.x{1};cell2mat(s.z(end))]')
title(['profiles ' info.test ' ' info.run], 'interpreter', 'none')
xb_write_plot(gcf, dirs.output, 'final_profiles')