2: f = PetscBinaryRead('perfectmovie','cell');
3: s = size(f); s = s(2);
4: clear M
5: for i=1:s,
6: cv = f{i};
7: cv = reshape(cv(2:5:5*96*96),96,96);
8: pcolor(cv);
9: caxis([0 1])
10: shading interp;
11: colorbar;
12: M((i-1)/1 + 1) = getframe;
13: end
14: movie(M)