function plotresults(info,dirs)

v = {'H' 'thetamean' 'Fy_mean'};
n = [1 3 5];
a = [0 10 20];

for i = 1:length(v)
    figure;
    for j = 1:length(n)
        for k = 1:length(a)
            run = sprintf('%1dbin_alpha%02d',n(j),a(k));
            
            disp(run);
            
            subplot(length(n),length(a),(j-1)*length(n)+k); hold on;
            
            try
                xbd = xb_read_dims(run);
                
                if regexp(v{i},'_mean$')
                    it = 0;
                else
                    it = xbd.nt-1;
                end
                
                xbo = xb_read_output(run,'vars',v{i},'start',it,'length',1);
                d   = squeeze(xs_get(xbo,v{i}));

                pcolor(xbd.x,xbd.y,d);
                
                shading flat;
                axis image;
                
                set(gca,'XTick',[],'YTick',[]);
            end
            
            box on;
            
            if j == 1
                title(sprintf('%d^o',a(k)));
            end
            
            if k == 1
                ylabel(sprintf('%1d bin',n(j)));
            end
        end
    end
    xb_write_plot(gcf,dirs.output,v{i});
end
            