create_gif_3D.m
d1b82f9414b6gather_output

create_gif_3D.m
View Options

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

			% Setup figure frame
			fig = figure('Color','white','Position', [100, 100, 400, 400]);
			plt = @(x) squeeze(reshape(x(:,:,1),[],1));
			scale_x = max(abs(plt(X(:,:,1))));
			scale_y = max(abs(plt(Y(:,:,1))));
			scale_z = max(abs(plt(Z(:,:,1))));
			plot3(plt(X)/scale_x,plt(Y)/scale_y,plt(Z)/scale_z,'.k','MarkerSize',MARKERSIZE);
			view(VIEW);
			axis tight manual % this ensures that getframe() returns a consistent size
			in = 1;
			nbytes = fprintf(2,'frame %d/%d',in,numel(FRAMES));
			for n = FRAMES % loop over selected frames
			plt = @(x) squeeze(reshape(x(:,:,n),[],1));
			scale_x = max(abs(plt(X)));
			scale_y = max(abs(plt(Y)));
			scale_z = max(abs(plt(Z)));
			plot3(plt(X)/scale_x,plt(Y)/scale_y,plt(Z)/scale_z,'.k','MarkerSize',MARKERSIZE);
			view(VIEW);
			xlabel(XNAME); ylabel(YNAME); zlabel(ZNAME);
			title(['$t \approx$', sprintf('%.3d',ceil(T(n)))...
			,', scaling = ',sprintf('%.1e',scale_x)]);
			grid on;
			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(FRAMES));
			in = in + 1;
			end
			disp(' ')
			disp(['Gif saved @ : ',GIFNAME])