Homec4science
Diffusion nonlinear diocotron solver (master) Browse

nonlindioc/matlab/plotfields.m
0b252c046024master

Last Change
Raw File

plotfields.m
View Options

resfile='/misc/nonlinear/results_testl6_high.h5';
n=5501;
M=nonlinhdf5(resfile);
t=(n-1)*M.dt;
r=M.rgrid;
thet=M.thetgrid;
thet=thet*2*pi/max(thet);
[R,THET]=meshgrid(r,thet);
X=R.*cos(THET);
Y=R.*sin(THET);
O=[0 0];
rint=M.ra;
rext=M.rb;
% External fields
phiext=M.phiext';
Erext=M.Erext';
Ethetext=M.Ethetext';
% Self-consistent fields
phi=M.phi(:,:,n)';
Er=M.Er(:,:,n)';
Ethet=M.Ethet(:,:,n)';
% Total fields
phitot=phiext+phi;
Ertot=Erext+Er;
Ethettot=Ethetext+Ethet;
if 0
%% External fields
figure('Position', [100, 100, 1600, 500])
subplot(1,3,1);
pcolor(X,Y,phiext)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['\phi_{ext} [V] (t = ', num2str(t),' s)']);
subplot(1,3,2);
pcolor(X,Y,Erext)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['E_{r,ext} [V/m] (t = ', num2str(t),' s)']);
subplot(1,3,3);
pcolor(X,Y,Ethetext)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['E_{\theta,ext} [V/m] (t = ', num2str(t),' s)']);
suptitle('External fields')
end
if 0
%% Self-consistent fields
figure('Position', [100, 100, 1600, 500])
subplot(1,3,1);
pcolor(X,Y,phi)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['\phi [V] (t = ', num2str(t),' s)']);
subplot(1,3,2);
pcolor(X,Y,Er)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['E_r [V/m] (t = ', num2str(t),' s)']);
subplot(1,3,3);
pcolor(X,Y,Ethet)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['E_{\theta} [V/m] (t = ', num2str(t),' s)']);
suptitle('Self-consistent fields')
end
if 1
%% Total fields
figure('Position', [100, 100, 1600, 500])
subplot(1,3,1);
pcolor(X,Y,phitot)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['\phi [V] (t = ', num2str(t),' s)']);
subplot(1,3,2);
pcolor(X,Y,Ertot)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['E_r [V/m] (t = ', num2str(t),' s)']);
subplot(1,3,3);
pcolor(X,Y,Ethettot)
shading interp
colorbar
rectangle('Position',[O(1)-rint O(2)-rint 2*rint 2*rint],'Curvature',[1 1])
hold on
rectangle('Position',[O(1)-rext O(2)-rext 2*rext 2*rext],'Curvature',[1 1])
axis square
xlabel('x [m]')
ylabel('y [m]')
title(['E_{\theta} [V/m] (t = ', num2str(t),' s)']);
suptitle('Total (external + self) fields')
end
c4science · Help