videodens.m
0b252c046024master

videodens.m
View Options

			resfile='/misc/nonlinear/results_test_bias1kV.h5';

			M=nonlinhdf5(resfile);

			r=M.rgrid;
			thet=M.thetgrid;
			thet=thet2pi/max(thet);
			[R,THET]=meshgrid(r,thet);
			X=R.*cos(THET);
			Y=R.*sin(THET);

			O=[0 0];
			ra=M.ra;
			rb=M.rb;
			rm=1.73e-2;
			rp=1.91e-2;
			n0=M.n0;

			n=-M.rho(:,:,:)/1.602e-19;
			t=M.t2d;

			video=VideoWriter('film_test_bias10kV');
			%video.FrameRate=10;
			video.Quality=100;
			open(video)

			figure('Position', [100, 100, 700, 600])
			for i=1:20:length(t)
			pcolor(X,Y,n(:,:,i)')
			shading interp
			colormap(flipud(hot))
			colorbar
			caxis([0 1.5*n0])
			%hold on
			%rectangle('Position',[O(1)-ra O(2)-ra 2ra 2ra],'Curvature',[1 1],'LineWidth',1)
			%hold on
			%rectangle('Position',[O(1)-rb O(2)-rb 2rb 2rb],'Curvature',[1 1],'LineWidth',1)
			%hold on
			%rectangle('Position',[O(1)-rm O(2)-rm 2rm 2rm],'Curvature',[1 1],'EdgeColor','k','LineWidth',0.5,'LineStyle','--')
			%hold on
			%rectangle('Position',[O(1)-rp O(2)-rp 2rp 2rp],'Curvature',[1 1],'EdgeColor','k','LineWidth',0.5,'LineStyle','--')
			xlabel('x [m]')
			ylabel('y [m]')
			title(['n_e (t = ', num2str(t(i)),' s)']);
			axis equal

			drawnow
			currframe=getframe(gcf);
			writeVideo(video,currframe)

			% GIF
			%if i == 1
			% imwrite(currframe.cdata(:,:,1), 'animation.gif', 'gif', 'Loopcount', inf, 'DelayTime', 0.1);
			%else
			% imwrite(currframe.cdata(:,:,1), 'animation.gif', 'gif', 'WriteMode', 'append', 'DelayTime', 0.1);
			%end
			end
			close(video)
			close(gcf)