%
% function [M] = show_movie(field,x,y,NREC,lax,palette,NFIG,NTIMES,FPS)
%
% return a movie matrix M.
%
    function [M] = show_movie(field,x,y,NREC,lax,palette,NFIG,NTIMES,FPS)

    M = moviein(NREC);
    figure(NFIG)
    for n=1:NREC
      pcolor(x,y,field(:,:,n)');
      title(['record N = ' num2str(n)],'Fontsize',18);
      caxis(lax);
      colormap(palette);
      colorbar;
      M(:,n) = getframe;
    end
    movie(M,NTIMES,FPS)