gaussscript.m

zi=-10;
xi=-10;
zf=10;
xf=10;
zstepfrac=0.1;
xstepfrac=0.1;
size=int32(((zf-zi)/(zstepfrac)+1));
eta = 3;
eps2=-1;
mu2=-1;
thetai='0.16pi';

 
g=[3.0]; 
for x=1:1
g(x) = round(g(x)*10^1)/(10^1);
end
sizesquare=int32(double(size)*double(size));
data=zeros(sizesquare, 4, 7);

for x=1:1
num2str(g(x), '%3.1f');
data(:,:,x) = load(strcat('data/',thetai,'rads',num2str(eta, '%3.1f'),'eta',num2str(g(x), '%3.1f'),'g',num2str(eps2, '%4.1f'),'eps2gaussdielecfieldmap.dat'));
end

for x=1:1

i=1;
j=1;
jcount=0;



eyarray=zeros(size,size);%this is the tranformed field modulus
eyrparray=zeros(size,size);%this is the transformed field real part
xarray=[xi:(xstepfrac):xf];
zarray=[zi:(zstepfrac):zf];


while (i <= size)
	while (j<=size)
		eyarray(i,j)=data((jcount* double(size) + j),3,x);
        eyrparray(i,j)=data((jcount* double(size) + j),4,x);
		j=j+1;
	end
	 jcount=jcount+1;
	 i=i+1;
	 j=1;
end

% 
% h=surf(xarray,zarray,eyarray);
% xlabel('x');
% ylabel('z');
% zlabel('Ey');
% title('thetai=0, d=1, mu1=1, mu2=1, eps1=1, eps2=10, eta=1');
% colorbar;
% set(h, 'edgecolor','none')
% % 
% %set(gca,'DataAspectRatio',[max(xarray) max(zarray) max(max(eyarray))]);
% %set(gca,'PlotBoxAspectRatio',[max(xarray) max(zarray) max(max(eyarray))]);
% 
% 
% print('-dpng', '0inc_eta1_eps11_eps210.png');

k=imagesc(xarray,zarray,eyarray);
xlabel('x/lambda');
ylabel('z/lambda');
%zlabel('Ey');
title(strcat('g=',num2str(g(1)),',thetai=',num2str(thetai),',mu1=1, mu2=',num2str(mu2),', eps1=1, eps2=',num2str(eps2),', eta=',num2str(eta) ));
colorbar;
line([xi xf],[eta eta],'linewidth',4,'Color', 'k');
line([xi xf],[0 0],'linewidth',4,'Color', 'k');
print('-dpng', strcat('plots/im_thetai',thetai,'_g',num2str(g(1)),'_eta',num2str(eta),'_eps11_eps2',num2str(eps2),'mu2=',num2str(mu2),'.png'));


k=imagesc(xarray,zarray,eyrparray);
xlabel('x/lambda');
ylabel('z/lambda');
%zlabel('Ey');
title(strcat('real part: g=',num2str(g(1)),',thetai=',num2str(thetai),',mu1=1, mu2=',num2str(mu2),', eps1=1, eps2=',num2str(eps2),', eta=',num2str(eta) ));
colorbar;
line([xi xf],[eta eta],'linewidth',4,'Color', 'k');
line([xi xf],[0 0],'linewidth',4,'Color', 'k');
print('-dpng', strcat('plots/real_im_thetai',thetai,'_g',num2str(g(1)),'_eta',num2str(eta),'_eps11_eps2',num2str(eps2),'mu2=',num2str(mu2),'.png'));
end