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

from obs251 over memeory stick mainly Gan

parent 6b25b023
......@@ -21,17 +21,17 @@
% attachment: same file
clear all
startdate = jd2000(2018, 05, 17); % HERE START DATE FOR non-online
enddate = jd2000(2018, 05, 17); % HERE END DATE FOR non-online
startdate = jd2000(2018, 04, 01); % HERE START DATE FOR non-online
enddate = jd2000(2018, 05, 31); % HERE END DATE FOR non-online
station = 'VNA1';
%station = 'VNA1';
%station = 'WNG1';
%station = 'BFO1';
%station = 'TDC1';
%station = 'NGK0';
station = 'NGK0';
IAGAtype = 'P'; % preliminary data
%IAGAtype = 'Q'; % quasi definitive data
%IAGAtype = 'P'; % preliminary data
IAGAtype = 'Q'; % quasi definitive data
%IAGAtype = 'P'; % quasi definitive data
recipient = {'gfzdata@swdc234.kugi.kyoto-u.ac.jp'};
......
......@@ -375,6 +375,7 @@ if exist(DI_file, 'file') == 2;
%
if strcmp(station(5), 'd') || strcmp(station(5), 's')...
|| strcmp(station(5), 'g') || strcmp(station(5), 'l')...
|| strcmp(station(5), 'j')...
&& exist(cdf_file, 'file') == 2
display('arrived here')
cdfID = cdflib.open(cdf_file);
......@@ -959,6 +960,48 @@ if exist(DI_file, 'file') == 2;
% Azimut = 201.6794? = 224.0883 gon = hier 24.0883
end
if strcmp(station(1:3), 'GAN') %
mire_nr = raw1{8,6};
pillar_nr = raw1{13,6};
deltaI = 0.0000 + 0.0000;
deltaD = 0.0000 + 0.0000;
if timeDI(1) < jd2000(2013, 12, 31)
error('no F difference in GAN before 2014 in Matlab')
elseif timeDI(1) < jd2000(2015, 2, 14, 5)
deltappm = 0.0 + 26.2; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 1, 1)
deltappm = 0.0 + 29.4; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 3, 1)
deltappm = 0.0 + 29.2; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 4, 1)
deltappm = 0.0 + 28.95; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 5, 1)
deltappm = 0.0 + 28.5; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 6, 1)
deltappm = 0.0 + 28.25; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 7, 1)
deltappm = 0.0 + 28.0; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2016, 8, 1)
deltappm = 0.0 + 27.75; % instrument_diff = 0, site diff OK
elseif timeDI(1) < jd2000(2017, 12, 23)
deltappm = 0.0 + 27.5; % instrument_diff = 0, site diff OK
else
deltappm = 0.0 + 26.8; % instr_diff = 0, site diff prelim 2018
end %if timeDI(1)
if mire_nr == 1
mire_azim = 17.8243; % Light post, 502 m from pillar 1
% 16 deg 02.51', see email Muslim 2018-05-30
% = 16.04183 deg = 17.8243 gon
elseif mire_nr == 2
mire_azim = 12.2248; % Light post, 492 m from pillar 1
% 11 deg 00.14', see email Muslim 2018-05-30
% = 11.00233 deg = 12.2248 gon
else
error('Wrong azimuth mark')
end
end
if strcmp(station(1:3), 'VSS') %
mire_nr = raw1{8,6};
pillar_nr = raw1{13,6};
......@@ -1761,8 +1804,15 @@ end
code(1:2) = '43'
elseif strcmp(upper(operator), 'JC')%Jorge Cruz 2018 SMA
code(1:2) = '44'
elseif strcmp(upper(operator), 'AND')%Horcio Andr 2018 SMA
elseif strcmp(upper(operator), 'AND')%Horcio Andr 2018 SMA
code(1:2) = '45'
elseif strcmp(upper(operator), 'AMU')||...
strcmp(upper(operator), 'MU') % Ahmed Muslim, GAN
code(1:2) = '46'
elseif strcmp(upper(operator), 'IZI')% Ibrahim Zijad, GAN
code(1:2) = '47'
elseif strcmp(upper(operator), 'NAX')% NAXMI, GAN
code(1:2) = '48'
end
%Danieal Armbruster vna 2014
%Johannes Lohse vna 2014
......@@ -1933,6 +1983,11 @@ end
Z10 = -25128; %17133.56 32257.05 -24913.32
end
if strcmp(station(1:4), 'GAN0') % IMPROVE!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
H10 = 37670; % ungefaehr
Z10 = -13330; % ungefaehr
end
if strcmp(station(1:4), 'VSS0') % more or less ok
H10 = 18420;
Z10 = -14345; %17133.56 32257.05 -24913.32
......
......@@ -19,8 +19,9 @@ clear all
close all
figure
station='GAN0j';
%station = 'SMA0g'; % c or d or g or s or l for datalogger
station = 'WNG5g'; % c or d or g or s or l for datalogger
%station = 'WNG5g'; % c or d or g or s or l for datalogger
%station = 'NGK0g'; % c or d or g or s or l for datalogger
%station = 'TDC1d'; % c or d or g or s or l for datalogger
%station = 'VNA1s'; % c or d or g or s or l for datalogger
......@@ -37,9 +38,9 @@ timeshift_array = [-300:1:300]';
number_DI = 1000; % read in all DI-measurements for this year
number_DI = 5; % read in the last x measurements in this year
number_DI = 0; % use date range from below
year = 2018;
start_mmdd = '0509';
end___mmdd = '0509';
year = 2017;
start_mmdd = '1201';
end___mmdd = '1231';
......@@ -84,7 +85,6 @@ end
%month = 1; removed 27/2-2012
%day = 7;
% % START Identify and convert Achims GFZ-DI-schemes
% % with name e.g. ng0_140102.xls
% %
......@@ -186,7 +186,8 @@ if strcmp(station(1:3), 'TDC') || strcmp(station(1:3), 'VNA') ...
|| strcmp(station(1:3), 'NGK') || strcmp(station(1:3), 'WNG') ...
|| strcmp(station(1:3), 'TTB') || strcmp(station(1:3), 'VRE')...
|| strcmp(station(1:3), 'SMA') ...
|| strcmp(station(1:3), 'SHE') || strcmp(station(1:3), 'VSS');
|| strcmp(station(1:3), 'SHE') || strcmp(station(1:3), 'VSS') ...
|| strcmp(station(1:3), 'GAN');
filename_criterion = [path_DI 'DI_' station(1:3) '_'...
num2str(year, '%4.4d') '*.xls'];
......
......@@ -39,9 +39,9 @@
clear all
pack
% supply
station = 'WNG1g'; % HIER IST DIE DATENLCKE
station = 'NGK0g'; % HIER IST DIE DATENLCKE
%with
supplementer = 'WNG2g'; % VON HIER DATEN HOLEN
supplementer = 'NGK6g'; % 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 = 'WNG2g'; % VON HIER DATEN HOLEN
% reading of SUP files!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! COULD BE
% VERY EASY
start_supp = jd2000(2014, 10, 20, 11 + 55 / 60) %Stunde Minute/60 INKLUSIVE
end___supp = jd2000(2014, 10, 20, 12 + 20 / 60) %Stunde Minute/60 EXKLUSIVE
start_supp = jd2000(2018, 05, 08, 10 + 58 / 60) %Stunde Minute/60 INKLUSIVE
end___supp = jd2000(2018, 05, 08, 11 + 05 / 60) %Stunde Minute/60 EXKLUSIVE
% chooses data investigation interval to determine linear relationship
% or user select a start and end data
......
......@@ -39,7 +39,7 @@ clear all
close all
pack
% supply
station = 'WNG1g';
station = 'NGK0g';
%with
supplementer1 = 'file'; % supplementer ist stationsname z. Bsp. 'THL4d'
% oder 'file', dann wird die station aus den
......@@ -47,12 +47,12 @@ supplementer1 = 'file'; % supplementer ist stationsname z. Bsp. 'THL4d'
% der SUP-Datei geholt., z. Bsp. 'NGK6g' aus
% 'definitive_NGK6g'
start_day = jd2000(2014, 01, 01);
end_day = jd2000(2017, 12, 31);
start_day = jd2000(2018, 05, 01);
end_day = jd2000(2018, 05, 31);
plotfigures = 0; % 1 = 4 Bilder pro Tag, nur bei max 3 Tagen? setzen.
plotfigures = 0; % 1 = 4 Bilder pro Interval, nur bei max 3 Tagen? setzen.
SUP_file = 'O:\jmat\cfg\SUP_WNG0g.2014.txt';
SUP_file = ['O:\jmat\cfg\SUP_' station '.2014.txt';]
......
......@@ -4,13 +4,14 @@ close all
%station = 'TDC1d';
station = 'NGK0g';
%station = 'GAN0j';
%station = 'WNG1g';
%station = 'VNA1s';
%station = 'BFO1s';
%calculate IAF files for timeinterval:
timeinterval = [2015; 1 ; 36]; %start year, start month, number of months
timeinterval = [2017; 1 ; 17]; %start year, start month, number of months
% option to apply time correction
time_correction_str = 'n';
......
......@@ -19,9 +19,10 @@
clear all
IAFversion = '2.10';
%station = 'GDH6d';
station = 'WNG1';
%station = 'TDC1';
%station = 'WNG1';
% station = 'GAN0';
%station = 'TDC1';
station = 'NGK0';
%station = 'SHE0';
%station = 'VNA1';
%station = 'BFO1';
......@@ -29,7 +30,7 @@ station = 'WNG1';
%calculate IAF files for timeinterval:
timeinterval = [2014; 9; 2]; %start year, start month, number of months
timeinterval = [2018; 5; 1]; %start year, start month, number of months
% works for 2011 data onwards--
......@@ -59,6 +60,7 @@ end
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
2018, 3, 16, 0 2018, 4, 5, 0
3000, 1, 1, 0, 3000, 1, 2, 0];%
end
if strcmp(station(1:3), 'NGK') % Hauptsystem in NGK0, use6 ist NGK1
......@@ -76,6 +78,9 @@ end
if strcmp(station(1:3), 'SHE')
use6 = [3000, 1, 1, 0, 3000, 1, 2, 0];
end
if strcmp(station(1:3), 'GAN')
use6 = [3000, 1, 1, 0, 3000, 1, 2, 0];
end
%what with TDC? has 1 instead of 4
......@@ -108,13 +113,13 @@ if strcmp(station(1:3), 'TDC') || strcmp(station(1:3), 'VNA')...
path_cdf4 = [path2 bs station(1:3) '1' bs 'cdf' bs];
elseif strcmp(station(1:3), 'NGK') || ...
strcmp(station(1:3), 'VSS') || strcmp(station(1:3), 'TTB')...
|| strcmp(station(1:3), 'SHE')
|| strcmp(station(1:3), 'SHE') || strcmp(station(1:3), 'GAN')
path_cdf4 = [path2 bs station(1:3) '0' bs 'cdf' bs];
end
path_cdf6 = [path2 bs station(1:3) '6' bs 'cdf' bs];
if strcmp(station(1:3), 'NGK') ...
|| strcmp(station(1:3), 'TTB')
|| strcmp(station(1:3), 'TTB') || strcmp(station(1:3), 'TTB')
path_cdf6 = [path2 bs station(1:3) '1' bs 'cdf' bs];
elseif strcmp(station(1:3), 'WNG')
path_cdf6 = [path2 bs station(1:3) '0' bs 'cdf' bs];
......@@ -170,7 +175,7 @@ for timeintervalmonth = 1:timeinterval(3)
elseif strcmp(station(1:3), 'NGK') ...
|| strcmp(station(1:3), 'VSS') || strcmp(station(1:3), 'TTB')...
|| strcmp(station(1:3), 'SHE')
|| strcmp(station(1:3), 'SHE') || strcmp(station(1:3), 'GAN')
min_file4 = [path_cdf4 station(1:3) '0_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
......@@ -911,8 +916,8 @@ for timeintervalmonth = 1:timeinterval(3)
if strcmp(station(1:3), 'BFO') % NEWEST ONE, USE TO COPY
% ' BFO'
if strcmp(station(1:3), 'BFO')
% ' BFO' to 'OFB ' = 4F464220
%converts number string but jrgm not sure if ' necessary
A1(j*5888+1) = hex2dec('4F464220');
......@@ -997,6 +1002,100 @@ for timeintervalmonth = 1:timeinterval(3)
A1(j*5888+5876) = int32(10*roundn(Gdmv,-1));
%Gmin % just for display
end %if BFO
if strcmp(station(1:3), 'GAN') % NEWEST ONE, USE TO COPY
% ' GAN' to 'NAG ' = 4E414720
% ' BFO' to 'OFB ' = 4F464220
%converts number string but jrgm not sure if ' necessary
A1(j*5888+1) = hex2dec('4E414720');
% year plus doy
A1(j*5888+2) = year*1000 + ...
jd2000(year, month, day) - jd2000(year-1, 12, 31);
% colat in degree * 1000 = 90 - (-0.695) deg = 1000 * 90.695
% = 90695
A1(j*5888+3) = int32(090695);
% eastlongitude in degree *1000 = 73.154*1000 = 73154
%
A1(j*5888+4) = int32(073154);
%elevation in m is 2
A1(j*5888+5) = int32(2);
if strcmp(IAFversion, '2.00')
%orientation of data in file XYZG = GZYX = '475A5958'
%(always G, even if no PPM)
A1(j*5888+6) = hex2dec('475A5958');
end
if strcmp(IAFversion, '2.10')
%orientation XYZG = GZYX = '475A5958' (if PPM)
%orientation ' XYZ' = 'ZYX ' = '5A595820'(if no PPM)
A1(j*5888+6) = hex2dec('475A5958');
end
%DATA Source ' GFZ' = 'ZFG ', '5A464720'I keep GFZ here
%DATA Source GAN:
% ' ETH' = 'HTE ' = 48544520
% as we produce the definitive data
A1(j*5888+7) = hex2dec('48544520');
% H conversion factor not used for HDZ, set to 10000 for XYZ
A1(j*5888+8) = int32(10000);
%Data quality IMAG = GAMI = '47414D49'
A1(j*5888+9) = hex2dec('47414D49');
% Instrumentation ' LC' = 'CL ' = '434C2020'
A1(j*5888+10) = hex2dec('434C2020');
%K9-value in nT %%%%%% SHALL ADOPT USEFUL VALUE FOR K9
% K9 is 500 for WNG, NGK, FUR
% GAN has no K-values, set K9 = 0
A1(j*5888+11) = int32(0);
% sampling rate in ms, GAN has 1 Hz
A1(j*5888+12) = int32(1000);
% sensor orientation HDZF = FZDH = '465A4448'
A1(j*5888+13) = hex2dec('465A4448');
% % publication date here YYMM = MMYY
% A1(j*5888+14) = hex2dec('59594D4D');
% publication date here actual date of file generation
% decimal 0 to 9 is hex 30 to 39
% find current date and store in publication date string
pudas = datestr(now, 'mmyy'); % e.g '1199'
%then we should write '09090101'
publ_date_hex_str = ['3' pudas(4) '3' pudas(3) '3'...
pudas(2) '3' pudas(1) ];
A1(j*5888+14) = hex2dec(publ_date_hex_str);
if strcmp(IAFversion, '2.00')
% reserved for future, since 2009? (2.0?) use 2 = format version
% 2.0
A1(j*5888+15) = int32(2);
end
if strcmp(IAFversion, '2.10')
% reserved for future, since 2010 (2.1) use 3 = format version 2.1
A1(j*5888+15) = int32(3);
end
% reserved for institution, NOW RESERVED FOR FUTURE USE
A1(j*5888+16) = int32(0);
% minute mean of G, if observatory has recording scalar don't do
% anything (keep values written to A1 from minute mean section),
% else is observatory as no recording scalar set Gmin to 88888.8
% and redo writting values to array A1:
% Gmin = 88888.8;
% A1(j*5888+4337:j*5888+5776) = int32(10*roundn(Gmin,-1));
% hourly and daily mean of G
if strcmp(IAFversion, '2.00')
Gmin(isnan(Gmin)) = 99999.9; %only NaN set to 9er
Ghmv1(:,1) = 99999.9; %all 24 values set to 9er
Gdmv = 99999.9; %daily mean
end
if strcmp(IAFversion, '2.10')
Gmin(isnan(Gmin)) = 99999.9; %only NaN set to 9er or 8er
Ghmv1(:,1) = 99999.9; %all 24 values set to 9er or 8er
Gdmv = 99999.9; % daily mean, 9er or 8er
end
A1(j*5888+4337:j*5888+5776) = int32(10*roundn(Gmin,-1));
A1(j*5888+5849:j*5888+5872) = int32(10*roundn(Ghmv1,-1));
A1(j*5888+5876) = int32(10*roundn(Gdmv,-1));
%Gmin % just for display
end %if GAN
clear tmin Hmin Dmin Zmin Xmin Ymin...
......
This diff is collapsed.
......@@ -128,6 +128,75 @@ function [H0alines, D0alines, Z0alines] = ...
if strcmp(station, 'GAN0j')
H0alines = {
5113 37671.30 'p' 'very prelim'
5114 37671.30 'd' '2014-01-01 very prelim'
5238 37670.30 'd' ''
5257 37672.5 'd' ''
5357 37672.5 'd' ''
5400 37673.2 'd' ''
5420 37673.0 'd' ''
5520 37672.4 'd' ''
5548 37672.4 'd' '2015-03-11 FGE rotated'
5548 37672.9 'd' '2015-03-11 FGE rotated'
5600 37672.8 'd' ''
5700 37672.4 'd' ''
5790 37672.6 'd' ''
5840 37672.5 'd' ''
6050 37669.3 'd' ''
6100 37669.3 'd' ''
6300 37668.5 'd' ''
6420 37668.3 'd' ''
6575 37667.5 'd' '2018-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 -4.50600 'p' 'very prelim'
5114 -4.50600 'd' '2014-01-01 very prelim'
5315 -4.51100 'd' ''
5350 -4.51000 'd' ''
5370 -4.50900 'd' ''
5385 -4.50900 'd' ''
5385 -4.51000 'd' ''
5420 -4.51000 'd' ''
5467 -4.51220 'd' ''
5495 -4.50820 'd' ''
5548 -4.51200 'd' '2015-03-11 FGE rotated'
5548 -4.93950 'd' '2015-03-11 FGE rotated'
6000 -4.94600 'd' ''
6150 -4.94600 'd' ''
6220 -4.94500 'd' ''
6280 -4.94500 'd' ''
6310 -4.94600 'd' ''
6450 -4.94600 'd' ''
6500 -4.94450 'd' ''
6575 -4.94450 'd' '2018-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 -13329.50 'p' 'very prelim'
5114 -13329.50 'p' '2014-01-01 very prelim'
5250 -13329.50 'p' ''
5257 -13326.30 'p' ''
5350 -13327.20 'p' ''
5360 -13327.20 'p' ''
5390 -13326.30 'p' ''
5410 -13326.30 'p' ''
5430 -13327.70 'p' ''
5445 -13327.00 'p' ''
5460 -13327.20 'p' ''
5485 -13328.70 'p' ''
5548 -13328.70 'p' ''
5548 -13327.60 'p' '2015-03-11 FGE rotated'
5660 -13327.00 'p' ''
6290 -13327.20 'p' ''
6400 -13325.20 'p' ''
6480 -13326.10 'p' ''
6550 -13326.10 'p' ''
6575 -13325.00 'p' '2018-01-01 very prelim'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......@@ -138,10 +207,13 @@ if strcmp(station, 'WNG0g')
%
%
H0alines = {
5113 18184.00 'p' ''
5114 18184.00 'p' '2014-01-01'
5435.75 18184.00 'p' ''
5435.75 18093.00 'p' 'datum_Sensor_ausgetauscht'
5113 18184.00 'p' ''
5114 18184.00 'p' '2014-01-01'
5435.75 18184.00 'p' ''
5435.75 18093.00 'p' 'datum_Sensor_ausgetauscht'
6575 18093.00 'p' '2018-01-01'
6600 18091.00 'p' '2018-01-26'
6665 18092.90 'p' '2018-03-17'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 0.63400 'p' ''
......@@ -154,6 +226,10 @@ if strcmp(station, 'WNG0g')
5114 45984.50 'p' '2014-01-01'
5435.75 45984.50 'p' 'datum_Sensor_ausgetauscht'
5435.75 46142.00 'p' ''
6210 46142.00 'p' '2017-01-01'
6574 46142.00 'p' '2017-12-31'
6634 46143.00 'p' '2018-03-01'
6665 46144.00 'p' '2018-04-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......@@ -331,6 +407,25 @@ if strcmp(station, 'WNG2g')
}; % REMEMBER for QD: default is the first day of the next month!!!
end
if strcmp(station, 'WNG2g')
%
% NOTE:
%
%
H0alines = {
4017 37650.00 'p' ''
4018 37650.00 'p' '2011-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
4017 -4.000 'p' ''
4018 -4.000 'p' '2011-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
4017 -13250.00 'p' ''
4018 -13250.00 'p' '2011-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
if strcmp(station, 'NGK0g')
%
......@@ -365,9 +460,11 @@ if strcmp(station, 'NGK0g')
6453 18929.3 'q' '2017-09-01'
6483 18929.3 'q' '2017-10-01'
6524 18929.3 'q' '2017-11-11'
6575 18929.3 'q' '2018-01-01'
6634 18929.5 'q' '2018-03-01'
6665 18929.5 'q' '2018-04-01'
6695 18929.5 'q' '2018-05-01'
6726 18929.5 'q' '2018-06-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 2.45166 'p' ''
......@@ -395,9 +492,11 @@ if strcmp(station, 'NGK0g')
6279.7 2.44900 'q' '2017-03-11 kleiner Sprung ADC'
6422 2.44600 'q' '2017-08-01'
6525 2.44450 'q' '2017-11-12'
6575 2.44500 'q' '2018-01-01'
6634 2.44650 'q' '2018-03-01'
6665 2.44700 'q' '2018-04-01'
6695 2.44900 'q' '2018-05-01'
6726 2.44900 'q' '2018-06-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 45426.17 'p' ''
......@@ -428,9 +527,11 @@ if strcmp(station, 'NGK0g')
6330 45425.55 'q' '2017-05-01'
6483 45425.55 'q' '2017-10-01'
6524 45425.55 'q' '2017-11-11'
6575 45425.65 'q' '2018-01-01'
6634 45425.75 'q' '2018-03-01'
6665 45425.90 'q' '2018-04-01'
6695 45425.80 'q' '2018-05-01'
6726 45425.50 'q' '2018-06-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......@@ -443,16 +544,20 @@ if strcmp(station, 'NGK1g')
H0alines = {
5113 18864.00 'p' ''
5114 18864.00 'p' '2014-01-01'
5479 18864.20 'p' ''
5479 18864.20 'p' '2015-01-01'
5844 18864.20 'p' '2016-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 2.22800 'p' ''
5114 2.22800 'p' '2014-01-01'
5479 2.22800 'p' ''
5113 2.22400 'p' ''
5114 2.22400 'p' '2014-01-01'
5479 2.22400 'p' '2015,01,01'
5844 2.22000 'p' '2016,01,01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 45419.00 'p' ''
5114 45419.00 'p' '2014-01-01'
5114 45419.00 'p' '2014-01-01'
5479 45418.50 'p' '2015-01-01'
5844 45418.50 'p' '2016-01-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......@@ -481,18 +586,24 @@ if strcmp(station, 'SMA0g')
%
H0alines = {
3653 27230.00 'p' '2010-01-01'
5114 27552.00 'p' '2014-01-01'
6680 27552.00 'p' '2018-04-16'
5114 27559.00 'p' '2014-01-01'
6680 27559.00 'p' '2018-04-16'
6704.375 27559.00 'p' '2018-05-10 Stromunterbrechung, neue Basis'
6704.375 27559.00 'p' '2018-05-10 Stromunterbrechung, neue Basis'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
3653 -5.305 'p' '2010-01-01'
5114 -4.99 'p' '2014-01-01'
6680 -4.99 'p' '2018-04-16'
3653 -5.3050 'p' '2010-01-01'
5114 -5.0033 'p' '2014-01-01'
6680 -5.0033 'p' '2018-04-16'
6704.375 -5.0033 'p' '2018-05-10 Stromunterbrechung, neue Basis'
6704.375 -4.9967 'p' '2018-05-10 Stromunterbrechung, neue Basis'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
3653 36075.00 'p' '2010-01-01'
5114 35780.00 'p' '2014-01-01'
6680 35780.00 'p' '2018-04-16'