[cwzebras]

for npt=1:nt
 
        offset = ( ( npt - 1 ) * nx * ny * nz ) * 4 + hoffset ;

        fseek(fid,offset,'bof') ;

        for j=1:ny

            for i=1:nx

                val(j,i)=fread(fid,1,fbit);

            end

        end

        pcolor(x,y,val)
        colorbar('vert');
        caxis(colrange);
        colormap colmap
        titlestr=['Snapshot of ' vnm ' at ' ... 
                  '{\fontsize{16}{\bf ' ... 
                  num2str(double(t(npt)),'%06.3f') '}}s'];
        title(titlestr);
        colorbar('vert');
        caxis(colrange);
        colormap(colmap);
        drawnow

        F  = getframe(gcf);
        im = frame2im(F);
        [im, map] = rgb2ind(im, 256);

        if npt == 1;
            imwrite(A, map, fout, 'gif', 'LoopCount', Inf, 'DelayTime', 0.1);
        else
            imwrite(A, map, fout, 'gif', 'WriteMode', 'append', 'DelayTime', 0.1);
        end

    end

上面的代码是相关的画图部分，输出的gif格式的图片是黑色的，但是运行时显示的是正确的结果（drawnow部分）。