function analyze_this(info, dirs)

% define runs
runs = {'profA', 'profB1', 'profB2', 'profC'};

fh = figure;
for i = 1:length(runs)
    subplot(2,2,i); hold on;
    
    runx = runs{i};
    if exist(runx, 'dir')
        cd(fullfile(dirs.output,runx));
        
        % read measurements
        load(fullfile(dirs.data, ['poststormJim_' runx '.mat']));
        plot(xd, yd, '-k', 'LineWidth', 2);
        
        % read dimensions
        d       = xb_read_dims;
        j       = ceil(d.globaly/2);
        x       = squeeze(d.x(:,j));

        % read output
        xbo     = xs_peel(xb_read_output(pwd, 'vars', {'zb'}, 'index', [0 d.nt]));
        
        ind = min(find(xbo.zb(1,j,:)>=-9));
        plot(x-x(ind), squeeze(xbo.zb(1,j,:)), ':k', 'LineWidth', 2);
        plot(x-x(ind), squeeze(xbo.zb(end,j,:)), '-r', 'LineWidth', 2);
        
        title(strrep(runx, 'prof', 'Profile'));
        xlabel('x (m)');
        ylabel('Elevation (m +NAVD88)');
        set(gca, 'XLim', [700 1450], 'YLim', [-3 3.5]);
        
        cd('..');
    end
end

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

% create conditions table
pname = fullfile(dirs.output, 'conditions.tex');
matrix2latex([.75 1; 4.1 3.9; 8.5 8.5], 'filename', pname, ...
    'caption', 'Hydrodynamic boundary conditions XBeach simulations', ...
    'rowlabel', {'Surge level [m +NAVD]', '$H_s$ [m]', '$T_p$ [s]'}, ...
    'collabel', {'Storm 1','Storm 2'}, ...
    'format', {'%s', '%2.1f', '%2.1f'});