function analyze_this_uniform(info, dirs)

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

% read dimensions
d       = xb_read_dims;
t       = [0 0.06 1.3 3.5 9.6];
[mn i]  = closest(t*3600, d.t);

% read input
xbi     = xb_read_input(pwd);
zs0     = max([0 xs_get(xbi, 'zs0') max(xs_get(xbi, 'zs0file.tide'))]);

if xs_exist(xbi, 'ne_layer')
    ne  = [d.x(1,:)',max(xs_get(xbi, 'depfile.depfile')-xs_get(xbi, 'ne_layer.ne_layer'),-2)'];
else
    ne  = [];
end

% read output
xbo1    = xb_read_output(pwd, 'vars', {'zb'}, 'index', i);
xbo2    = xb_read_output(pwd, 'vars', {'zb'});
xbo3    = xb_read_output(pwd, 'vars', {'zs' 'zb' 'H' 'u' 'ue' 'urms' 'hh_var'});

% analyze output
xbb     = xb_get_morpho(xbo2, 'level', zs0);
xbw     = xb_get_hydro(xbo3, 'Trep', 1/xs_get(xbi, 'bcfile.fp'));

% plot output
xb_plot_profile(xbo1,                               ...
    'measured',     [m.x',cat(1,m.z{:})'-4.2],      ...
    'nonerodible',  ne                                  );

axis([150 210 -2 2.5]); 

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

xb_plot_morpho(xbb);

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

xb_plot_hydro(xbw, ...
    'umean', [h.xm' h.umm']);

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

% write skill table
xb_write_skilltable(                                        ...
    xs_join(xbo1, xbb),                                     ...
    {[m.x',cat(1,m.z{:})'-zs0] [h.xm' h.umm']},                 ...
    'vars', {'zb' 'umean'},  ...
    'file', fullfile(dirs.output, 'skills.tex')                 );

xb_plot_profile(xbo1,'measured', [m.x',cat(1,m.z{end})'-4.2], 'nonerodible', ne)
title(['profiles ' info.test ' ' info.run])
xb_write_plot(gcf, dirs.output, 'final_profiles')