create_gif.m
f916da899a9bPPB110

create_gif.m
View Options

			title1 = GIFNAME;
			FIGDIR = BASIC.RESDIR;
			GIFNAME = [FIGDIR, GIFNAME,'.gif'];

			% Set colormap boundaries
			hmax = max(max(max(FIELD)));
			hmin = min(min(min(FIELD)));

			flag = 0;
			if hmax == hmin
			disp('Warning : h = hmin = hmax = const')
			else
			% Setup figure frame
			fig = figure('Color','white','Position', [100, 100, 400, 400]);
			pcolor(X,Y,FIELD(:,:,1)); % to set up
			colormap jet
			axis tight manual % this ensures that getframe() returns a consistent size
			if INTERP
			shading interp;
			end
			in = 1;
			nbytes = fprintf(2,'frame %d/%d',in,numel(FIELD(1,1,:)));
			for n = FRAMES % loop over selected frames
			scale = max(max(abs(FIELD(:,:,n))));
			pclr = pcolor(X,Y,FIELD(:,:,n)/scale); % frame plot
			if INTERP
			shading interp;
			end
			set(pclr, 'edgecolor','none'); axis square;
			% caxis([min(min(FIELD(:,:,n))),max(max(FIELD(:,:,n)))]);
			xlabel(XNAME); ylabel(YNAME); %colorbar;
			title([FIELDNAME,', $t \approx$', sprintf('%.3d',ceil(T(n)))...
			,', scaling = ',sprintf('%.1e',scale)]);
			drawnow
			% Capture the plot as an image
			frame = getframe(fig);
			im = frame2im(frame);
			[imind,cm] = rgb2ind(im,32);
			% Write to the GIF File
			if in == 1
			imwrite(imind,cm,GIFNAME,'gif', 'Loopcount',inf);
			else
			imwrite(imind,cm,GIFNAME,'gif','WriteMode','append', 'DelayTime',DELAY);
			end
			% terminal info
			while nbytes > 0
			fprintf('\b')
			nbytes = nbytes - 1;
			end
			nbytes = fprintf(2,'frame %d/%d',n,numel(FIELD(1,1,:)));
			in = in + 1;
			end
			disp(' ')
			disp(['Gif saved @ : ',GIFNAME])
			end