PlotParticleTrajectory.m
bf4048137560geom_test

PlotParticleTrajectory.m
View Options

			function PlotParticleTrajectory(M,partnumber,t)
			%UNTITLED3 Summary of this function goes here
			% Detailed explanation goes here
			f=figure();
			time=((t(1):t(end))-1)M.dtdouble(M.it2);
			%sgtitle(sprintf('Particles %d',partnumber));
			for part=partnumber
			ax1=subplot(1,2,1);
			hold on
			plot(ax1,M.Z(part,t),M.R(part,t),'x','Linewidth',1.2,'Displayname',sprintf('part=%d',part),'MarkerIndices',1:10:length(t));
			ax2=subplot(1,2,2);
			hold on
			plot(ax2,M.R(part,t).cos(M.THET(part,t)),M.R(part,t).sin(M.THET(part,t)),'x-.','Linewidth',1.2,'Displayname',sprintf('part=%d',part),'MarkerIndices',1:10:length(t))

			% ax=subplot(1,2,2);
			% hold on
			% V=M.VR(part,t);
			% plot(ax,time,V','Displayname',sprintf('part=%d',part))
			end
			xlabel(ax1,'Z')
			ylabel(ax1,'R')
			sgtitle(sprintf('Position between t=[%3.2f-%3.2f] ns',time(1)1e9,time(end)1e9))
			xlim(ax1,[M.zgrid(1) M.zgrid(end)])
			ylim(ax1,[M.rgrid(1) M.rgrid(41)])
			legend(ax1)
			grid(ax1,'on');
			xlabel(ax2,'X')
			ylabel(ax2,'Y')
			axis(ax2,'equal')
			legend(ax2)
			grid(ax2,'on');

			% ax=subplot(1,2,2);
			% xlabel(ax,'t [s]')
			% ylabel(ax,'VR [m/s]')
			% title(ax,'Radial velocity')
			% grid on;

			f.PaperOrientation='landscape';
			[~, name, ~] = fileparts(M.file);
			print(f,sprintf('%sTrajectories',name),'-dpdf','-fillpage')

			end