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)