%%Figure 6 apparent frequency via distance compare with spectrum
%use example waveform from 300km-600km
spec=abs(fft(wfbank(:,30:60)));
[a,b]=max(spec(51:100,:));
fc=b/5*1000+1e4;
for ii=1:31
    [ afg, afG,wafg,afpeak ] = calculate_apparent_freuqency( wfbank(:,ii+29), fc(ii), 2e3 ,151);
% [ afg, afG,wafg,afpeak ] = calculate_apparent_freuqency( wfbank(:,ii+29), 15e3, 5e3 );
    af(ii)=afpeak+fc(ii);
% af(ii)=afpeak+15e3;
end
fftwf=log10(abs(fft(wfbank(:,26:65))*2/5000)*1000);
h=figure;imagesc(fftwf(40:110,:)'); shading flat;set(gca,'ydir','normal');hold on
colormap(flipud(cbrewer('div','RdYlBu',32)));caxis([-1.5 0.4])
colorbar;
h=colorbar;
ylabel(h,'log10([A]=mV/m)','fontsize',20)
plot((af-8e3)/200,5:35,'color','r','linewidth',4)
plot([10 10],[5 35],'--','color','b','linewidth',2)
plot([10 60],[5 5],'--','color','b','linewidth',2)
plot([60 60],[5 35],'--','color','b','linewidth',2)
plot([10 60],[35 35],'--','color','b','linewidth',2)
set(gca,'fontsize',16)
xlabel('[f] = kHz','FontSize', 20);ylabel('[d] = km','FontSize', 20);
set(gca,'xtick',[10 35 60])
set(gca,'XTickLabel',{'10','15','20'})
set(gca,'YTickLabel',{'300','350','400','450','500','550','600','650'})
  print fig6.tif -dtiff