% просмотр фрагментов в файле pev

clear, close all

run ..\..\ini_parameters\fge_ini.m     % инициализационные параметры

file='v2025255_20250913_1428LT.pev';

file_fge=['d:\Paratunka\data\fge_gsm\momwe\fge\2025\pev\' file];

load frg_mat_png\FRG_20250913_dHdDdZ_ini.mat Tj dHj dDj dZj

% исходные значения FGE, 2 Гц
[Tv2,dH2,dD2,dZ2,Ttr,TR]=var_fge_read_fun(file_fge);

% kjump=find(diff(Tv2)*86400<0); % индексы скачков
% N=length(kjump);

% любые скачки
kj=find(~between_my(diff(Tv2)*86400,0.49,0.51));
kja=[1; kj+1; length(Tv2)+1]; % скачки с первой и последней точкой массива
N=length(kja);

figure(1)
clf, hold on 
% for k=1:length(kja)-1
for k=1:N-1
    k1=kja(k);k2=kja(k+1)-1;
    disp([k1 k2 sum(~isnan(dH2(k1:k2)))])
    hhh=plot(Tv2(k1:k2),dH2(k1:k2),'.-');
    text(Tv2(k1),nanmean(dH2)+k*2,sprintf('%d',k),'Color',hhh.Color)
end
datetick, hold off

figure(2)
plot(dH2,'.-')

% return

Tc=[]; dHc=[]; dDc=[]; dZc=[];
k1=171961; k2=172800;
Tc=[Tc; Tv2(k1:k2)-1]; dHc=[dHc; dH2(k1:k2)]; dDc=[dDc; dD2(k1:k2)]; dZc=[dZc; dZ2(k1:k2)];
k1=254; k2=17880;
Tc=[Tc; Tv2(k1:k2)]; dHc=[dHc; dH2(k1:k2)]; dDc=[dDc; dD2(k1:k2)]; dZc=[dZc; dZ2(k1:k2)];

% dTd=fix(nanmean(Tc)-fix(nanmean(Tj)));
% Tc=Tc-dTd;

figure(10)
plot(Tc,dHc,'.-');datetick

% return

% сдвиг между кривыми до 1 сут.
dTd=fix(nanmean(Tc)-fix(nanmean(Tj)));
Tv=Tc-dTd; dHv=dHc; dDv=dDc; dZv=dZc;

figure(10)
% plot(Tv,dHv,'.-',Tj,dHj-nanmean(dHj)-120),datetick
plot(Tv,dHv-nanmean(dHv),'.-',Tj,dHj-nanmean(dHj)-1),datetick
legend('FGE','FRG')

% return

% определяем сдвиг в пределах суток
Tfge=datenum('13-Sep-2025 00:36:20');
Tfrg=datenum('13-Sep-2025 00:37:33');
dTh1=Tfge-Tfrg;
% Tfge=datenum('11-Sep-2025 23:11:08');
% Tfrg=datenum('11-Sep-2025 22:36:19');
% dTh2=Tfge-Tfrg; % после 03UT
% 
Tv=Tv-dTh1;
% k2=Tv>datenum(2025,09,11,20,55,00);
% Tv(k2)=Tv(k2)-dTh2;

plot(Tv,dHv-nanmean(dHv),'.-',Tj,dHj-nanmean(dHj)-1),datetick
legend('FGE','FRG')

% return

hhh=figure(2);
set(hhh,'Position',get(0,'Screensize'));
plot(Tv,dHv-nanmean(dHv),Tj,dHj-nanmean(dHj)),datetick
title([datestr(nanmean(Tv),'dd.mm.yyyy') '   FRG vs FGE   dH'])
ylabel('centered dH, nT'), legend('FGE','FRG')
xlabel([datestr(min([min(Tv) min(Tj)]),'dd.mm.yyyy HH:MM:SS') ...
       ' --- ' ...
        datestr(max([max(Tv) max(Tj)]),'dd.mm.yyyy HH:MM:SS')])
file_fig=['fge_from_pev_' datestr(min(Tv),'yyyymmdd_HHMMSS') '-' datestr(max(Tv),'yyyymmdd_HHMMSS') '_H.png'];
fig=gcf; fig.PaperPositionMode = 'auto';
print(file_fig, '-dpng', '-r0')

hhh=figure(3);
set(hhh,'Position',get(0,'Screensize'));
plot(Tv,dDv-nanmean(dDv),Tj,dDj-nanmean(dDj)),datetick
title([datestr(nanmean(Tv),'dd.mm.yyyy') '   FRG vs FGE   dD'])
ylabel('centered dD, nT'), legend('FGE','FRG')
xlabel([datestr(min([min(Tv) min(Tj)]),'dd.mm.yyyy HH:MM:SS') ...
       ' --- ' ...
        datestr(max([max(Tv) max(Tj)]),'dd.mm.yyyy HH:MM:SS')])
file_fig=['fge_from_pev_' datestr(min(Tv),'yyyymmdd_HHMMSS') '-' datestr(max(Tv),'yyyymmdd_HHMMSS') '_D.png'];
fig=gcf; fig.PaperPositionMode = 'auto';
print(file_fig, '-dpng', '-r0')
    
hhh=figure(4); set(hhh,'Position',get(0,'Screensize'));
plot(Tv,dZv-nanmean(dZv),Tj,dZj-nanmean(dZj)),datetick
title([datestr(nanmean(Tv),'dd.mm.yyyy') '   FRG vs FGE   dZ'])
ylabel('centered dZ, nT'), legend('FGE','FRG')
xlabel([datestr(min([min(Tv) min(Tj)]),'dd.mm.yyyy HH:MM:SS') ...
       ' --- ' ...
        datestr(max([max(Tv) max(Tj)]),'dd.mm.yyyy HH:MM:SS')])
file_fig=['fge_from_pev_' datestr(min(Tv),'yyyymmdd_HHMMSS') '-' datestr(max(Tv),'yyyymmdd_HHMMSS') '_Z.png'];
fig=gcf; fig.PaperPositionMode = 'auto';
print(file_fig, '-dpng', '-r0')

commandwindow

% выводим строку для формирования имени файлов
str_to_png=['fge_from_pev_' datestr(min(Tv),'yyyymmdd_HHMMSS') '-' datestr(max(Tv),'yyyymmdd_HHMMSS') '_C.png'];
disp(str_to_png)
comments=[];
str_to_mat=['fge_from_pev_' datestr(min(Tv),'yyyymmdd_HHMMSS') '-' datestr(max(Tv),'yyyymmdd_HHMMSS') '_dHdDdZ.mat Tv dHv dDv dZv'];
disp(str_to_mat)
