%%Figure 5 apparent frequency via time
%use example waveform from 300km
load averaged_night_time.mat
[ afg, afG,wafg, afpeak ] = calculate_apparent_freuqency( wfbank(:,30), 10.8e3, 2e3,151 );
figure
subplot(2,1,2)
plot(-1000/1000:1/1000:3999/1000,(afg+10.8e3)/1000,'--','color','k','linewidth',2)
hold on
plot(-1000/1000:1/1000:3999/1000,(wafg+10.8e3)/1000,'color','r','linewidth',2)
plot(-1000/1000:1/1000:3999/1000,abs(gtrace(gfilter(wfbank(:,30),1e6,8.8e3,12.8e3)))*1.5e6/1000,'b','linewidth',2)
xlabel('[t] = ms','FontSize', 20);
ylabel('[fi] = kHz','FontSize', 20);
set(gca,'fontsize',16)
ylim([-5 25])
xlim([-1 2])
grid


% load averaged_night_time.mat
[ afg, afG,wafg, afpeak ] = calculate_apparent_freuqency( wfbank(:,30), 10e3, 8e3,151 );
% figure
subplot(2,1,1)
plot(-1000/1000:1/1000:3999/1000,(afg+10.8e3)/1000,'--','color','k','linewidth',2)
hold on
plot(-1000/1000:1/1000:3999/1000,(wafg+10.8e3)/1000,'color','r','linewidth',2)
plot(-1000/1000:1/1000:3999/1000,abs(gtrace(gfilter(wfbank(:,30),1e6,2e3,18e3)))*5e5/1000,'b','linewidth',2)
% xlabel('[t] = \mus','FontSize', 20);
ylabel('[fi] = kHz','FontSize', 20);
set(gca,'fontsize',16)
ylim([-10 25])
xlim([-1 2])
grid
print fig5.1.tif -dtiff