%%Figure 3 waveform spectrum
specb=log10(abs(fft(wfbank)*2/5000)*1000);
h=figure;imagesc(specb(1:500,:)'); shading flat;set(gca,'ydir','normal');hold on
% colormap(flipud(cbrewer('div','RdYlBu',32)));
colormap hot;
caxis([-1.5 1.2]);
plot([50 50],[30 60],'--','color','b','linewidth',2)
plot([50 100],[30 30],'--','color','b','linewidth',2)
plot([100 100],[30 60],'--','color','b','linewidth',2)
plot([50 100],[60 60],'--','color','b','linewidth',2)
colorbar;
h=colorbar;
ylabel(h,'log10([A]=mV/m)','fontsize',20)
xlabel('[f] = kHz','FontSize', 20);ylabel('[d] = km','FontSize', 20);
set(gca,'XTickLabel',{'20','40','60','80','100'})
set(gca,'YTickLabel',{'100','200','300','400','500','600','700','800','900','1000'})
set(gca,'fontsize',16)
 print fig3.tif -dtiff