%%Figure 2 complex waveform
%use example 300km
fso=exp(-1j*2*pi*1e4*(0:5000-1)'/1e6);
cy2_18=gfilter(wfbank(:,30).*fso,1e6,0,8e3);
figure
% scatter3(imag(cy2_18),-1000:3999,real(cy2_18),4,abs(cy2_18)*1000,'filled')
scatter3(imag(cy2_18)*1000,-1000/1000:1/1000:3999/1000,real(cy2_18)*1000,4,-1:0.001:3.999,'filled')
caxis([-0.5 1])
hold on
plot3(imag(cy2_18)*1000,-1000/1000:1/1000:3999/1000,real(cy2_18)*1000,'k')
view(60,25)
scatter3(-0.017*ones(5000,1)*1000,-1000/1000:1/1000:3999/1000,real(cy2_18)*1000,2,abs(cy2_18))
plot3(-0.017*ones(5000,1)*1000,-1000/1000:1/1000:3999/1000,real(cy2_18)*1000)
plot3(-0.017*ones(5000,1)*1000,-1000/1000:1/1000:3999/1000,abs(cy2_18)*1000,'r')
scatter3(imag(cy2_18)*1000,-1000/1000:1/1000:3999/1000,-0.013*ones(5000,1)*1000,2,abs(cy2_18))
plot3(imag(cy2_18)*1000,-1000/1000:1/1000:3999/1000,-0.013*ones(5000,1)*1000)
plot3(abs(cy2_18)*1000,-1000/1000:1/1000:3999/1000,-0.013*ones(5000,1)*1000,'r')

plot3(-0.017*ones(5000,1)*1000,-1000/1000:1/1000:3999/1000,0.023*ones(5000,1)*1000,'k')
plot3(-0.017*ones(5000,1)*1000,-1000/1000:1/1000:3999/1000,-0.013*ones(5000,1)*1000,'k--')
plot3(-0.017*1000:(0.023+0.017)/4999*1000:0.023*1000,4000/1000*ones(5000,1),0.023*ones(5000,1)*1000,'k')
plot3(-0.017*1000:(0.023+0.017)/4999*1000:0.023*1000,4000/1000*ones(5000,1),-0.013*ones(5000,1)*1000,'k--')
plot3(0.02295*ones(5000,1)*1000,4000/1000*ones(5000,1),-0.013*1000:(0.023+0.013)/4999*1000:0.023*1000,'k')
plot3(-0.017*ones(5000,1)*1000,4000/1000*ones(5000,1),-0.013*1000:(0.023+0.013)/4999*1000:0.023*1000,'k--')


colorbar;
h=colorbar;
set(h,'position',[.8 .33 .035 .34])
title(h,'[t]=ms','fontsize',16);
xlabel('[Ai]=\muV/m','FontSize', 20);ylabel('[t]=ms','FontSize', 20);zlabel('[Ar]=\muV/m','FontSize', 20);
set(get(gca, 'XLabel'), 'Rotation', -32);
set(get(gca, 'YLabel'), 'Rotation', 14);
set(gca,'fontsize',16)
ylim([-1 4])
xlim([-0.017*1000  0.023*1000])
zlim([-0.013*1000 0.023*1000])
  print fig2.tif -dtiff