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

Version von obs251 z bsp WNG5, definitive WNG

parent d32e44f6
......@@ -21,8 +21,8 @@
% attachment: same file
clear all
startdate = jd2000(2018, 05, 04); % HERE START DATE FOR non-online
enddate = jd2000(2018, 05, 06); % HERE END DATE FOR non-online
startdate = jd2000(2018, 05, 14); % HERE START DATE FOR non-online
enddate = jd2000(2018, 05, 14); % HERE END DATE FOR non-online
station = 'VNA1';
%station = 'WNG1';
......
......@@ -192,6 +192,153 @@ end
if strcmp(station, 'WNG1g') %works only for daily files!!! % JMAT 10.5.2018
if timein(1) >= jd2000(2014, 08, 10) && timein(1) < jd2000(2014, 11, 8)
display('Temperatur-Korrektur fuer WNG, Aug. 2014 bis Nov 2014.')
%not sure if the following works
if max(diff(timein)) > 0.00005
display(['more than 2 consecutive seconds missing, check if array'...
'is suitable for smoothing temperature'])
end
%
%
%smooth temperature records
%
%
%T1dec = T1in(1:20:end); %decimate temp
%T1dec = nanmean(reshape(T1in, 60, 1440)', 2);
timedec = [timein(1); nanmean(reshape(timein, 120, 720)', 2); timein(end)];
T1dec = [nanmean(T1in(1:30)); ...
nanmean(reshape(T1in, 120, 720)', 2); nanmean(T1in(end-30:end))];
T1smoothdec = smooth(T1dec, 5);
T1smooth = interp1(timedec, T1smoothdec, timein);
% T2dec = [nanmean(T2in(1:30)); ...
% nanmean(reshape(T2in, 120, 720)', 2); nanmean(T2in(end-30:end))];
% T2smoothdec = smooth(T2dec, 5);
% T2smooth = interp1(timedec, T2smoothdec, timein);
%
%
% use the adopted temperature correction
%
%
% this is adopted temp correction, see notes:
% stable temperature at just before and after the worst event is
% 20 C, for this zero correction.
%Other temperature variations were not so visible in baelines.
% uncorrected baselines rise with 0.00255 degrees per C
% H = 18170 nT, this is 0.153 minutes of arc per C
% or 0.8 nT per C
% that means that variations are too small for higher temperatures
% and too large for lower temperatures.
% variation corretction is: HEout = HEin - (25-T1smooth)*(+0.8)
%
% compensate for temperature corection:
% HEout = HEin + (25-T1smooth)*(+0.8)
%
% note:
% 'direction' 1 == correct for temperature
% % -1 == compensate for temperature corection
if direction > 0
HEout = HEin + (T1smooth - 20.0)*(+0.8);
end
if direction < 0
HEout = HEin - ( T1smooth - 20.0)*(+0.8);
end
temp_cor_structure.HE_T1_factor = +0.8;
temp_cor_structure.HE_T1_offset = -20.0;
end
end
% % erster Versuch WNG1, zu langer Zeitraum
% if strcmp(station, 'WNG1g') %works only for daily files!!! % JMAT 10.5.2018
% if timein(1) >= jd2000(2013, 12, 31) && timein(1) < jd2000(2015, 3, 26)
%
% display('Temperatur-korrektur fuer WNG, Jan 20154 bis Maerz 2015.')
%
% %not sure if the following works
% if max(diff(timein)) > 0.00005
% display(['more than 2 consecutive seconds missing, check if array'...
% 'is suitable for smoothing temperature'])
% end
%
% %
% %
% %smooth temperature records
% %
% %
% %T1dec = T1in(1:20:end); %decimate temp
% %T1dec = nanmean(reshape(T1in, 60, 1440)', 2);
% timedec = [timein(1); nanmean(reshape(timein, 120, 720)', 2); timein(end)];
%
% T1dec = [nanmean(T1in(1:30)); ...
% nanmean(reshape(T1in, 120, 720)', 2); nanmean(T1in(end-30:end))];
% T1smoothdec = smooth(T1dec, 5);
% T1smooth = interp1(timedec, T1smoothdec, timein);
%
% % T2dec = [nanmean(T2in(1:30)); ...
% % nanmean(reshape(T2in, 120, 720)', 2); nanmean(T2in(end-30:end))];
% % T2smoothdec = smooth(T2dec, 5);
% % T2smooth = interp1(timedec, T2smoothdec, timein);
%
%
% %
% %
% % use the adopted temperature correction
% %
% %
%
% % this is adopted temp correction, see notes:
% % stable temperature at 26.3.2015 is 25 C, for this zero correction
% % uncorrected baselines rise with 0.00255 degrees per C
% % H = 18170 nT, this is 0.153 minutes of arc per C
% % or 0.8 nT per C
% % that means that variations are too small for higher temperatures
% % and too large for lower temperatures.
% % variation corretction is: HEout = HEin - (25-T1smooth)*(+0.8)
% %
% % compensate for temperature corection:
% % HEout = HEin + (25-T1smooth)*(+0.8)
% %
% % note:
% % 'direction' 1 == correct for temperature
% % % -1 == compensate for temperature corection
%
%
% if direction > 0
% HEout = HEin + (T1smooth - 25.0)*(+0.8)
%
% end
% if direction < 0
% HEout = HEin - ( T1smooth - 25.0)*(+0.8)
%
% end
% temp_cor_structure.HE_T1_factor = +0.8;
% temp_cor_structure.HE_T1_offset = -25.0;
%
%
% end
% end
%OLD STUFF, KEPT FOR DOCUMENTATION:
% if strcmp(station, 'TDC1d') %works only for daily files!!!
% if timein(1) >= jd2000(2011, 2, 1)
......
This diff is collapsed.
......@@ -19,8 +19,8 @@ clear all
close all
figure
station = 'SMA0g'; % c or d or g or s or l for datalogger
%station = 'WNG1g'; % c or d or g or s or l for datalogger
%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 = '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
......@@ -38,8 +38,8 @@ 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 = '0419';
end___mmdd = '0419';
start_mmdd = '0509';
end___mmdd = '0509';
......
......@@ -39,9 +39,9 @@
clear all
pack
% supply
station = 'WNG1g'; % HIER IST DIE DATENLCKE
station = 'NGK0g'; % HIER IST DIE DATENLCKE
%with
supplementer = 'WNG0g'; % VON HIER DATEN HOLEN
supplementer = 'NGK1g'; % 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 = 'WNG0g'; % VON HIER DATEN HOLEN
% reading of SUP files!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! COULD BE
% VERY EASY
start_supp = jd2000(2018, 03, 17, 00 + 00 / 60) %Stunde Minute/60 INKLUSIVE
end___supp = jd2000(2018, 04, 04, 23 + 58 / 60) %Stunde Minute/60 EXKLUSIVE
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
% 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 = 'NGK0g';
station = 'WNG1g';
%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(2017, 12, 01);
end_day = jd2000(2017, 12, 31);
start_day = jd2000(2014, 01, 16);
end_day = jd2000(2018, 01, 01);
plotfigures = 0; % 1 = 4 Bilder pro Tag, nur bei max 3 Tagen? setzen.
SUP_file = 'O:\jmat\cfg\SUP_NGK0g.2014.txt';
SUP_file = 'O:\jmat\cfg\SUP_WNG1g.2014.txt';
......
......@@ -10,7 +10,7 @@ station = 'WNG1g';
%calculate IAF files for timeinterval:
timeinterval = [2018; 1 ;3]; %start year, start month, number of months
timeinterval = [2014; 1 ; 48]; %start year, start month, number of months
% option to apply time correction
time_correction_str = 'n';
......
......@@ -19,9 +19,9 @@
clear all
IAFversion = '2.10';
%station = 'GDH6d';
%station = 'NGK0';
%station = 'TDC1';
station = 'WNG1';
%station = 'TDC1';
%station = 'WNG1';
%station = 'SHE0';
%station = 'VNA1';
%station = 'BFO1';
......@@ -29,8 +29,8 @@ station = 'WNG1';
%calculate IAF files for timeinterval:
timeinterval = [2018; 2; 2]; %start year, start month, number of months
timeinterval = [2018; 3; 1];
timeinterval = [2014; 1; 12]; %start year, start month, number of months
% works for 2011 data onwards--
% use 0 UT 1/1-3000 for furture date
......@@ -57,8 +57,8 @@ if strcmp(station(1:3), 'VNA')
use6 = [3000, 1, 1, 0, 3000, 1, 2, 0];
end
if strcmp(station(1:3), 'WNG')
use6 = [2013, 7, 1, 0, 2018, 3, 15, 0 % needs finetuning
2018, 4, 6, 0, 3000, 1, 2, 0]; % needs finetuning
use6 = [2014, 5, 15, 0, 2014, 5, 21, 0 % needs finetuning
3000, 1, 1, 0, 3000, 1, 2, 0]; % needs finetuning
end
if strcmp(station(1:3), 'NGK')
use6 = [3000, 1, 1, 0, 3000, 1, 2, 0];
......@@ -103,18 +103,20 @@ end
path_cdf4 = [path2 bs station(1:3) '4' bs 'cdf' bs];
if strcmp(station(1:3), 'TDC') || strcmp(station(1:3), 'VNA')...
|| strcmp(station(1:3), 'BFO')
|| strcmp(station(1:3), 'BFO') || strcmp(station(1:3), 'WNG')
path_cdf4 = [path2 bs station(1:3) '1' bs 'cdf' bs];
elseif strcmp(station(1:3), 'NGK') || strcmp(station(1:3), 'WNG')...
|| strcmp(station(1:3), 'VSS') || strcmp(station(1:3), 'TTB')...
elseif strcmp(station(1:3), 'NGK') || ...
strcmp(station(1:3), 'VSS') || strcmp(station(1:3), 'TTB')...
|| strcmp(station(1:3), 'SHE')
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), 'WNG')...
if strcmp(station(1:3), 'NGK') ...
|| 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];
end
month_str = ['JAN'; 'FEB'; 'MAR'; 'APR'; 'MAY'; 'JUN';...
......@@ -158,13 +160,17 @@ for timeintervalmonth = 1:timeinterval(3)
min_file4 = [path_cdf4 station(1:3) '4_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
if strcmp(station(1:3), 'TDC') || strcmp(station(1:3), 'VNA')...
|| strcmp(station(1:3), 'BFO')
|| strcmp(station(1:3), 'BFO')|| strcmp(station(1:3), 'WNG')
min_file4 = [path_cdf4 station(1:3) '1_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
elseif strcmp(station(1:3), 'NGK') || strcmp(station(1:3), 'WNG')...
elseif strcmp(station(1:3), 'NGK') ...
|| strcmp(station(1:3), 'VSS') || strcmp(station(1:3), 'TTB')...
|| strcmp(station(1:3), 'SHE')
min_file4 = [path_cdf4 station(1:3) '0_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
end
......@@ -172,11 +178,13 @@ for timeintervalmonth = 1:timeinterval(3)
min_file6 = [path_cdf6 station(1:3) '6_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
if strcmp(station(1:3), 'NGK') || strcmp(station(1:3), 'WNG')...
if strcmp(station(1:3), 'NGK') ...
|| strcmp(station(1:3), 'TTB')
min_file6 = [path_cdf6 station(1:3) '1_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
elseif strcmp(station(1:3), 'WNG')
min_file6 = [path_cdf6 station(1:3) '0_' num2str(year, '%4.4d')...
num2str(month, '%2.2d') num2str(day, '%2.2d') '.mat'];
end
......
......@@ -55,14 +55,16 @@ clear all
% remember to check scale values!
% 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 = 'VNA1s', startime= jd2000(2013, 11, 01);
%station = 'TTB0g'; startime= jd2000(2015, 11, 01);
% station = 'TTB0g'; startime= jd2000(2015, 11, 01);
% station = 'TTB1l'; startime= jd2000(2008, 06, 01);
% station = 'NGK0g', startime= jd2000(2014, 1, 1);
% station = 'BFO1s', startime= jd2000(2014, 1, 1);
% station = 'VSS0g', startime= jd2000(2016, 1, 1);
%station = 'SMA0g', startime= jd2000(2016, 1, 1);
% station = 'SMA0g', startime= jd2000(2016, 1, 1);
% TIME SPAN,
%startime= jd2000(2014, 12, 03);
endtime = jd2000(2018, 12, 31);
......@@ -78,9 +80,9 @@ endscvtime = jd2000(2017, 9, 19);
% endscvtime = 5450,
% type of data to be produced
baselinetype = 'preliminary';
% baselinetype = 'preliminary';
% baselinetype = 'QD';
% baselinetype = 'definitive';
baselinetype = 'definitive';
%baselinetype = 'test';
write_IBFV200_for_year = 2017; %this I want to automize!!!!!!!!
......@@ -157,11 +159,11 @@ operator_str = [...
];
if strcmp(station(1:3), 'NGK')
bluecross = int8(30); % hier kann Zahl gendert werden
redcross = int8(10);
redcross = int8(42);
end %if
if strcmp(station(1:3), 'WNG')
bluecross = int8(16);
redcross = int8(10);
redcross = int8(42);
end %if
if strcmp(station(1:3), 'VNA')
bluecross = int8(10);
......@@ -282,12 +284,22 @@ for counteri = 1:size(D0)
end %for findnegi
%above replaces D0(D0<0) = NaN;
D0 = floor(D0/100) + (D0 - 100*floor(D0/100))/60;
% am 8.5.2018 geaendert und wieder rausgenommen, noch testen JMAT
% for counteri = 1:size(D0)
% if D0(counteri) >= 1
% D0(counteri) = floor(D0(counteri)/100)...
% + (D0(counteri) - 100*floor(D0(counteri)/100))/60;
% else
% D0(counteri) = (D0 - 100*floor(D0(counteri)/100))/60;
% end
% end
D0 = floor(D0/100) + (D0 - 100*floor(D0/100))/60; %wiederhergestellt
%new on 2014-01-09
for counteri = 1:size(D0)
if isnan(D0(counteri)) == 0 && D0(counteri) > 180
D0 = D0 - 360;
D0(counteri) = D0(counteri)-360;%am 8.5.2018 geaendert, testen JMAT
end %if
end %for findnegi
% above replaces: D0(D0 > 180) = D0 - 360;
......@@ -451,11 +463,11 @@ end
% triple-plot to magnify and get usefull numbers (not fully tested yet)
figure
subplot(3,1,1)
subplot(4,1,1)
plot(time, H0 - 10*floor(min(H0)/10), '.')
title(['add ' num2str(10*floor(min(H0)/10)) ' to get H0 in nT'])
title([station ', add ' num2str(10*floor(min(H0)/10)) ' to get H0 in nT'])
%figure
subplot(3,1,2)
subplot(4,1,2)
if min(D0) < 0
plot(time, D0 - 0.1*ceil(max(D0)/0.1), '.')
title(['add ' num2str(0.1*ceil(max(D0)/0.1)) ' to get D0 in degree'])
......@@ -464,7 +476,7 @@ else
title(['add ' num2str(0.1*floor(min(D0)/0.1)) ' to get D0 in degree'])
end
%figure
subplot(3,1,3)
subplot(4,1,3)
if min(Z0) < 0
plot(time, Z0 - 10*ceil(max(Z0)/10), '.')
title(['add ' num2str(10*ceil(max(Z0)/10)) ' to get Z0 in nT'])
......@@ -472,6 +484,9 @@ else
plot(time, Z0 - 10*floor(min(Z0)/10), '.')
title(['add ' num2str(10*floor(min(Z0)/10)) ' to get Z0 in nT'])
end
subplot(4,1,4)
plot(time, T1, '.')
title('Temperature T1')
xlabel('Julian days since 01.01.2000')
% three plots with adopted
......@@ -530,43 +545,46 @@ ylabel('Z0 in nT')
% plot(time, T2, '.')
% PLOT for scale vuale test
% save('WNG1_base_temp') Save workspace / absolute data here
scvplotHabsmeas = Habsmeas(time<endscvtime);
scvplotDabsmeas = Dabsmeas(time<endscvtime);
scvplotZabsmeas = Zabsmeas(time<endscvtime);
scvplotH0 = H0(time<endscvtime);
scvplotD0 = D0(time<endscvtime);
scvplotZ0 = Z0(time<endscvtime);
scvplottime = time(time<endscvtime);
scvplotHabsmeas = Habsmeas(time>startscvtime);
scvplotDabsmeas = Dabsmeas(time>startscvtime);
scvplotZabsmeas = Zabsmeas(time>startscvtime);
scvplotH0 = H0(time>startscvtime);
scvplotD0 = D0(time>startscvtime);
scvplotZ0 = Z0(time>startscvtime);
scvplotAzm = Azm(time<endscvtime);
scvplotAzm = Azm(time>startscvtime);
scvplottime = time(time>startscvtime);
% PLOT for scale vuale test
figure
plot(scvplotHabsmeas, scvplotH0, '+', 'MarkerEdgeColor','r')
title('H0 vs. H')
figure
plot(scvplotDabsmeas, scvplotD0, '+', 'MarkerEdgeColor','r')
title('D0 vs. D')
figure
plot(scvplotZabsmeas, scvplotZ0, '+', 'MarkerEdgeColor','r')
title('Z0 vs. Z')
figure
plot(scvplotDabsmeas, scvplotAzm, '+', 'MarkerEdgeColor','r')
title('D0 vs. Mire')
% scvplotHabsmeas = Habsmeas(time<endscvtime);
% scvplotDabsmeas = Dabsmeas(time<endscvtime);
% scvplotZabsmeas = Zabsmeas(time<endscvtime);
% scvplotH0 = H0(time<endscvtime);
% scvplotD0 = D0(time<endscvtime);
% scvplotZ0 = Z0(time<endscvtime);
% scvplottime = time(time<endscvtime);
% scvplotHabsmeas = Habsmeas(time>startscvtime);
% scvplotDabsmeas = Dabsmeas(time>startscvtime);
% scvplotZabsmeas = Zabsmeas(time>startscvtime);
% scvplotH0 = H0(time>startscvtime);
% scvplotD0 = D0(time>startscvtime);
% scvplotZ0 = Z0(time>startscvtime);
% scvplotAzm = Azm(time<endscvtime);
% scvplotAzm = Azm(time>startscvtime);
% scvplottime = time(time>startscvtime);
%
% figure
% plot(scvplotHabsmeas, scvplotH0, '+', 'MarkerEdgeColor','r')
% title('H0 vs. H')
%
% figure
% plot(scvplotDabsmeas, scvplotD0, '+', 'MarkerEdgeColor','r')
% title('D0 vs. D')
%
% figure
% plot(scvplotZabsmeas, scvplotZ0, '+', 'MarkerEdgeColor','r')
% title('Z0 vs. Z')
%
% figure
% plot(scvplotDabsmeas, scvplotAzm, '+', 'MarkerEdgeColor','r')
% title('D0 vs. Mire')
% figure
% plot(scvplotHabsmeas, scvplotHabsmeas-scvplotH0, '+', 'MarkerEdgeColor','r')
......@@ -625,6 +643,14 @@ for i = 1:size(c)
SCV_Xa = 1.0;
SCV_Ya = 1.0;
SCV_Za = 1.0;
elseif strcmp(station, 'WNG5g')
SCV_Xa = 1.0;
SCV_Ya = 1.0;
SCV_Za = 1.0;
elseif strcmp(station, 'WNG0g')
SCV_Xa = 1.0;
SCV_Ya = 1.0;
SCV_Za = 1.0;
elseif strcmp(station, 'WNG1g')
SCV_Xa = 1.0;
SCV_Ya = 1.0;
......
......@@ -56,6 +56,10 @@ function [H0alines, D0alines, Z0alines] = ...
% SCV_Xa = 1.0;
% SCV_Ya = 1.0;
% SCV_Za = 1.0;
% elseif strcmp(station, 'WNG5g')
% SCV_Xa = 1.0;
% SCV_Ya = 1.0;
% SCV_Za = 1.0;
% elseif strcmp(station, 'WNG1g')
% SCV_Xa = 1.0;
% SCV_Ya = 1.0;
......@@ -128,7 +132,54 @@ function [H0alines, D0alines, Z0alines] = ...
if strcmp(station, 'WNG0g')
%
% NOTE:
%
%
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'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 0.63400 'p' ''
5114 0.63400 'p' '2014-01-01'
5435.75 0.63400 'p' 'datum_Sensor_ausgetauscht'
5435.75 2.41600 'p' ''
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 45984.50 'p' ''
5114 45984.50 'p' '2014-01-01'
5435.75 45984.50 'p' 'datum_Sensor_ausgetauscht'
5435.75 46142.00 'p' ''
}; % REMEMBER for QD: default is the first day of the next month!!!
end
% ALT:
% if strcmp(station, 'WNG0g')
% %
% % NOTE:
% %
% %
% H0alines = {
% 5114 18092.00 'p' 'test01012014'
% 6600 18092.00 'p' 'test'
% }; % REMEMBER for QD: default is the first day of the next month!!!
% D0alines = {
% 5114 2.4167 'p' 'test'
% 6600 2.4167 'p' 'test'
% }; % REMEMBER for QD: default is the first day of the next month!!!
% Z0alines = {
% 5114 46143.00 'p' 'test'
% 6600 46143.00 'p' 'test'
% }; % REMEMBER for QD: default is the first day of the next month!!!
% end
......@@ -138,82 +189,55 @@ if strcmp(station, 'WNG1g')
%
%
H0alines = {
5113 18064.00 'p' ''
5114 18064.00 'p' '2014-01-01'
5875 18063.60 'q' '2016-02-01'
5904 18063.60 'q' '2016-03-01'
5965 18063.60 'q' '2016-05-01'
5996 18063.60 'q' '2016-06-01'
6026 18063.00 'q' '2016-07-01'
6057 18063.00 'q' '2016-08-01'
6088 18063.00 'q' '2016-09-01'
6118 18063.00 'q' '2016-10-01'
6149 18063.50 'q' '2016-11-01'