function mask_NaN=remove_radio_generation_fun(tau,dt_session,dt_min,dt_max)
% генерация маски из NaN для удаления помех при ЛЧМ-зондировании
%
%        mask_NaN=remove_radio_generation_fun(tau,dt_session,dt_min,dt_max);
%
% tau - периодичность измерений
% dt_session - периодичность сеансов, в сек.
% dt_min - начало чистки (относительно 00 сек) в сек
% dt_max - конец чистки (относительно 00 сек) в сек
% mask_NaN - вектор из 0, где NaN заполнены данные во время ЛЧМ-сеансов
%
% 1) повторяемость сеансов - st_session (минут)
% 2) NaN ставится на dt_min...dt_max сек от начала каждого сеанса
% 3) учет маски - сложением с данными измерений

% 06.05.2013 - С.Хомутов, Паратунка
% 09.05.2013 - параметры чистки задаются как параметры ф-ции при вызове
% 22.03.2014 - сделано dt_nan(dt_nan<0)=dt_nan(dt_nan<0)+dt_session;
% 05.08.2020 - секунды интервала зондирования округляются
%               dt_nan=fix((dt_min:dt_max)');

% if dt_min<0,
%     dt_nan=[0:dt_max dt_session+dt_min:dt_session-1]'; % секунды ЛЧМ-зондирования
% else
%     dt_nan=(dt_min:dt_max)'; % секунды ЛЧМ-зондирования
% end

dt_nan=fix((dt_min:dt_max)'); % секунды ЛЧМ-зондирования
% if any(dt_nan<=0),
%     dt_nan=dt_nan+dt_session;
% end
dt_nan(dt_nan<0)=dt_nan(dt_nan<0)+dt_session;

ts_mask=(0:86400/tau-1)'*tau; % равномерное время за сутки
mask_NaN=zeros(length(ts_mask),1);
% for ksec=0:24,
for ksec=1:length(dt_nan),
    k=find(rem(fix(ts_mask),dt_session)==dt_nan(ksec));
    mask_NaN(k)=NaN;
end
