function [T_mean,y_mean,y_std,y_n]=daily_mean_fun(T,y)
%
% Вычисление среднего за сутки значениЯ компоненты полЯ
%        [T_mean,y_mean,y_std,y_n]=daily_mean_fun(T,y)
%
% T - сериальная дата
% y - измеренные значениЯ (NaN - пропуски)

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

T1=fix(min(T)); T2=fix(max(T)); % начальная и конечная даты
T_mean=T1:T2; T_mean=T_mean(:); % равномерная сетка
Nd=length(T_mean);

y_mean=ones(Nd,1)*NaN; y_std=y_mean; y_n=y_mean;
for kd=1:Nd
    k=between_my(T,T_mean(kd),T_mean(kd)+(86400-0.1)/86400);
    if sum(k)>1,
        y_mean(kd)=nanmean(y(k));
        y_std(kd)=nanstd(y(k));
        y_n(kd)=sum(k);
    end
end
T_mean=T_mean+0.5;
