提取ERA5再分析降水数据的年、月值，不需要考虑闰年和平年

slc83328317

提取ERA5再分析降水数据的年、月值，不需要考虑闰年和平年的代码。

% clc;clear   
% load('C:\Users\admin\Desktop\paper\ERA5 EPEs\data\日期参考.mat')  %日期参考文件，用来提取年月日数据,本次是提取总降水，不是极端（第一列是年，第二列是月，第三列是日，用来参考的）
% % load('C:\Users\admin\Desktop\90%极端降水日数据mat\天数\Allday_frequency.mat')
% % load('C:\Users\admin\Desktop\90%极端降水日数据mat\天数\Daytime_frequency.mat')
% load('C:\Users\admin\Desktop\90%极端降水日数据mat\天数\Nighttime_frequency.mat')
% % load('C:\Users\admin\Desktop\90%极端降水日数据mat\总量\Allday_amount.mat')
% % load('C:\Users\admin\Desktop\90%极端降水日数据mat\总量\Daytime_amount.mat')
% % load('C:\Users\admin\Desktop\90%极端降水日数据mat\总量\Nighttime_amount.mat')
% SLC=zeros(58,41,41);
%
% for m=1:41
%     for n=1:41
%
% for i=1961:2018
%
%
% D=Nighttime_frequency(:,m,n);
%     
%     
% A=find(Num(:,1)==i,1,'first'); %提取下界1961的第一个数
% B=find(Num(:,1)==i,1,'last');
% data=D(A:B,:,:);
% Date=Num(A:B,:,:); %
% M=find(Date(:,2)==12,1,'first');%
% N=find(Date(:,2)==12,1,'last');%
% New_data=data(M:N,:,:);%
%
% % data=sum(data);
% data1=sum(New_data);
% SLC(i-1960,m,n)=data1;
% disp(i)
% end
%     end
% end
%
% ERA_nighttime_fre_mon12=SLC;