clc
clear
close all
for m = 2:13
  
    
    figure()
    
    for deal_year = 2018:2021
        
        
         if ~exist(['.\dataset\daily_average\' num2str(deal_year) '_' num2str(m-1) '.mat'],'file')
            continue
         end

        load(['.\dataset\daily_average\' num2str(deal_year) '_' num2str(m-1) '.mat'])

        x_time = datenum(char(pairs.tweek.UTCTIME));
        start_time = [num2str(deal_year) '-01-01 00:00:00'];
        end_time = [num2str(deal_year+1) '-01-01 00:00:00'];


        if deal_year == 2020
         
            year_month = [0 31 29 31 30 31 30 31 31 30 31 30 31];

        else

            year_month = [0 31 28 31 30 31 30 31 31 30 31 30 31];

        end
        
        
        yy = [pairs.month.Ne_Mid];
        err = [pairs.month.Ne_err];
            
          if deal_year == 2018
              errorbar([1:year_month(m)]-0.5, yy(1:year_month(m)), err(1:year_month(m)),'-o','color',[1 0 0],'Linewidth', 5, 'MarkerSize', 20, 'Capsize', 20);
          elseif deal_year == 2019
              errorbar([1:year_month(m)]-0.5, yy(1:year_month(m)), err(1:year_month(m)),'-o','color',[0 1 0],'Linewidth', 5, 'MarkerSize', 20, 'Capsize', 20);
          elseif deal_year == 2020
              errorbar([1:year_month(m)]-0.5, yy(1:year_month(m)), err(1:year_month(m)),'-o','color',[0 0 1],'Linewidth', 5, 'MarkerSize', 20, 'Capsize', 20);
          else
              errorbar([1:year_month(m)]-0.5, yy(1:year_month(m)), err(1:year_month(m)),'-o','color',[1 0.7 0],'Linewidth', 5, 'MarkerSize', 20, 'Capsize', 20);
          end

        hold on

    end

    % 字体大小
    figure_FontSize = 40;

    set(gcf,'OuterPosition', [-7,33,1936,1056])
    set(gca,'Position', [.20 .20 .55 .75]);

    set(gca,'XLim',[0 31]);
    set(gca,'XTick',[0:5:31]);
    set(gca,'XTickLabel',{0:5:30});
    set(gca,'xminortick','on');
    ax = gca;
    ax.XAxis.MinorTickValues = 0:31 ;
    xlabel('Local Time(Day)')

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



    set(findobj('FontSize',10),'FontSize',figure_FontSize);
    set(gca,'Linewidth',2);
    set(gca,'TickLength', [0.025 0.0125]);
    box on
    
    
    grid on 
    grid minor
    
   
end