function analyze_this(info, dirs)

% read dimensions
d   = xb_read_dims;

% read output
xbo = xs_peel( ...
    xb_read_output(pwd, 'vars', {'zb' 'u' 'v'}, 'start', d.globaltime-1, 'length', 1));

% plot output
fh = figure;

surf(d.x,d.y,squeeze(xbo.zb)); hold on;

shading flat;
axis equal;
caxis([-10 2]);
daspect([10 10 1]);
colorbar;

dnx = 1;
dny = 3;

quiver( ...
    d.x(1:dnx:end,1:dny:end), ...
    d.y(1:dnx:end,1:dny:end), ...
    squeeze(xbo.u(1,1:dnx:end,1:dny:end)), ...
    squeeze(xbo.v(1,1:dnx:end,1:dny:end)),1,'w');

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