

clear; clc;

load('fig8.mat');

figure;
set(gca,'color',[1 1 1]);
set(gca,'defaultAxesFontSize',20);
box on; grid on;
set(gcf,'Position',[300 250 750 300],'Unit','pixels');
for i=1:length(fig8.meavlf(:,1))  % for all station
    subplot(2,4,i);
    if i==1 || i==4 || i==6 || i==8
        plot(fig8.meat,fig8.meavlf(i,:),'m','lineWidth',1); hold on;        
    else
        plot(fig8.meat,fig8.meavlf(i,:),'b','lineWidth',1); hold on;  
    end
    plot(fig8.modelt,fig8.modelvlf(i,:),'k','lineWidth',2);
    datetick('x','mm-dd')
    xlabel('Time (UT)');
    ylabel('Amplitude (dB)')
    set(gca,'XLim',[fig8.xind+13/24 fig8.xind+15/24]);
    set(gca,'XTick',[fig8.xind+13/24:1/24:fig8.xind+15/24]);
    set(gca,'XTickLabel',{'13:00' '14:00' '15:00'});
    set(gca,'xminortick','on')
    xlabel('Time (UT)'); ylabel('Amplitude (dB)');
    box on; grid on; grid minor;
    set(gcf,'Position',[300 250 750 300],'Unit','pixels');
end