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

fennecs/matlab/dispespicParts.m
a2847b9464fageom_testpardiso

Last Change
Raw File

dispespicParts.m
View Options

function dispespicParts(M)
%dispespicParts Show on an interactive plot the phase space of the given specie in M
% Detailed explanation goes here
fieldstep=1;
f2=uifigure('Name','Particles data');
mf2=uipanel(f2,'Position',[5 50 f2.Position(3)-10 f2.Position(4)-55]);
mf2.AutoResizeChildren='off';
m2=uipanel(f2,'Position',[5 5 f2.Position(3)-10 40]);
sgtitle(mf2,sprintf('step=%d t=%0.5e s',fieldstep*M.it2,M.tpart(fieldstep)))
sld2 = uislider(m2,'Position',[10 30 0.6*m2.Position(3) 3]);
sld2.Value=fieldstep;
sld2.Limits=[1 length(M.tpart)];
edt2 = uieditfield(m2,'numeric','Limits',[1 length(M.tpart)],'Value',1);
edt2.Position=[sld2.Position(1)+sld2.Position(3)+25 5 40 20];
edt2.RoundFractionalValues='on';
Printbt2=uibutton(m2,'Position',[edt2.Position(1)+edt2.Position(3)+10 5 40 20],'Text', 'Save');
sld2.ValueChangingFcn={@updatefigpartdata,edt2,mf2};
edt2.ValueChangedFcn={@updatefigpartdata,sld2,mf2};
Printbt2.ButtonPushedFcn={@plotPartButtonPushed};
PlotEspic2dpartdata(mf2,M,fieldstep);
function PlotEspic2dpartdata(fig,M,fieldstep)
%PlotEspic2dgriddata Plot the 2d data of espic2d at time step fieldstep
sgtitle(fig,sprintf('step=%d t=%0.5e s',(fieldstep-1)*M.it2,M.tpart(fieldstep)))
ax1=subplot(3,2,1,'Parent',fig);
plot(ax1,M.Z(1:M.nbparts(fieldstep),fieldstep),M.R(1:M.nbparts(fieldstep),fieldstep),'.');
xlim(ax1,[M.zgrid(1) M.zgrid(end)])
ylim(ax1,[M.rgrid(1) M.rgrid(end)])
xlabel(ax1,'Z [m]')
ylabel(ax1,'R [m]')
title(ax1,'Position')
grid(ax1, 'on')
ax2=subplot(3,2,2,'Parent',fig);
plot(ax2,M.Vpar(1:M.nbparts(fieldstep),fieldstep),M.Vperp(1:M.nbparts(fieldstep),fieldstep),'.');
xlabel(ax2,'V_{par} [m/s]')
ylabel(ax2,'V_{perp} [m/s]')
axis(ax2, 'equal')
grid(ax2, 'on')
ax3=subplot(3,2,3,'Parent',fig);
plot(ax3,M.Z(1:M.nbparts(fieldstep),fieldstep),M.VZ(1:M.nbparts(fieldstep),fieldstep),'.');
xlim(ax3,[M.zgrid(1) M.zgrid(end)])
xlabel(ax3,'Z [m]')
ylabel(ax3,'VZ [m/s]')
grid(ax3, 'on')
ax4=subplot(3,2,4,'Parent',fig);
plot(ax4,M.VR(1:M.nbparts(fieldstep),fieldstep),M.R(1:M.nbparts(fieldstep),fieldstep),'.')
ylabel(ax4,'R [m]')
xlabel(ax4,'VZ [m/s]')
ylim(ax4,[M.rgrid(1) M.rgrid(end)])
grid(ax4, 'on')
ax5=subplot(3,2,5,'Parent',fig);
plot(ax5,M.Z(1:M.nbparts(fieldstep),fieldstep),M.VR(1:M.nbparts(fieldstep),fieldstep),'.');
xlim(ax5,[M.zgrid(1) M.zgrid(end)])
xlabel(ax5,'Z [m]')
ylabel(ax5,'VR [m/s]')
grid(ax5, 'on')
ax6=subplot(3,2,6,'Parent',fig);
plot(ax6,M.VR(1:M.nbparts(fieldstep),fieldstep),M.R(1:M.nbparts(fieldstep),fieldstep),'.');
ylim(ax6,[M.rgrid(1) M.rgrid(end)])
ylabel(ax6,'R [m]')
xlabel(ax6,'VR [m/s]')
grid(ax6, 'on')
end
function updatefigpartdata(control, event, Othercontrol, fig)
if strcmp(event.EventName,'ValueChanged')
fieldstep=floor(control.Value);
control.Value=fieldstep;
else
fieldstep=floor(event.Value);
end
Othercontrol.Value=fieldstep;
sgtitle(fig,sprintf('t=%0.5e s',M.tpart(fieldstep)))
%% update Position plot
ax1=fig.Children(end);
ax1.Children.XData=M.Z(1:M.nbparts(fieldstep),fieldstep);
ax1.Children.YData=M.R(1:M.nbparts(fieldstep),fieldstep);
view(ax1,2)
%% update VZ VPERP
fig.Children(end-1).Children(1).XData=M.Vpar(1:M.nbparts(fieldstep),fieldstep);
fig.Children(end-1).Children(1).YData=M.Vperp(1:M.nbparts(fieldstep),fieldstep);
axis(fig.Children(end-1),'equal')
%% update Z VZ
fig.Children(end-2).Children(1).XData=M.Z(1:M.nbparts(fieldstep),fieldstep);
fig.Children(end-2).Children(1).YData=M.VZ(1:M.nbparts(fieldstep),fieldstep);
%% update VZ VR
fig.Children(end-3).Children(1).XData=M.VZ(1:M.nbparts(fieldstep),fieldstep);
fig.Children(end-3).Children(1).YData=M.R(1:M.nbparts(fieldstep),fieldstep);
%% update R VR
fig.Children(end-4).Children(1).XData=M.Z(1:M.nbparts(fieldstep),fieldstep);
fig.Children(end-4).Children(1).YData=M.VR(1:M.nbparts(fieldstep),fieldstep);
%% update Z VR
fig.Children(end-5).Children(1).XData=M.VR(1:M.nbparts(fieldstep),fieldstep);
fig.Children(end-5).Children(1).YData=M.R(1:M.nbparts(fieldstep),fieldstep);
drawnow limitrate
end
function plotPartButtonPushed(btn,ax)
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
f=figure();
PlotEspic2dpartdata(f,M,sld.Value);
f.PaperOrientation='portrait';
[~, name, ~] = fileparts(M.file);
print(f,sprintf('%sParts%d',name,sld.Value),'-dpdf','-fillpage')
end
end
c4science · Help