function analyze_this(info, dirs)

xbd = xb_read_dims;
xbo = xs_peel(xb_read_output(pwd,'vars',{'H' 'u' 'v'}));
xbm = xs_peel(xb_read_output(pwd,'vars',{'H_mean' 'u_mean' 'v_mean'}));

%% plot last time step

figure; hold on;

pcolor(xbd.x,xbd.y,squeeze(xbo.H(end,:,:)));
shading interp; caxis([0 2]); colorbar;

quiver(xbd.x,xbd.y,squeeze(xbo.u(end,:,:)),squeeze(xbo.v(end,:,:)),2,'k');
axis equal;

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

%% plot mean

figure; hold on;

pcolor(xbd.x,xbd.y,squeeze(xbm.H_mean(end,:,:)));
shading interp; caxis([0 2]); colorbar;

quiver(xbd.x,xbd.y,squeeze(xbm.u_mean(end,:,:)),squeeze(xbm.v_mean(end,:,:)),2,'k');
axis equal;

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