Homec4science
Diffusion FENNECS lite (3Dmpimode) Browse

fennecslite/matlab/@fennecslitehdf5/displayconfiguration.m
6b37315661933Dmpimode

Last Change
Raw File

displayconfiguration.m
View Options

function f=displayconfiguration(obj,fieldsteps,xlimits,ylimits)
%displayconfiguration plot the configuration of the simulation
% domain withe boundaries, the magnetic field lines the
% electric equipotential lines and the electron density
% averaged in time between t2d(fieldsteps(1)) and t2d(fieldsteps(end))
fieldstart=fieldsteps(1);
fieldend=fieldsteps(end);
if nargin<3
xlimits=inf*[-1 1];
end
if nargin<4
ylimits=inf*[-1 1];
end
dens=mean(obj.N(:,:,fieldstart:fieldend),3);
geomw=obj.geomweight(:,:,1);
maxdens=max(dens(:));
geomw(geomw<0)=0;
geomw(geomw>0)=maxdens;
dens(geomw<=0)=0;
geomw(geomw>0)=NaN;
f=figure('Name', sprintf('%s fields',obj.name));
ax1=gca;
title(sprintf('Configuration'))
%dens(dens<=1e13)=NaN;
%% electron density
h=contourf(ax1,obj.zgrid*1000,obj.rgrid*1000,dens,50,'Displayname','n_e [m^{-3}]', 'linestyle','none');
hold on;
colormap(flipud(hot));
%% Magnetic field lines
Blines=obj.rAthet;
levels=linspace(min(Blines(obj.geomweight(:,:,1)>0)),max(Blines(obj.geomweight(:,:,1)>0)),15);
[~,h1]=contour(ax1,obj.zgrid*1000,obj.rgrid*1000,Blines,real(levels),'-.','color','k','linewidth',1.1,'Displayname','Magnetic field lines');
%% Equipotential lines
Pot=mean(obj.pot(:,:,fieldstart:fieldend),3);
Pot(obj.geomweight(:,:,1)<0)=NaN;
%levels=8;%[-3.4 -5 -10 -15 -20 -25];%7;
potcolor='b';%[0.3660 0.6740 0.1880];
[c1,h2]=contour(ax1,obj.zgrid*1000,obj.rgrid*1000,Pot/1e3,[-18:4:-2],'--','color',potcolor,'ShowText','on','linewidth',1.2,'Displayname','Equipotentials [kV]');
clabel(c1,h2,'Color',potcolor)
clabel(c1,h2, 'labelspacing', 450);
clabel(c1,h2, 'fontsize', 13);
% Grey outline shows metallic parts
[c1,hContour]=contourf(ax1,obj.zgrid*1000,obj.rgrid*1000,geomw, [0 0]);
drawnow;
% set the axia limits
xlim(ax1,[obj.zgrid(1)*1000 obj.zgrid(end)*1000])
if(obj.conformgeom)
ylim([ax1 ],[obj.rgrid(1)*1000 obj.rgrid(rgridend)*1000])
else
ylim([ax1],[obj.rgrid(1)*1000 obj.rgrid(end)*1000])
end
legend([h1,h2],{'Magnetic field lines','Equipotentials [kV]'},'location','northeast')
xlabel(ax1,'z [mm]')
ylabel(ax1,'r [mm]')
c = colorbar(ax1);
c.Label.String= 'Electron density [m^{-3}]';
view(ax1,2)
set(ax1,'fontsize',14)
grid on;
hFills=hContour.FacePrims;
[hFills.ColorType] = deal('truecoloralpha'); % default = 'truecolor'
try
hFills(1).ColorData = uint8([150;150;150;255]);
for idx = 2 : numel(hFills)
hFills(idx).ColorData(4) = 0; % default=255
end
catch
end
[~, name, ~] = fileparts(obj.file);
% with this you could show the outline of the maxwellian source
% if obj.maxwellsrce.present
% rlen=diff(obj.maxwellsrce.rlim);
% zlen=diff(obj.maxwellsrce.zlim);
% rectangle('Position',[obj.maxwellsrce.zlim(1) obj.maxwellsrce.rlim(1) zlen rlen]*1000,'Edgecolor','g','Linewidth',2,'Linestyle','--')
% end
% in case of coaxial configuration, extend the display domain
% and add grey rectangles to show metallic boundaries
if( obj.walltype >=2 && obj.walltype<=4)
rectangle('Position',[obj.zgrid(1) obj.r_b obj.zgrid(end)-obj.zgrid(1) 0.001]*1e3,'FaceColor',[150 150 150]/255,'Edgecolor','none')
ylimits=ylim;
ylim([ylimits(1),ylimits(2)+1])
end
if sum(obj.geomweight(:,1,1))==0
rectangle('Position',[obj.zgrid(1) obj.r_a-0.001 obj.zgrid(end)-obj.zgrid(1) 0.001]*1e3,'FaceColor',[150 150 150]/255,'Edgecolor','none')
ylimits=ylim;
ylim([ylimits(1)-1,ylimits(2)])
end
if nargin>2
xlim(xlimits)
end
if nargin>3
ylim(ylimits)
end
f.PaperUnits='centimeters';
%axis equal
papsize=[14 5 ];
%%
axis equal
%ylim([60 85])
%xlim([-40 80])
legend('location','northeast')
legend off
Pos=f.Position;
Pos(3)=1.8*Pos(3);
f.Position=Pos;
obj.savegraph(f,sprintf('%s/%sFields',obj.folder,obj.name),papsize);
end
c4science · Help