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

fennecs/matlab/@fennecshdf5/displayNErprof.m
4ebefa948a9cmaster

Last Change
Raw File

displayNErprof.m
View Options

function displayNErprof(obj,zpos,fieldstep)
% Displays the electric field and density in 1D at the axial position zpos
% and time obj.t2d(timestep)
% The vessel radial boundaries are shown by vertical dashed lines
f=figure;
geomw=obj.geomweight(:,zpos,1);
geomw(geomw>=0)=1;
geomw(geomw<0)=NaN;
ax1=gca;
p=plot(ax1,obj.rgrid*1e3,obj.N(:,zpos,fieldstep),'linewidth',1.5);
xlim(ax1,[obj.rgrid(1) obj.rgrid(end)]*1e3)
xlabel(ax1,'r [mm]')
%title(ax1,'Density')
ylabel(ax1,'n_e[m^{-3}]');
%c.Limits=[0 max(M.N(:))];
hold(ax1, 'on')
M=obj;
[~,rcenterid]=max(M.geomweight(:,zpos,1));
[~,id1]=min(abs(M.geomweight(1:rcenterid,zpos,1)));
[~,id2]=min(abs(M.geomweight(rcenterid:end,zpos,1)));
id2=id2+rcenterid;
rlim1=M.rgrid(id1)*[1 1]*1e3;
rlim2=M.rgrid(id2)*[1 1]*1e3;
ylimits=ylim;
plot(ax1,rlim1,ylimits,'k--','linewidth',1.5,'Displayname','Boundaries');
plot(ax1,rlim2,ylimits,'k--','linewidth',1.5,'Displayname','Boundaries');
yyaxis(ax1,'right')
hold(ax1, 'on')
Er=obj.Er(:,zpos,fieldstep).*geomw;
Ez=obj.Ez(:,zpos,fieldstep).*geomw;
Erxt=obj.Erxt(:,zpos,1).*geomw;
Ezxt=obj.Ezxt(:,zpos,1).*geomw;
p1=plot(ax1,obj.rgrid*1e3,Er,'linewidth',1.5);
p2=plot(ax1,obj.rgrid*1e3,Ez,'linewidth',1.5);
p3=plot(ax1,obj.rgrid*1e3,Erxt,'linewidth',1.5);
p4=plot(ax1,obj.rgrid*1e3,Ezxt,'linewidth',1.5);
ylabel(ax1,'E [V/m]')
if max(abs([Er(:); Ez(:)]))>0
ylim(ax1,[ -max(abs([Er(:); Ez(:)])) max(abs([Er(:); Ez(:)]))])
end
legend(ax1,[p p1 p2 p3 p4],{'n_e','E_r','E_z','Vacuum E_r','Vacuum E_z'},'location','northwest')
xlim([M.rgrid(id1) M.rgrid(id2)]*1e3)
obj.savegraph(f,sprintf('%s/%s_NEproft%i_z%i',obj.folder,obj.name,fieldstep,zpos));
end
c4science · Help