function analyze_this_2d(info, dirs)

% read output
xbo     = xb_read_output(pwd, 'vars', {'zb'});

transects = struct( ...
    'pr04', 10, ...
    'pr18', 80, ...
    'pr20', 95          );

f = fieldnames(transects);

% plot overview
x       = load('x_2d.grd');
y       = load('y_2d.grd');
z       = load('bed_2d.dep');

figure; pcolor(x,y,z); shading flat; axis equal; hold on;
for i = 1:length(f)
    plot(x(transects.(f{i}),:),y(transects.(f{i}),:),'-k');
end

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

% plot output
for i = 1:length(f)
    
    t       = f{i};
    nr      = transects.(t);
    
    % read measurements
    x_comp  = mean(load(fullfile(dirs.data, t, 'x_comp.grd')),1)';
    x_meas  = mean(load(fullfile(dirs.data, t, 'x_meas.grd')),1)';
    z_comp  = mean(load(fullfile(dirs.data, t, 'z_comp.dat')),2);
    z_meas  = mean(load(fullfile(dirs.data, t, 'z_meas.dep')),1)';
    
    xbt = xb_get_transect(xbo, 'transect', nr);

    xb_plot_profile(xbt, ...
        'measured',     [x_meas z_meas], ...
        'xbeach',       [x_comp z_comp]    );

    xb_write_plot(gcf, dirs.output, ['profile_' t]);
end