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

fennecs/matlab/@fennecshdf5/displaygeometry.m
762e842ba17cmaster

Last Change
Raw File

displaygeometry.m
View Options

function displaygeometry(obj,ax,xlimits,ylimits)
%DISPLAYGEOMETRY Display the simulation geometry in the axis ax
% if ax is not specified, a new figure is generated\
savefile=false;
if nargin<2
fig=figure('Name', sprintf('%s geometry',obj.name));
ax=axes(fig);
savefile=true;
end
geomw=obj.geomweight(:,:,1);
%% Magnetic field lines
Blines=obj.rAthet;
levels=linspace(min(Blines(obj.geomweight(:,:,1)>0)),max(Blines(obj.geomweight(:,:,1)>0)),60);
[~,h1]=contour(ax,obj.zgrid*1000,obj.rgrid*1000,Blines,real(levels),'-.','color','k','linewidth',1.5,'Displayname','Magnetic field lines');
hold on
%% Equipotential lines
Pot=obj.potxt(:,:,1);
Pot(obj.geomweight(:,:,1)<0)=NaN;
potcolor='b';
[c1,h2]=contour(ax,obj.zgrid*1000,obj.rgrid*1000,Pot/1e3,'--','color',potcolor,'linewidth',1,'Displayname','Equipotentials [kV]');
%clabel(c1,h2,'Color',potcolor)
%clabel(c1,h2, 'labelspacing', 200);
% % Grey outline shows metallic parts
% if obj.spl_bound.nbsplines<1
% [c]=contourc(obj.zgrid*1000,obj.rgrid*1000,geomw, [0 0]);
% [x,y,z]=C2xyz(c);
%
% else
% x=cell(obj.spl_bound.nbsplines,1);
% y=cell(obj.spl_bound.nbsplines,1);
% for i=1:obj.spl_bound.nbsplines
% z=obj.spl_bound.boundary(i).coefs(:,1)*1000;
% r=obj.spl_bound.boundary(i).coefs(:,2)*1000;
% if z(1)>obj.zgrid(end)*1000 ||z(end)<obj.zgrid(1)*1000
% z=[z(1); z];
% r=[obj.rgrid(end)*1000; r];
% end
% if z(1)<obj.zgrid(1)*1000 || z(end)>obj.zgrid(end)*1000
% z=[z(1); z];
% r=[obj.rgrid(1)*1000; r];
% end
% if r(1)>obj.rgrid(end)*1000 || r(end)<obj.rgrid(1)*1000
% r=[r(1); r];
% z=[obj.zgrid(end)*1000; z];
% end
% if r(1)<obj.rgrid(1)*1000 || r(end)>obj.rgrid(end)*1000
% r=[r(1); r];
% z=[obj.zgrid(1)*1000; z];
% end
% if r(1)>obj.rgrid(end)*1000
% end
% x{i}=z;
% y{i}=r;
% end
% end
% for i=1:length(x)
% z=[x{i}];
% r=[y{i}];
% patch(z,r,[150,150,150]./255,'edgecolor','none')
% end
% set the axia limits
xlim(ax,[obj.zgrid(1)*1000 obj.zgrid(end)*1000])
if(obj.conformgeom)
ylim(ax ,[obj.rgrid(1)*1000 obj.rgrid(rgridend)*1000])
else
ylim(ax,[obj.rgrid(1)*1000 obj.rgrid(end)*1000])
end
[~, name, ~] = fileparts(obj.file);
% 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>3
xlim(xlimits)
end
if nargin>4
ylim(ylimits)
end
xlabel(ax,'z [mm]')
ylabel(ax,'r [mm]')
% Plot of the geometry and of the different boundaries
% considered for the calculation of the currents
geomw=obj.geomweight(:,:,1);
geomw(geomw>1e-6)=NaN;
geomw(geomw<1e-6)=max(Pot(:)/1e3);
[c1,hContour]=contourf(ax,obj.zgrid*1000,obj.rgrid*1000,geomw, max(Pot(:)/1e3)*[1 1]);
hold on
dirw=obj.dirichletweight(:,:,1);
mask=obj.geomweight(:,:,1)<=1e-6&dirw>=1e-6;
dirw(mask)=min(Pot(:)/1e3);
dirw(~mask)=NaN;
[c1,hContour2]=contourf(ax,obj.zgrid*1000,obj.rgrid*1000,dirw,min(Pot(:)/1e3)*[1 1]);
grid on;
%legend([h1,h2],{'Magnetic field lines','Equipotentials [kV]'},'location','northeast','Orientation','horizontal')
colormap(ax,([255,165,0;255,255,255;150,150,150])./255)
% We add the coloring for the solid parts of the geometry
% Metallic parts
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
drawnow;
% Dielectrics
hFills=hContour2.FacePrims;
[hFills.ColorType] = deal('truecoloralpha'); % default = 'truecolor'
try
hFills(1).ColorData = uint8([255;165;0;255]);
for idx = 2 : numel(hFills)
hFills(idx).ColorData(4) = 0; % default=255
end
catch
end
drawnow;
axis equal
if savefile
obj.savegraph(fig,sprintf('%s/%sgeometry',obj.folder,obj.name),[16 14]);
end
end
c4science · Help