function [] = plotsteps()
load stepblocks3mer
%load stepblocks3mer_cov30
n = 6;
b1 = 80;
b2 = 6;
j = [1 4 7 10 2 5 8 11 3 6 9 12 13 14 16 17];
C = ['bbbbrrrrmmmmcccc'];
for i = 1:16
subaxis(6,6,2*j(i)-1, 'SpacingVert',0.02,'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([reshape(ma3(i,:,:),6,6) zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b1 b1]);
axis ij square, colorbar, shading flat;
title(strcat(StepA(i,:), ' averages'), 'Color', C(i));
subaxis(6,6,2*j(i), 'SpacingVert',0.02, 'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([reshape(sa3(i,:,:),6,6) zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b2 b2]);
axis ij square, colorbar, shading flat;
title(strcat(StepA(i,:), ' std (', int2str(na3(i)),' instances)'), 'Color', C(i));
end;
% -----average XAY block-----
subaxis(6,6,2*17+1, 'SpacingVert',0.02,'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([ma3all zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b1 b1]);
axis ij square, colorbar, shading flat;
title('All XAY averages','Color','k');
subaxis(6,6,2*18, 'SpacingVert',0.02, 'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([sa3all zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b2 b2]);
axis ij square, colorbar, shading flat;
title(strcat('All XAY std (', int2str(na3all),' instances)'),'Color','k');
%-----------------------------
%if 0
figure()
for i = 1:16
subaxis(6,6,2*j(i)-1, 'SpacingVert',0.02,'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([reshape(mg3(i,:,:),6,6) zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b1 b1]);
axis ij square, colorbar, shading flat;
title(strcat(StepG(i,:), ' averages'), 'Color', C(i));
subaxis(6,6,2*j(i), 'SpacingVert',0.02, 'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([reshape(sg3(i,:,:),6,6) zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b2 b2]);
axis ij square, colorbar, shading flat;
% title(strcat(StepG(i,:), ' std'));
title(strcat(StepG(i,:), ' std (', int2str(ng3(i)),' instances)'), 'Color', C(i));
end;
% -----average XAY block-----
subaxis(6,6,2*17+1, 'SpacingVert',0.02,'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([mg3all zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b1 b1]);
axis ij square, colorbar, shading flat;
title('All XGY averages','Color','k');
subaxis(6,6,2*18, 'SpacingVert',0.02, 'SpacingHoriz',0.02,'MarginTop',0.02,'MarginBottom',0.02,'MarginLeft',0.02,'MarginRight',0.02);
pcolor([sg3all zeros(n,1);zeros(1,n+1)]);
caxis manual,caxis([-b2 b2]);
axis ij square, colorbar, shading flat;
title(strcat('All XGY std (', int2str(ng3all),' instances)'),'Color','k');
%-----------------------------
%end;
end