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

fennecs/matlab/dispespicWell.m
d8035b6d442amaster

Last Change
Raw File

dispespicWell.m
View Options

function dispespicWell(M,logdensity,showgrid,fixed)
%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
fixed=fi(fixed);
f=uifigure('Name',sprintf('Well 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;
Minwell=0;
Maxwell=0;
plotaxes=gobjects(2);
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};
PlotEspic2dgriddata(mf,M,fieldstep);
function plotPlayButtonPushed(btn,ax)
stop=false;
for i=edt.Value:sld.Limits(2)
edt.Value=i;
sld.Value=i;
updatesubplotsdata(i,mf);
pause(0.01)
if stop
stop=false;
break;
end
end
end
function PauseButtonPushed(btn,ax)
stop = true;
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)))
plotaxes(1)=subplot(1,2,1,'Parent',fig);
sf=surface(plotaxes(1),M.zgrid,M.rgrid,M.N(:,:,fieldstep),'edgecolor','none');
if(showgrid)
set(sf,'edgecolor','black')
end
xlim(plotaxes(1),[M.zgrid(1) M.zgrid(end)])
ylim(plotaxes(1),[M.rgrid(1) M.rgrid(end)])
xlabel(plotaxes(1),'z [m]')
ylabel(plotaxes(1),'r [m]')
title(plotaxes(1),'Density')
c = colorbar(plotaxes(1));
c.Label.String= 'n[m^{-3}]';
%c.Limits=[0 max(M.N(:))];
if(isboolean(fixed) && fixed)
climits=caxis(plotaxes(1));
MaxN=climits(2);
elseif(~isboolean(fixed))
MaxN=fixed.data;
caxis(plotaxes(1),[-Inf MaxN]);
end
view(plotaxes(1),2)
if logdensity
set(plotaxes(1),'zscale','log')
set(plotaxes(1),'colorscale','log')
end
plotaxes(2)=subplot(1,2,2,'Parent',fig);
well=-PotentialWell(M.pot(:,:,fieldstep),M.Athet(:,:),M.rgrid,M.zgrid,false);
sf=surface(plotaxes(2),M.zgrid(2:end),M.rgrid(2:end),well(2:end,2:end),'edgecolor','none');
if(showgrid)
set(sf,'edgecolor','black')
end
xlim(plotaxes(2),[M.zgrid(1) M.zgrid(end)])
ylim(plotaxes(2),[M.rgrid(1) M.rgrid(end)])
xlabel(plotaxes(2),'z [m]')
ylabel(plotaxes(2),'r [m]')
title(plotaxes(2),'Well')
c = colorbar(plotaxes(2));
c.Label.String= 'depth [V]';
%c.Limits=[0 max(M.N(:))];
if(isboolean(fixed) && fixed)
climits=caxis(plotaxes(2));
Minwell=climits(1);
Maxwell=climits(2);
elseif(~isboolean(fixed))
climits=caxis(plotaxes(2));
Minwell=climits(1);
Maxwell=climits(2);
end
view(plotaxes(2),2)
colormap(plotaxes(2),'jet');
linkaxes(plotaxes);
end
function plotGridButtonPushed(btn,ax)
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;
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,double((fieldstep-1)*M.it1)*M.dt))
%% update Position histogram
dens=M.N(:,:,fieldstep);
well=-PotentialWell(M.pot(:,:,fieldstep),M.Athet(:,:),M.rgrid,M.zgrid,false);
well=well(2:end,2:end);
plotaxes(1).Children.ZData=dens;
plotaxes(1).Children.CData=dens;
if(isboolean(fixed) && fixed)
climits=caxis(plotaxes(1));
MaxN=climits(2);
nmax=max(dens(:));
if(nmax>MaxN)
MaxN=nmax;
end
caxis(plotaxes(1),[0 MaxN]);
elseif(~isboolean(fixed))
MaxN=fixed.data;
caxis(plotaxes(1),[-Inf MaxN]);
end
plotaxes(2).Children.ZData=well;
plotaxes(2).Children.CData=well;
if(isboolean(fixed) && fixed || ~isboolean(fixed))
climits=caxis(plotaxes(2));
Maxwell=max([well(:);climits(2)]);
Minwell=min([well(:);climits(1)]);
caxis(plotaxes(2),[Minwell Maxwell]);
end
drawnow limitrate
end
end
c4science · Help