function analyze_this_2d(info, dirs)

% read dimensions
d       = xb_read_dims;
nr      = str2num(info.run(2:4));

% read measurements
load(fullfile(dirs.data, 'profA.mat'));
load(fullfile(dirs.data, 'P.mat'));

% plot overview
figure; pcolor(P.X,P.Y,P.zb(:,:,1));
shading flat; hold on;
plot(P.s113a_grid(:,1),P.s113a_grid(:,2),'-k');
plot(P.s114a_grid(:,1),P.s114a_grid(:,2),'-k');
plot(P.s115a_grid(:,1),P.s115a_grid(:,2),'-k');
plot(P.s116a_grid(:,1),P.s116a_grid(:,2),'-k');
plot(P.s117a_grid(:,1),P.s117a_grid(:,2),'-k');
xlabel('distance [m]'); ylabel('distance [m]');
title('transect overview');

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

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

% plot output
for i = 113:117

    m   = P.(['s' num2str(nr) 'a_grid']);
    z   = nan(size(zb,1),1,size(m,1));
    
    for j = 1:size(zb,1)
        z(j,1,:) = griddata(d.x,d.y,zb(j,:,:),m(:,1),m(:,2));
    end
    
    xbp = xs_set([], 'DIMS', xs_set([], 'DIMS.globalx_DATA', sqrt(m(:,1).^2+m(:,2).^2)), 'zb', z);

    xb_plot_profile(xbp, ...
        'measured',     profA(nr).november, ...
        'nonerodible',  profA(nr).dijk, ...
        'xbeach',       P{nr}.XB2D, ...
        'durosta',      profA(nr).durosta                            );

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