% Вычисление базисных в цикле для magdas

% 10.11.2012 - для нового азимута Миры
%            - градиенты считываются из файла текущих dFgrad
% 18.02.2013 - для нового азимута (по набл. в феврале 2013 г.)
%            - используется вычисленные градиент dFgrad_obs
%            - учитывается зависимость H(dD)
%            - dFgrad выводится в файл
% 10.03.2013 - проба с учетом температурной зависимости (параметры - в FGE_ini.m)
% 02.07.2013 - вариант с коррекцией зависимости Z=f(H,D)
% 22.01.2020 - может использоваться F как GSM-90, так и POS-1
% 30.09.2022 - дополнено вызовом fge_ini.m для инициации параметров GSM-90

clear
run ..\..\ini_parameters\fge_ini.m
run ..\..\ini_parameters\pos1_ini.m
run ..\..\ini_parameters\magdas_ini.m

% оставляем реально используемый скалярный магнитометр
% scal_mag='GSM-90';
scal_mag='POS-1';

file_blv_all='blv_all_md_TR.dat';
file_blv='blv_md_TR.dat';
file_sens='blv_sensor_md_TR.dat';

T1=datenum(2025,09,13);  % начальная дата
T2=datenum(2025,09,13);  % конечная дата
Nday=T2-T1+1;           % кол-во суток для обработки

for kday=1:Nday
    T=T1-1+kday;        % текущая сериальная дата
    date0=strrep(datestr(T,24),'/','.');
    file=[path_data_DI date0(7:10) '\DIF_' strrep(datestr(T,26),'/','') '.dat'];
    if ~exist(file), continue, end  % нет файла, переходим к следующей
    disp(['------------------------------------------------------- ' date0])
    dif_hms_cycle
    pause(2)
    disp('-------------------------------------------------------')
    base_calc_magdas_tr_cycle
    pause(2)
end
