Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (geom_testpardiso) Browse

fennecs/matlab/PlotParticleTrajectory.m
f7ad604a87f5geom_testpardiso

Last Change
Raw File

PlotParticleTrajectory.m
View Options

function PlotParticleTrajectory(M,partnumber,t,text)
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
if (nargin<4)
text='';
else
text=['_',text];
end
f=figure();
%time=((t(1):t(end))-1)*M.dt*double(M.it2);
time=M.tpart(t(1):t(end));
%sgtitle(sprintf('Particles %d',partnumber));
ax1=subplot(1,2,1);
hold(ax1);
ax2=subplot(1,2,2);
hold(ax2);
Msize=8;
i=1;
for part=partnumber
p1(i)=plot(ax1,M.Z(part,t),M.R(part,t),'x-.','Linewidth',1.1,...
'Displayname',sprintf('part=%d',part),'MarkerIndices',1,...
'Markersize',Msize);
plot(ax1,M.Z(part,t(end)),M.R(part,t(end)),'^','Linewidth',1.1,...
'Displayname',sprintf('part=%d',part),'MarkerIndices',1,...
'Markersize',Msize,'color',p1(i).Color)
x=M.R(part,t).*cos(M.THET(part,t));
y=M.R(part,t).*sin(M.THET(part,t));
%pp=spline(x,y);
p2=plot(ax2,x,y,'x-.','Linewidth',1.1,...
'Displayname',sprintf('part=%d',part),'MarkerIndices',1,...
'Markersize',Msize);
plot(ax2,M.R(part,t(end)).*cos(M.THET(part,t(end))),M.R(part,t(end)).*sin(M.THET(part,t(end))),'^',...
'Linewidth',1.1,'Displayname',sprintf('part=%d',part),'MarkerIndices',1,...
'Markersize',Msize,'color',p2.Color)
i=i+1;
end
xlabel(ax1,'z [m]')
ylabel(ax1,'r [m]')
sgtitle(sprintf('Position between t=[%3.2f-%3.2f] ns',time(1)*1e9,time(end)*1e9))
xlim(ax1,[M.zgrid(1) M.zgrid(end)])
%xlim(ax1,[-.17 .17])
ylim(ax1,[M.rgrid(1) M.rgrid(end)])
%ylim(ax1,[0.005 0.025])
legend(ax1,p1(1:end))
grid(ax1,'on');
xlabel(ax2,'x [m]')
ylabel(ax2,'y [m]')
t=2*pi*linspace(0,1,100);
xin=M.rgrid(1)*cos(t);
yin=M.rgrid(1)*sin(t);
plot(ax2,xin,yin,'k-')
xout=M.rgrid(end)*cos(t);
yout=M.rgrid(end)*sin(t);
plot(ax2,xout,yout,'k-')
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.fullpath);
f.PaperUnits='centimeters';
f.PaperSize=[16 9];
print(f,sprintf('%sTrajectories%s',name,text),'-dpdf','-fillpage')
end
c4science · Help