% рисуем градиент в абсол. пав-не (между DI-столбом и точкой установки
% скалярного магнитометра
% Для статьи

clear, close all

Tout=datenum(2016,07,01);   % дата, с которой рисуем и выводим

tmp=load('d:\Paratunka\ini_parameters\gradF_Obs_Abs_hut.dat');
th=tmp(:,4); th(isnan(th))=0;
Tdfo=datenum(tmp(:,1:3))+th/24;
dFo=tmp(:,7); obs=tmp(:,8);

% убираем скачок из-за смены dIdD на POS-1 и смены места 16 нТл
k=Tdfo<datenum(2015,06,29);
dFo(k)=dFo(k)-15.8;

% выбраковка
k=obs>0; Tdfo_err=Tdfo(~k); dFo_err=dFo(~k); % выбросы

% сплайн-аппроксимация на каждые сутки
% p=0.0001;
p=0.001;
% Tdfa=min(Tdfo):max(Tdfo); Tdfa=Tdfa(:);
Tdfa=min(fix(Tdfo)):max(fix(Tdfo)); Tdfa=Tdfa(:);
dFa=csaps(Tdfo(k),dFo(k),p,Tdfa);

ko=Tdfo>datenum(2016,07,01) & (obs>1 & obs<9);
plot(Tdfo(ko),dFo(ko),'k.'),datetick
ylabel('dFgrad, nT')
