Commit 168c79ee authored by Jürgen Matzka's avatar Jürgen Matzka
Browse files

version von obs251, neu SUP GFZ

parent b203f8e4
......@@ -21,8 +21,8 @@
% attachment: same file
clear all
startdate = jd2000(2018, 05, 14); % HERE START DATE FOR non-online
enddate = jd2000(2018, 05, 14); % HERE END DATE FOR non-online
startdate = jd2000(2018, 05, 15); % HERE START DATE FOR non-online
enddate = jd2000(2018, 05, 16); % HERE END DATE FOR non-online
station = 'VNA1';
%station = 'WNG1';
......
......@@ -39,9 +39,9 @@
clear all
pack
% supply
station = 'NGK0g'; % HIER IST DIE DATENLCKE
station = 'WNG1g'; % HIER IST DIE DATENLCKE
%with
supplementer = 'NGK1g'; % VON HIER DATEN HOLEN
supplementer = 'WNG2g'; % VON HIER DATEN HOLEN
% supplementing interval, it is recommended to use a start time on the
% previous day instead of hour 00 and an end date on the day following day
......@@ -51,8 +51,8 @@ supplementer = 'NGK1g'; % VON HIER DATEN HOLEN
% reading of SUP files!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! COULD BE
% VERY EASY
start_supp = jd2000(2017, 03, 11, 15 + 15 / 60) %Stunde Minute/60 INKLUSIVE
end___supp = jd2000(2017, 03, 11, 17 + 20 / 60) %Stunde Minute/60 EXKLUSIVE
start_supp = jd2000(2015, 07, 17, 6 + 46 / 60) %Stunde Minute/60 INKLUSIVE
end___supp = jd2000(2015, 07, 17, 6 + 58 / 60) %Stunde Minute/60 EXKLUSIVE
% chooses data investigation interval to determine linear relationship
% or user select a start and end data
......
......@@ -47,8 +47,8 @@ supplementer1 = 'file'; % supplementer ist stationsname z. Bsp. 'THL4d'
% der SUP-Datei geholt., z. Bsp. 'NGK6g' aus
% 'definitive_NGK6g'
start_day = jd2000(2014, 01, 16);
end_day = jd2000(2018, 01, 01);
start_day = jd2000(2014, 01, 01);
end_day = jd2000(2018, 12, 31);
plotfigures = 0; % 1 = 4 Bilder pro Tag, nur bei max 3 Tagen? setzen.
......@@ -94,6 +94,9 @@ clear B timeB
endi = i;
end
end
if isnan(starti) && start_day < timestartB(1)
starti = 1;
end %if
jd2date(timestartB(starti:endi,1))
......
......@@ -3,14 +3,14 @@ clear all
close all
%station = 'TDC1d';
%station = 'NGK0g';
station = 'WNG1g';
station = 'NGK0g';
%station = 'WNG1g';
%station = 'VNA1s';
%station = 'BFO1s';
%calculate IAF files for timeinterval:
timeinterval = [2014; 1 ; 48]; %start year, start month, number of months
timeinterval = [2015; 1 ; 36]; %start year, start month, number of months
% option to apply time correction
time_correction_str = 'n';
......
......@@ -29,7 +29,7 @@ station = 'WNG1';
%calculate IAF files for timeinterval:
timeinterval = [2014; 1; 12]; %start year, start month, number of months
timeinterval = [2014; 1; 1]; %start year, start month, number of months
% works for 2011 data onwards--
......@@ -56,11 +56,13 @@ end
if strcmp(station(1:3), 'VNA')
use6 = [3000, 1, 1, 0, 3000, 1, 2, 0];
end
if strcmp(station(1:3), 'WNG')
use6 = [2014, 5, 15, 0, 2014, 5, 21, 0 % needs finetuning
3000, 1, 1, 0, 3000, 1, 2, 0]; % needs finetuning
if strcmp(station(1:3), 'WNG') % Hauptsystem ist WNG1, use6 ist WNG0
use6 = [2014, 1, 17, 0, 2014, 1, 24, 0
2014, 5, 15, 0, 2014, 5, 21, 0
2015, 7, 2, 0, 2015, 7, 6, 0
3000, 1, 1, 0, 3000, 1, 2, 0];%
end
if strcmp(station(1:3), 'NGK')
if strcmp(station(1:3), 'NGK') % Hauptsystem in NGK0, use6 ist NGK1
use6 = [3000, 1, 1, 0, 3000, 1, 2, 0];
end
if strcmp(station(1:3), 'BFO')
......
......@@ -56,8 +56,8 @@ clear all
% TDC works only after 2011 01 01
% station = 'TDC1d'; startime= jd2000(2011, 1, 1);
% station = 'WNG0g', startime= jd2000(2014, 1, 1);
station = 'WNG1g', startime= jd2000(2014, 1, 1);
% station = 'WNG2g', startime= jd2000(2014, 1, 1);
% station = 'WNG1g', startime= jd2000(2014, 1, 1);
station = 'WNG2g', startime= jd2000(2014, 1, 1);
% station = 'VNA1s', startime= jd2000(2013, 11, 01);
% station = 'TTB0g'; startime= jd2000(2015, 11, 01);
% station = 'TTB1l'; startime= jd2000(2008, 06, 01);
......@@ -3442,14 +3442,154 @@ if strcmp(station(1:3), 'NGK')
end % year 2016
if write_IBFV200_for_year == 2017
adopted_S = 5.12; %
% IBFV200_filename = 'NGK2015.blv';
IBFV200_firstline = 'HDZF HHHHH FFFFF NGK 2017';
IBFV200_comments = [
'COMMENTS FOR 2017 AND DIARY ';
'Data was processed with a suite of Matlab ';
'programs developped by Juergen Matzka and Achim ';
'Morschhauser. ';
'A new feaure since 2016 is manual flagging ';
'of lightning signals. The signals remain in the ';
'raw data, but they are not used for calculating ';
'the definitive minute means. ';
'On January 13, the Zenith ocular of the main ';
'theodolite Zeiss 010B ser.nr. 154180 on pillar 8';
'was replaced by a normal ocular. ';
' ';
'PROCEDURES: ';
'Absolutes were regularly taken by Katrin Tornow ';
'and Stefan Rettig. ';
'A total of two or three absolutes per week are ';
'usual. ';
'Absolute measurements with inconsistent sensor ';
'offset or misalignement have been omitted. ';
'Some absolute measurements have been corrected ';
'for obvious misreadings when correction gives ';
'consistent sensor offset and misalignement. ';
' ';
'ADOPTION OF BASELINE: ';
'H baseline by linear segments ';
'D baseline by linear segments ';
'Z baseline by linear segments ';
' ';
'VARIOMETER DATA AND TIME STAMP ACCURACY: ';
'Data is taken from NG0 and digitsed by the 0.5Hz';
'data logger MAGDALOG developed at Niemegk. ';
'Time stamping is controlled by DCF. ';
' ';
'FILTER ';
'Minute means filtered according to INTERMAGNET ';
'filter. Note, that before 2016, the INTERMAGNET ';
'filter was not use, see readme file. ';
' ';
'SUPPLANTED DATA: ';
'Short intervals were supplanted when human ';
'activity like gardening has disturbed NG0 data. ';
' ';
'MAJOR DATA GAPS VARIOMETER: ';
'none ';
' ';
'MAJOR DATA GAPS SCALAR MAGNETOMETER: ';
'none ';
' ';
'TEMPERATURE PROBLEMS: ';
'none ';
' ';
'VARIOMETER QUALITY: ';
'good ';
' ';
'RECORDING ABSOLUTE SCALAR QUALITY: ';
'good ';
' ';
'END '];
IBFV200_discontinuity_index = [];
IBFV200_discontinuity = zeros(1,...
sum(eomday(write_IBFV200_for_year, 1:12)));
IBFV200_discontinuity(IBFV200_discontinuity_index) = 1;
end % year 2017
end % NGK
if strcmp(station(1:3), 'WNG')
if write_IBFV200_for_year == 2014
adopted_S = -3.34; % becomes -3.6 from 2017-01-01
IBFV200_firstline = 'HDZF HHHHH FFFFF WNG 2014';
IBFV200_comments = [
'COMMENTS FOR 2014 AND DIARY ';
'Data was processed with a suite of Matlab ';
'programs developped by Juergen Matzka and Achim ';
'Morschhauser. This allows to deliver definitive ';
'data filtered according to INTERMAGNET standard.';
'A new feaure since 2014 is manual flagging ';
'of lightning signals. The signals remain in the ';
'raw data, but they are not used for calculating ';
'the definitive minute means. ';
'A number of challenge were faced in 2014: ';
'Absolute measurements were of varying quality. ';
'The main theodolite was replaced in March 2015 ';
'and found to be damaged. ';
'The temperature of the main and back up ';
'variometer varied between 15 and 28 degree C. ';
' ';
'PROCEDURES: ';
'Absolutes were regularly taken by Erika Koenig ';
'and sometimes by Hans-Joachim Linthe. Occasional';
'absolutes were taken by Juergen Matzka ';
'Absolute measurements with inconsistent sensor ';
'offset or misalignement have been omitted. ';
'Some absolute measurements have been corrected ';
'for obvious misreadings when correction gives ';
'consistent sensor offset and misalignement. ';
' ';
'ADOPTION OF BASELINE: ';
'H baseline by linear segments ';
'D baseline by linear segments ';
'Z baseline by linear segments ';
' ';
'VARIOMETER DATA AND TIME STAMP ACCURACY: ';
'Data is taken from WN1 and digitsed by the 0.5Hz';
'data logger MAGDALOG developed at Niemegk. ';
'Time stamping is controlled by DCF. ';
' ';
'FILTER ';
'Minute means filtered according to INTERMAGNET ';
'filter. Note, that before 2014, the INTERMAGNET ';
'filter was not use, see readme file. ';
' ';
'SUPPLANTED DATA: ';
'Short intervals were supplanted when human ';
'activity like gardening has disturbed NG0 data. ';
' ';
'MAJOR DATA GAPS VARIOMETER: ';
'none ';
' ';
'MAJOR DATA GAPS SCALAR MAGNETOMETER: ';
'none ';
' ';
'TEMPERATURE PROBLEMS: ';
'none ';
' ';
'VARIOMETER QUALITY: ';
'good ';
' ';
'RECORDING ABSOLUTE SCALAR QUALITY: ';
'good ';
' ';
'END '];
IBFV200_discontinuity_index = [];
IBFV200_discontinuity = zeros(1,...
sum(eomday(write_IBFV200_for_year, 1:12)));
IBFV200_discontinuity(IBFV200_discontinuity_index) = 1;
end % year 2014
end % WNG
%
% observed baseline of scalar F
% usually missing value
......
......@@ -193,8 +193,8 @@ if strcmp(station, 'WNG1g')
5114 18060.00 'd' '2014-01-01'
5300 18059.00 'd' ''
5350 18059.00 'd' ''
5361 18062.00 'd' ''
5410 18062.00 'd' ''
5361 18060.50 'd' ''
5410 18060.50 'd' ''
5417 18059.00 'd' ''
5479 18059.00 'd' '2015-01-01'
5500 18059.00 'd' ''
......@@ -214,7 +214,7 @@ if strcmp(station, 'WNG1g')
6050 0.644000 'd' ''
6210 0.641000 'd' '2017-01-01'
6250 0.637000 'd' ''
6460 0.637000 'd' ''
6460 0.631000 'd' ''
6483 0.630000 'd' ''
6575 0.630000 'd' '2018-01-01'
6665 0.63000 'q' '2018-04-01'
......@@ -222,8 +222,8 @@ if strcmp(station, 'WNG1g')
Z0alines = {
5113 45847.00 'd' ''
5114 45847.00 'd' '2014-01-01'
5280 45847.00 'd' ''
5390 45845.10 'd' ''
5250 45847.00 'd' ''
5390 45846.00 'd' ''
5479 45847.40 'd' '2015-01-01'
5500 45847.70 'd' ''
5600 45846.80 'd' ''
......@@ -245,6 +245,71 @@ if strcmp(station, 'WNG1g')
}; % REMEMBER for QD: default is the first day of the next month!!!
end
% Ansatz gut fuer dF, aber nicht gut im VGL mit Niemegk, 16.5.2018
% if strcmp(station, 'WNG1g')
% %
% % NOTE:
% %
% %
% H0alines = {
% 5113 18060.00 'd' ''
% 5114 18060.00 'd' '2014-01-01'
% 5300 18059.00 'd' ''
% 5350 18059.00 'd' ''
% 5361 18062.00 'd' ''
% 5410 18062.00 'd' ''
% 5417 18059.00 'd' ''
% 5479 18059.00 'd' '2015-01-01'
% 5500 18059.00 'd' ''
% 5600 18063.70 'd' ''
% 5844 18063.70 'd' '2016-01-01'
% 6210 18064.00 'd' '2017-01-01'
% 6575 18064.00 'd' '2018-01-01'
% 6665 18065.30 'q' '2018-04-01'
% }; % REMEMBER for QD: default is the first day of the next month!!!
% D0alines = {
% 5113 0.684000 'd' ''
% 5114 0.684000 'd' '2014-01-01'
% 5253 0.684000 'd' ''
% 5253 0.639000 'd' '2014-05-20 Sprung in D in Datenluecke'
% 5479 0.640500 'd' '2015-1-01'
% 5844 0.645000 'd' '2016-01-01'
% 6050 0.644000 'd' ''
% 6210 0.641000 'd' '2017-01-01'
% 6250 0.637000 'd' ''
% 6460 0.637000 'd' ''
% 6483 0.630000 'd' ''
% 6575 0.630000 'd' '2018-01-01'
% 6665 0.63000 'q' '2018-04-01'
% }; % REMEMBER for QD: default is the first day of the next month!!!
% Z0alines = {
% 5113 45847.00 'd' ''
% 5114 45847.00 'd' '2014-01-01'
% 5280 45847.00 'd' ''
% 5390 45845.10 'd' ''
% 5479 45847.40 'd' '2015-01-01'
% 5500 45847.70 'd' ''
% 5600 45846.80 'd' ''
% 5775 45845.80 'd' ''
% 5810 45846.80 'd' ''
% 5844 45846.80 'd' '2016-01-01'
% 5996 45846.80 'd' ''
% 6040 45845.60 'd' ''
% 6130 45845.60 'd' ''
% 6180 45846.50 'd' ''
% 6210 45846.50 'd' '2017-01-01'
% 6330 45846.50 'd' ''
% 6380 45845.60 'd' ''
% 6500 45845.60 'd' ''
% 6575 45846.20 'q' '2018-01-01'
% 6606 45846.20 'q' '2018-02-01'
% 6634 45846.20 'q' '2018-03-01'
% 6665 45847.20 'q' '2018-04-01'
% }; % REMEMBER for QD: default is the first day of the next month!!!
% end
if strcmp(station, 'WNG2g')
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment