clc
clear
close all

for start_month = 1:12

    start_time = '2020-01-01 18:00:00';
    end_time = '2020-01-02 6:00:00';
    
    figure()

    for year = 2018:2021

        date_day_mid = [num2str(year) '-' num2str(start_month,'%02d')];
        load(['.\dataset\hourly_average\' num2str(year) '\' date_day_mid '.mat'])



          if year == 2018
              errorbar([pairs.day.time], [pairs.day.Ne_Mid], [pairs.day.Ne_err],'-o','color',[1 0 0],'Linewidth', 5, 'MarkerSize', 30, 'Capsize', 20);
          elseif year == 2019
              errorbar([pairs.day.time], [pairs.day.Ne_Mid], [pairs.day.Ne_err],'-o','color',[0 1 0],'Linewidth', 5, 'MarkerSize', 30, 'Capsize', 20);
          elseif year == 2020
              errorbar([pairs.day.time], [pairs.day.Ne_Mid], [pairs.day.Ne_err],'-o','color',[0 0 1],'Linewidth', 5, 'MarkerSize', 30, 'Capsize', 20);
          else
              errorbar([pairs.day.time], [pairs.day.Ne_Mid], [pairs.day.Ne_err],'-o','color',[1 0.7 0],'Linewidth', 5, 'MarkerSize', 30, 'Capsize', 20);
          end

        hold on

    end

    set(gcf,'OuterPosition', [-7,33,1936,1056])
    set(gca,'Position', [.20 .20 .55 .75]);
    set(gca,'XLim',[datenum(start_time) datenum(end_time)])
    set(gca, 'XTick', [datenum(start_time):4/24:datenum(end_time)]);
    set(gca,'XTickLabel',{'18:00' '22:00' '02:00' '06:00'});
   
    set(gca,'xminortick','on');
    ax = gca;
    ax.XAxis.MinorTickValues = datenum(start_time):1/24:datenum(end_time);
    xlabel('Local Time(Hour)')

    set(gca, 'YLim', [20 32]);
    set(gca, 'YTick', [20:5:30]);
    set(gca,'YTickLabel',{'85' '90' '95'});
    set(gca,'Linewidth',2);
    set(gca,'yminortick','on');
    bx = gca;
    bx.YAxis.MinorTickValues = 20:32;
    ylabel('Electron Density(el/cm^3)')

    set(gca,'TickLength', [0.05 0.0125]);
    box on

    grid on
    grid minor


    figure_FontSize = 40;
    set(findobj('FontSize',10),'FontSize',figure_FontSize);
    
    
end