function [c_m,c_s,c_n]=one_hour_median(mjd0,mjd,c);
% Вычисление медианного значениЯ за 1 час по 1-минутным данным
%
% [c_m,c_s,c_n]=one_hour_median(mjd0,mjd,c)
%
% mjd0 - юлианскаЯ дата моментов, длЯ которых будет считатьсЯ медиана
%        (mjd0(k)-0.5/24 ... mjd0(k)+0.5/24)
%  mjd - юлианскаЯ дата моментов, длЯ которых имеютсЯ измерениЯ (1-минутные)
%    c - измеренные значениЯ
%
%   c_m - медиана за час (==Inf, если количество измерений <60)
%   c_s - стандартное отклонение
%   c_n - количество измерений длЯ данного часа (<60 длЯ Inf)
%

n=length(mjd0); % кол-во часовых значений
k=c~=Inf; mjd=mjd(k); c=c(k);   % исключаем пропуски в данных (=Inf)

c_m=zeros(n,1); c_s=zeros(n,1); c_n=zeros(n,1);
for k=1:n,
    if mod(k,100)==0, 
        fprintf('%7d',k), if mod(k,500)==0, fprintf('\n'), end, 
    end
    kk=between(mjd,mjd0(k)-0.5/24,mjd0(k)+0.5/24,'11');
    c_n(k)=sum(kk);     % кол-во измерений, принЯтых в расчет
    if sum(kk)>=59,     % берем в обработку только полный час
        c_m(k)=median(c(kk)); % медианна за час
        c_s(k)=std(c(kk)); % STD
    else
        c_m(k)=Inf; % среднее за час
        c_s(k)=Inf; % STD
    end        
end
    