figure
for i=1:somRows*somCols
  subplot(somRows,somCols,i)
  h=pcolor(squeeze(ptrns(i,:,:)));
  set(h,'edgecolor','none')
  caxis([-1.5 1.5])
end
ht=suptitle('SOM patterns');
set(ht,'interpreter','none')
print(['/home/nsteiger/Dropbox/enso-volcanoes/sompatterns_' num2str(somRows*somCols) '_phyda_1-1925.png'],'-dpng','-r300');

figure
plot(syrs,bmus,'o')
title('SOM BMUs')
print(['/home/nsteiger/Dropbox/enso-volcanoes/sombmus_' num2str(somRows*somCols) '_phyda_1-1925.png'],'-dpng','-r300');




