clear all
close all

%% load constants
c=2.998e8;

%% parameters
q3i=[200 100 50 20 5 1];
q1=1;
q2=100*q1;

for ii=1:length(q3i)
q3=q3i(ii);
a=1;

% parameters for contours
N=20;
Nx=150;
Nz=300;
x=linspace(-N*q1,N*q1,Nx);
y=x;
z=linspace(-60*q1,60*q1,Nz);
[Y,Z,X]=meshgrid(y,z,x);

[Etheta1] = FDb(q1,q2,q3,a,-81.2/2/c,X,Y,Z);
[Etheta2] = FDb(q1,q2,q3,a,0*q1/c,X,Y,Z); 
[Etheta3] = FDb(q1,q2,q3,a,81.2/2*q1/c,X,Y,Z);

Etheta1(Z>-20.1)=0;
Etheta2(Z>19.9)=0;
Etheta2(Z<-19.9)=0;
Etheta3(Z<20.1)=0;
    %% E-field - real pulse
v1=real(Etheta1); 
v1=v1/max(max(max(abs(v1))));
v2=real(Etheta2); 
v2=v2/max(max(max(abs(v2))));
v3=real(Etheta3); 
v3=v3/max(max(max(abs(v3))));

figure;box on;
vt=v1+v2+v3;
imagesc(y,z,vt(:,:,length(x)/2));axis equal;axis tight;
load Mymap;colormap(Mymap),colorbar;set(gca,'Clim',[-1 1]);
xlim([y(1) y(end)]);ylim([z(1) z(end)]);
end