% основные инициализационные параметры

global path_all

% добавляем основную рабочую папку для обработки Паратунки
% в конец путей MatLab (зависит от конкретного PC)
path(pathdef)   % восстанавливаем пути по умолчанию
% addpath 'D:\Paratunka\FGE_GSM' -end
addpath 'd:\Paratunka\FGE_GSM' -end

% всем путям предшествует глобальный каталаг (зависит от PC)
path_all='d:\Paratunka\';

% каталог с данными абсол. наблюдений
path_data_DI=[path_all 'Data\Abs_Data\'];  % далее - YYYY\DIF_YYYYMMDD.dat

% каталоги с исходными данными FGE и GSM 
path_data_fge=[path_all 'data\FGE_GSM\momwe\fge\'];  % вариометр FGE
path_data_gsm=[path_all 'data\FGE_GSM\momwe\gsm\'];  % протонный GSM в вариац. пав-не
% path_data_fabs=[path_all 'di\data\'];        % протонный в абсол. пав-не

% среднегодовые значения поля (по оценкам в октябре 2012 г.)
H_year=21670; Z_year=47000; F_year=51750; 
D_year=-(5+55/60); I_year=65+15/60; % в угл.мин.
% уровни для рисования графиков
FF=51000; HH=21000; ZZ=46000; DD=-6; II=65; XX=21000; YY=-2000;

% периодичность измерения FGE и GSM (секунд)
dt_fge=0.5;
dt_gsm=5;

% кол-во точек, которые берём в интервале без помехи от GSM
% и систем. поправки из-за влияния GSM-90 (по 11-12.08.2018)
% H,Z в нТл, D в угл.мин.
% исправленное = исходное + поправка
% np_fge=0; dHsel=+0.0; dDsel=+0.0; dZsel=0.0;      % длЯ теста - наводку от GSM-90 не чистим
% np_fge=0; dHsel=+0.06; dDsel=+0.0024; dZsel=+0.04; % все данные - c 02.08.2018 !!!
% np_fge=0; dHsel=+0.10; dDsel=+0.025; dZsel=+0.30; % все данные - c 29.09.2018 !!!
% np_fge=1; dHsel=+0.0; dDsel=+0.0; dZsel=0.0;      % 1 точку в момент измерения GSM
% np_fge=3; dHsel=+0.0; dDsel=+0.0; dZsel=0.0;      % 1 точку в момент измерения GSM и еще 2 следующих

% np_fge=0; dHsel=+0.0; dDsel=+0.0; dZsel=0.0;      % все данные за 26-30.09.2024 (GSm-90 не работал)
np_fge=4; dHsel=+0.0; dDsel=+0.0; dZsel=0.0;      % 1 точку в момент измерения GSM и еще 3 следующих

% принятый азимут Миры
% A_astr=dms_deg(209,52,54.0); % Main pillar "fDI" - старые данные
% A_astr=dms_deg(209,48,35.0); % по астрон. набл. Солнца со штатива 10.09.2012
A_astr=dms_deg(209,48,56); % по измерениям 8 и 12 февраля 2013 г.

% систематическая (инструментальная) погрешность ММП-203М2
% F(real) = F(Obs)+dFsys
% dFsys=-2.0; % по результатам сверки 12.09.2012

% градиент F между точкой регистрации протонного магнитометра и точкой DI-наблюдений
%dFgrad=+68.5; % F(DI) = F(F) + dFgrad	- из 08_06_03.xls (D14)

% амплитуда единичного выброса в данных GSM (nT), 
% который будет заменен на NaN
dF_spike_lim_GSM90=5;

% Температурные коэффициенты FGE (H,Z - в нТл/С, D - в град./C)
% Hcorr=Hobs-KtH*(TR-TR0)
KtH=-0.08; KtZ=-0.20; KtD=-0.13/60;
% Температура приведения (нулевой поправки)
TR0=10;

% Матрица коэфф. лин. зависимости между компонентами
% dZ = R(3,1)*dH + R(3,2)*dD
% компоненты - вариации в нТл, оценка по сравнению с dIdD 29.06.2013
% и по F(FGE)-F(GSM) - минимизация
%     H        D       Z
R=[+1.0000 +0.0000 +0.0000;
   +0.0000 +1.0000 +0.0000;
   +0.0130 -0.0050 +1.0000];
% R=[+1.0000 +0.0000 +0.0000;
%    +0.0000 +1.0000 +0.0000;
%    +0.0000 +0.0000 +1.0000];
% Средние вариации (по данным за 01-31 мая 2013 г)
dHmean=-130; dDmean=-460; dZmean=+205;
