function k_NaN=groups_of_nan(x,y)
% выделение фрагментов непрерывно идущих NaNs и запись индексов границ
% этих фрагментов
%
% k_NaN=groups_of_nan(x,y)
%
% x,y - время и анализируемый ряд

% 17.02.2020 - S.Y. Khomutov, Observatory "Paratunka", Kamchatka, Russia

k=isnan(y);
kk=find(abs(diff(k))==1);
if isnan(y(1))  % NaNs в начале 
    if isnan(y(length(y))) % NaNs в конце
        kkk=[1; kk; length(x)];
    else
        kkk=[1; kk];
    end
else            % норм. измерение в начале
    if isnan(y(length(y))) % NaNs в конце
        kkk=[kk; length(x)];
    else
        kkk=[kk];
    end
end
k_NaN=[]; 
for l=1:2:length(kkk/2)
    k_NaN=[k_NaN; kkk(l)+1 kkk(l+1)]; 
end