MY=F3(:,1);
Mminy=min(MY);
MM=F3(:,2)+12*(MY-Mminy);
MT=F3(:,8);

MSumma=splitapply(@sum,MT,findgroups(MM));

F3=F2(CSum>1000,:);%F3=F2(CSum>1500,:);

S=histcounts(F3(:,1));

txtfiles = dir( '*.txt' );
files={txtfiles(~[txtfiles.isdir]).name};
ss=1:size(files,2);
Nss=ss(1,end);
Nss1=Nss/2;
N=1;
for Nn=1:N;
for fnc = 1:ss(1,Nn);
filenameP1=[num2str(txtfiles(fnc).name) ];
ncidP1=importdata(files{fnc});
%end
A=ncidP1.data;
%A=importdata('res/d0424469.txt');
data1=A;
F=data1(:,[1,2,3,6,7,8,10]);
%r=find(all(F<9999,2) & data1(:,4)>=0);
r=find(all(F<9999,2) & data1(:,4)>=0 & data1(:,2)>=4); % с 4 мес
F1=F(r,:);
F1(:,8)=F1(:,4).*(F1(:,4)-F1(:,7)).*(F1(:,5)<2.5);
Y=F1(:,1);
miny=min(Y);
M=F1(:,2)+12*(Y-miny);
T=F1(:,8);
Summa=splitapply(@sum,T,findgroups(M));
f=find(diff(M)~=0)+1;
K=T(f);
Summa(2:end)=Summa(2:end)+Summa(1:end-1).*(K~=0);
UM=unique(M);
C=fix(UM/12);
Month=mod(UM,12);
Year=miny+C;
mz=find(Month==0);
Year(mz)=Year(mz)-1;
Month(mz)=12;
ff=sort([1;f;find(T==0);length(T)+1]);
dff=diff(ff);
A=mat2cell(T,dff,1);
CSum=cell2mat(cellfun(@cumsum,A,'UniformOutput',false));
F2=[F1 CSum];
F3=F2(CSum>1000,:);
X=(F3(:,1));
Ygod=F3(:,1);
miny1=min(Ygod);maxy2=max(Ygod);
G=(maxy2-miny1)+1;
S2=histcounts(X,G);
%dlmwrite(' filenameP1', S2,'\t');
dlmwrite([num2str(filenameP1) 'G.txt'],S2);
end
end
 
 
Az=Ygod;
d=diff(Az)~=0;
cs=cumsum([0 d']);
g=groupcounts(cs');% вектор из 59 значений, за каждый год 
god1=miny;
% для объяснения алгоритма
   g1=g(1);% первый год получили количество значений 
   g2=g(2);% второй год общее количество значений
   g3=g(3);
   g4=g(4);
 
% нарезка интервалов месяцев
   XX1=(F3(1:g1,2));
   XX2=(F3(g1+1:g1+g2,2));
% нарезка интервалов лет
 
GOD1=(F3(1:g1,1));
   GOD2=(F3(g1+1:g1+g2,1));
% получили годы все
 
 
 AzXX1=XX1; % за каждый год получить сколько есть значений в  каждом месяце и уникальные месяцы, но если один месяц выпадает, то значение уже не совпадает с уникальным месяцем
 dXX1=diff(AzXX1)~=0;
  csXX1=cumsum([0 dXX1']);
   gXX1=groupcounts(csXX1');
    Cuniq1 = unique(XX1);
     C1 = cat(1,Cuniq1,gXX1);
     
     
     AzXX2=XX2;
dXX2=diff(AzXX2)~=0;
 csXX2=cumsum([0 dXX2']);
  gXX2=groupcounts(csXX2');
   Cuniq2 = unique(XX2);
    C2 = cat(1,Cuniq2,gXX2);
 CunikALL = cat(1,Cuniq1    ,Cuniq2 ,Cuniq3 ,Cuniq4 ,Cuniq5 ,Cuniq6 ,Cuniq7 ,Cuniq8 ,Cuniq9 ,Cuniq10    ,Cuniq11    ,Cuniq12    ,Cuniq13    ,Cuniq14    ,Cuniq15    ,Cuniq16    ,Cuniq17    ,Cuniq18    ,Cuniq19    ,Cuniq20    ,Cuniq21    ,Cuniq22    ,Cuniq23    ,Cuniq24    ,Cuniq25    ,Cuniq26    ,Cuniq27    ,Cuniq28    ,Cuniq29    ,Cuniq30    ,Cuniq31    ,Cuniq32    ,Cuniq33    ,Cuniq34    ,Cuniq35    ,Cuniq36    ,Cuniq37    ,Cuniq38    ,Cuniq39    ,Cuniq40    ,Cuniq41    ,Cuniq42    ,Cuniq43    ,Cuniq44    ,Cuniq45    ,Cuniq46    ,Cuniq47    ,Cuniq48    ,Cuniq49    ,Cuniq50    ,Cuniq51    ,Cuniq52    ,Cuniq53    ,Cuniq54    ,Cuniq55    ,Cuniq56    ,Cuniq57    ,Cuniq58);
 gXXALL = cat(1,gXX1    ,gXX2   ,gXX3   ,gXX4   ,gXX5   ,gXX6   ,gXX7   ,gXX8   ,gXX9   ,gXX10  ,gXX11  ,gXX12  ,gXX13  ,gXX14  ,gXX15  ,gXX16  ,gXX17  ,gXX18  ,gXX19  ,gXX20  ,gXX21  ,gXX22  ,gXX23  ,gXX24  ,gXX25  ,gXX26  ,gXX27  ,gXX28  ,gXX29  ,gXX30  ,gXX31  ,gXX32  ,gXX33  ,gXX34  ,gXX35  ,gXX36  ,gXX37  ,gXX38  ,gXX39  ,gXX40  ,gXX41  ,gXX42  ,gXX43  ,gXX44  ,gXX45  ,gXX46  ,gXX47  ,gXX48  ,gXX49  ,gXX50  ,gXX51  ,gXX52  ,gXX53  ,gXX54  ,gXX55  ,gXX56  ,gXX57  ,gXX58);
GUnic=cat(1,CunikALL,gXXALL);
k4 = find(A111==4);
AN=GUnic(:,2);
% также и с позицией по году
GOD1=(F3(1:g1,1));
   GOD2=(F3(g1+1:g1+g2,1));
AzGOD1=GOD1;
     dGOD1=diff(AzGOD1)~=0;
     csGOD1=cumsum([0 dGOD1']);
     gGOD1=groupcounts(csGOD1');
     CuniqG = unique(GOD1);
     CGOD1 = cat(2,CuniqG,gGOD1);
  AzGOD2=GOD2;
  dGOD2=diff(AzGOD2)~=0;
  csGOD2=cumsum([0 dGOD2']);
  gGOD2=groupcounts(csGOD2');
  CuniqG2 = unique(GOD2);
  C2 = cat(2,CuniqG2,gGOD2);
 
  CunikALLGOD = cat(1,CuniqG    ,CuniqG2    ,CuniqG3    ,CuniqG4    ,CuniqG5    ,CuniqG6    ,CuniqG7    ,CuniqG8    ,CuniqG9    ,CuniqG10   ,CuniqG11   ,CuniqG12   ,CuniqG13   ,CuniqG14   ,CuniqG15   ,CuniqG16   ,CuniqG17   ,CuniqG18   ,CuniqG19   ,CuniqG20   ,CuniqG21   ,CuniqG22   ,CuniqG23   ,CuniqG24   ,CuniqG25   ,CuniqG26   ,CuniqG27   ,CuniqG28   ,CuniqG29   ,CuniqG30   ,CuniqG31   ,CuniqG32   ,CuniqG33   ,CuniqG34   ,CuniqG35   ,CuniqG36   ,CuniqG37   ,CuniqG38   ,CuniqG39   ,CuniqG40   ,CuniqG41   ,CuniqG42   ,CuniqG43   ,CuniqG44   ,CuniqG45   ,CuniqG46   ,CuniqG47   ,CuniqG48   ,CuniqG49   ,CuniqG50   ,CuniqG51   ,CuniqG52   ,CuniqG53   ,CuniqG54   ,CuniqG55   ,CuniqG56   ,CuniqG57   ,CuniqG58);