function C=remove_radio_noise_fun(T,C,tau);
% заполнение символами NaN помех при радиозондировании в измерениях FGE
% каждые 5 минут, начало в 00 сек, конец в 56 сек
% частота измерений 2 Гц (0.5 сек)
%
%     C=remove_radio_noise_fun(T,C);
%
% T   - сериальная дата измерений
% C   - измерения
% tau - интервал между измерениями FGE
%
% 1) ряд данных - не более суток
% 2) первое значение времени используется для генерирования опорного
%   времени
% 3) маска в виде 1 и NaN для ряда 86400*2 значений - в файле Remove_Noise

% 28.04.2013 - С.Хомутов, Паратунка
% 05.05.2013 - маска NaN-ов генерируется с помощью
%               remove_radio_generation_fun(), а не из файла
% 02.02.2018 - маска с NaN может генерироваться только для части суток
% 11.01.2020 - при вычисление ts_obs rem(T,1) заменено на T-nanmean(T)
% 05.08.2020 - значение продолжительности сессии в сек округляется до
%               целого - dt_session=round(tmp(:,13)*60);

% tau=0.5;    % периодичность измерений FGE

% инициализация постоянных длЯ обработки
run ..\Ini_parameters\fge_ini.m

file_session=[path_all 'ini_parameters\Radio_Sessions.dat'];
tmp=load(file_session);
T1_LFM=datenum(tmp(:,1:6)); T2_LFM=datenum(tmp(:,7:12));
dt_session=round(tmp(:,13)*60);
dt_min=tmp(:,14); dt_max=tmp(:,15);

ts_mask=(0:86400/tau-1)'*tau; % равномерное время за сутки
T0=fix(nanmean(T)); % текущая дата
ts_obs=round((T-T0)*86400*10)/10; % время измерений в сек

for k=1:length(dt_session), % по всем параметрам чистки
    if any(between(T,T1_LFM(k),T2_LFM(k))), % если попадает в интервал
        % генерация маски с 0 и NaN
        disp('... Чистка помех при ЛЧМ-зондировании ...')
        mask_NaN=remove_radio_generation_fun(tau,...
                             dt_session(k),dt_min(k),dt_max(k));
        % оставляем маску с NaN только для отмеченной части суток (02.02.2018)
        k1=between_my(ts_mask,(T1_LFM(k)-T0)*86400,(T2_LFM(k)-T0)*86400);
        mask=zeros(length(mask_NaN),1); mask(k1)=mask_NaN(k1);
        mask_NaN=mask;
        [kmask,kobs]=ismember_digit(ts_mask,ts_obs,9);
        C(kobs)=C(kobs)+mask_NaN(kmask);
    end
end
