function [mjd_mean,y_mean,y_std,y_n]=one_day_mean(mjd,y);
%
% Вычисление среднего за сутки значениЯ компоненты полЯ
%        [mjd_mean,y_mean,y_std,y_n]=one_day_mean(mjd,y);
%
% mjd - модифицированнаЯ юлианскаЯ дата
%   y - измеренные значениЯ (NaN - пропуски)

% 03.05.2006 - с учетом пропусков, заполненными как NaN

kd=find(diff(fix(mjd))>0);    % ищем индексы смены суток
kd=kd(:); kd=[kd;length(mjd)];% добавлЯем последние сутки
n=length(kd);            % кол-во суток с измерениЯми

fprintf('Всего суток длЯ обработки ... %d\n', n);

mjd_mean=zeros(n,1); y_mean=zeros(n,1); y_std=zeros(n,1); y_n=zeros(n,1);
for k=1:n,
    if rem(k,50)==0, fprintf(' %5d',k); end
    if rem(k,500)==0 | k==n, fprintf('\n'); end
    kk=fix(mjd)==fix(mjd(kd(k)));
    mjd_mean(k)=nanmean(mjd(kk));   % mjd длЯ вычислениЯ статистики
    y_mean(k)=nanmean(y(kk));       % среднее значение с пропусками NaN
    y_std(k)=nanstd(y(kk));
    y_n(k)=sum(~isnan(y(kk)));
end
