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

fennecs/matlab/dispespicFields.m
dd5453d97d74master

Last Change
Raw File

dispespicFields.m
View Options

function dispespicFields(M,logdensity,showgrid,fixed,parper)
%dispespicFields Allows to display the time evolution of the density, electric potential and electric fields
% M is of class espic2dhdf5 and contains the simulation results
fieldstep=1;
if nargin <2
logdensity=false;
showgrid=false;
fixed=false;
end
if nargin <3
showgrid=false;
fixed=false;
end
if nargin <4
fixed=false;
end
if nargin <5
parper=false;
end
fixed=fi(fixed);
f=uifigure('Name',sprintf('Grid data %s',M.name));
mf=uipanel(f,'Position',[5 50 f.Position(3)-10 f.Position(4)-55]);
mf.AutoResizeChildren='off';
m=uipanel(f,'Position',[5 5 f.Position(3)-10 40]);
sgtitle(mf,sprintf('step=%d t=%0.5e s',fieldstep*M.it1,M.t2d(fieldstep)))
sld = uislider(m,'Position',[10 30 0.6*m.Position(3) 3]);
sld.Value=fieldstep;
sld.Limits=[1 size(M.t2d,1)];
edt = uieditfield(m,'numeric','Limits',[1 size(M.t2d,1)],'Value',1);
edt.Position=[sld.Position(1)+sld.Position(3)+25 5 40 20];
edt.RoundFractionalValues='on';
MaxN=0;
Printbt=uibutton(m,'Position',[edt.Position(1)+edt.Position(3)+10 5 40 20],'Text', 'Save');
Play=uibutton(m,'Position',[Printbt.Position(1)+Printbt.Position(3)+10 5 40 20],'Text', 'Play');
Pause=uibutton(m,'Position',[Play.Position(1)+Play.Position(3)+10 5 40 20],'Text', 'Pause');
%Playbt=uibutton(m,'Position',[Printbt.Position(1)+Printbt.Position(3)+10 5 40 20],'Text', 'Play/Pause');
stop=false;
sld.ValueChangingFcn={@updatefigdata,edt,mf};
edt.ValueChangedFcn={@updatefigdata,sld,mf};
Printbt.ButtonPushedFcn={@plotGridButtonPushed};
Play.ButtonPushedFcn={@plotPlayButtonPushed};
Pause.ButtonPushedFcn={@PauseButtonPushed};
set(f,'KeyPressFcn',{ @onKeyDown,sld,edt,mf})
PlotEspic2dgriddata(mf,M,fieldstep);
function plotPlayButtonPushed(btn,ax)
stop=false;
i=sld.Value;
while ~stop
edt.Value=i;
sld.Value=i;
updatesubplotsdata(i,mf);
pause(0.01)
i=sld.Value;
i=i+10;
if(i>sld.Limits(2))
stop=true;
end
end
end
function PauseButtonPushed(btn,ax)
stop = true;
end
function onKeyDown(src,event,slider,editfield, fig)
direction=0;
if strcmp(event.Key,'leftarrow')
direction=-1;
elseif strcmp(event.Key,'rightarrow')
direction=+1;
elseif strcmp(event.Key,'uparrow')
direction=+10;
elseif strcmp(event.Key,'downarrow')
direction=-10;
end
if(direction~=0)
currval=slider.Value;
slider.Value=max(slider.Limits(1),min(currval+direction,slider.Limits(2)));
updatefigdata(slider, event, editfield ,fig)
end
end
function PlotEspic2dgriddata(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.it1,M.t2d(fieldstep)))
ax1=subplot(2,2,1,'Parent',fig);
sf=surface(ax1,M.zgrid,M.rgrid,M.N(:,:,fieldstep),'edgecolor','none');
if(showgrid)
set(sf,'edgecolor','black')
end
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')
c = colorbar(ax1);
c.Label.String= 'n[m^{-3}]';
%c.Limits=[0 max(M.N(:))];
if(isboolean(fixed) && fixed)
climits=caxis(ax1);
MaxN=climits(2);
elseif(~isboolean(fixed))
MaxN=fixed.data;
caxis(ax1,[-Inf MaxN]);
end
view(ax1,2)
hold(ax1, 'on')
contour(ax1,M.zgrid,M.rgrid,M.geomweight(:,:,1),[0 0],'r-.','linewidth',1.5,'Displayname','Boundaries');
if logdensity
set(ax1,'zscale','log')
set(ax1,'colorscale','log')
end
ax2=subplot(2,2,2,'Parent',fig);
pot=M.pot(:,:,fieldstep);
pot(M.geomweight(:,:,1)<0)=0;
contourf(ax2,M.zgrid,M.rgrid,pot);
%plot(ax2,M.zgrid,data.pot(:,5))
xlabel(ax2,'z [m]')
ylabel(ax2,'r [m]')
colormap(ax2,'jet')
c = colorbar(ax2);
c.Label.String= '\Phi [V]';
title(ax2,'Es potential')
grid(ax2, 'on')
hold(ax2, 'on')
contour(ax2,M.zgrid,M.rgrid,M.geomweight(:,:,1),[0 0],'w-.','linewidth',1.5,'Displayname','Boundaries');
Blines=zeros(size(M.Athet));
for i=1:length(M.zgrid)
Blines(i,:)=M.Athet(i,:).*M.rgrid';
end
zindex=floor(length(M.zgrid/2));
levels=logspace( log10(min(min(Blines(:,2:end)))), log10(max(max(Blines(:,:)))),8);
contour(ax2,M.zgrid,M.rgrid,Blines',real(levels),'m-.','linewidth',1.5,'Displayname','Magnetic field lines');
ax3=subplot(2,2,3,'Parent',fig);
if parper
Ez=M.Epar(fieldstep);
else
Ez=M.Ez(:,:,fieldstep);
end
Ez(M.geomweight(:,:,1)<0)=0;
contourf(ax3,M.zgrid,M.rgrid,Ez)
xlabel(ax3,'z [m]')
ylabel(ax3,'r [m]')
c = colorbar(ax3);
c.Label.String= 'E_{z} [V/m]';
if parper
title(ax3,'Parallel Electric field')
else
title(ax3,'Axial Electric field')
end
grid(ax3, 'on')
ax4=subplot(2,2,4,'Parent',fig);
if parper
Er=M.Eperp(fieldstep);
else
Er=M.Er(:,:,fieldstep);
end
Er(M.geomweight(:,:,1)<0)=0;
contourf(ax4,M.zgrid,M.rgrid,Er)
xlabel(ax4,'z [m]')
ylabel(ax4,'r [m]')
c = colorbar(ax4);
c.Label.String= 'E_{r} [V/m]';
if parper
title(ax4,'Perpendicular Electric field')
else
title(ax4,'Radial Electric field')
end
grid(ax4, 'on')
linkaxes([ax1,ax2,ax3,ax4]);
end
function plotGridButtonPushed(btn,ax)
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
f=figure();
PlotEspic2dgriddata(f,M,sld.Value);
f.PaperOrientation='landscape';
[~, name, ~] = fileparts(M.file);
print(f,sprintf('%sGrid%d',name,sld.Value),'-dpdf','-fillpage')
end
function updatefigdata(control, event, Othercontrol, fig)
if strcmp(event.EventName,'ValueChanged')
fieldstep=floor(control.Value);
control.Value=fieldstep;
elseif strcmp(event.EventName,'KeyPress')
fieldstep=floor(control.Value);
control.Value=fieldstep;
else
fieldstep=floor(event.Value);
end
Othercontrol.Value=fieldstep;
updatesubplotsdata(fieldstep, fig);
end
function updatesubplotsdata(fieldstep, fig)
sgtitle(fig,sprintf('step=%d t=%0.5e s',(fieldstep-1)*M.it1,M.t2d(floor(fieldstep))))
%% update Position histogram
ax1=fig.Children(end);
dens=M.N(:,:,fieldstep);
ax1.Children(end).ZData=dens;
ax1.Children(end).CData=dens;
if(isboolean(fixed) && fixed)
climits=caxis(ax1);
MaxN=climits(2);
nmax=max(dens(:));
if(nmax>MaxN)
MaxN=nmax;
end
caxis(ax1,[0 MaxN]);
elseif(~isboolean(fixed))
MaxN=fixed.data;
caxis(ax1,[-Inf MaxN]);
end
% view(ax1,2)
%% update ES Potential
pot=M.pot(:,:,fieldstep);
pot(M.geomweight(:,:,1)<0)=0;
fig.Children(end-2).Children(end).ZData=pot;
%% update Radial electric field
if parper
Ez=M.Epar(fieldstep);
else
Ez=M.Ez(:,:,fieldstep);
end
Ez(M.geomweight(:,:,1)<0)=0;
fig.Children(end-4).Children(1).ZData=Ez;
%% update Axial electric field
if parper
Er=M.Eperp(fieldstep);
else
Er=M.Er(:,:,fieldstep);
end
Er(M.geomweight(:,:,1)<0)=0;
fig.Children(end-6).Children(1).ZData=Er;
%drawnow limitrate
end
end
c4science · Help