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

Neuwingst mira, daily stuff

parent f484095b
......@@ -21,8 +21,8 @@
% attachment: same file
clear all
startdate = jd2000(2018, 06, 29); % HERE START DATE FOR non-online
enddate = jd2000(2018, 07, 01); % HERE END DATE FOR non-online
startdate = jd2000(2018, 07, 10); % HERE START DATE FOR non-online
enddate = jd2000(2018, 07, 10); % HERE END DATE FOR non-online
station = 'VNA1';
%station = 'WNG1';
......@@ -30,8 +30,8 @@ station = 'VNA1';
%station = 'TDC1';
%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'};
......
......@@ -1180,13 +1180,13 @@ if exist(DI_file, 'file') == 2;
if strcmp(station(1:4), 'WNG5') || strcmp(station(1:4), 'WNG6')
% Das ist Neuwingst.
% Wir wollen spaeter auch die Variometer in Neu-Wingst mit
% den Absolutmessungen in Altwingst und umgekehrt vergleichen,
% deshalb bekommt Neu-Wingst die Sockelnummer 3 (1 ist
% Altwingst und 2 ist nicht belegt, 8 wird fuer Testbetrieb in
% Niemegk benutzt.
pillar_nr = raw1{13,6};
mire_nr = raw1{8,6}
if timeDI(1) < jd2000(2018, 5, 1, 0)
......@@ -1194,24 +1194,42 @@ if exist(DI_file, 'file') == 2;
end
if pillar_nr == 1
% Sockel aus Altwingst.
error('Altwingst Mire und Pfeilerdiff noch hier kopieren')
elseif pillar_nr == 3
% Hauptpfeiler in Neuwingst
display('Neuwingst Hauptpfeiler')
deltaI = 0.0000 + 0.0000;
deltaD = 0.0000 + 0.0000;
deltappm = 0.0 % muss noch angepasst weden!!!!!!!!!!!!!!
if mire_nr == 1
display('Neuwingst Haupt-Mire')
mire_azim = 0.0 %gon, noch anpassen!!!!!!!!!!!!!!!!!!!
elseif mire_nr == 2
display('Neuwingst Variometer-Haube-Mire')
mire_azim = 20.0 %gon, noch anpassen!!!!!!!!!!!!!!!!!!!
elseif mire_nr == 3
display('Neuwingst Elektronikhuette-Mire')
mire_azim = 200.0 %gon, noch anpassen!!!!!!!!!!!!!!!!!!!
else
error('This mire doesnt exist')
end % Mire
elseif pillar_nr == 8 % NUR FR TESTBETRIEB IN NIEMEGK
% Niemegk Hauptpfeiler
display('niemegk sockel 8')
deltaI = 0.0000 + 0.0000;
deltaD = 0.0000 + 0.0000;
deltappm = 22.93 % EIN GSM bei Testhuette in Niemegk
if mire_nr == 1
mire_azim = 101.30988 % gon, water tower WT from pillar 8
......@@ -1224,11 +1242,10 @@ if exist(DI_file, 'file') == 2;
end % Mire
else
error('This pillar doesnt exist')
end %Pillar number
end % Neuwingst
end %Pillar number
end % Neuwingst
if strcmp(station(1:4), 'WNG0') || strcmp(station(1:4), 'WNG1')...
|| strcmp(station(1:4), 'WNG2')...
|| strcmp(station(1:4), 'WNG3')
......
......@@ -39,9 +39,9 @@
clear all
pack
% supply
station = 'NGK0g'; % HIER IST DIE DATENLCKE
station = 'TTB1g'; % HIER IST DIE DATENLCKE
%with
supplementer = 'NGK6g'; % VON HIER DATEN HOLEN
supplementer = 'TTB1l'; % 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 = 'NGK6g'; % VON HIER DATEN HOLEN
% reading of SUP files!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! COULD BE
% VERY EASY
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
start_supp = jd2000(2018, 05, 26, 00 + 00 / 60) %Stunde Minute/60 INKLUSIVE
end___supp = jd2000(2018, 05, 26, 12 + 27 / 60) %Stunde Minute/60 EXKLUSIVE
% chooses data investigation interval to determine linear relationship
% or user select a start and end data
......
......@@ -11,7 +11,7 @@ station = 'NGK0g';
%calculate IAF files for timeinterval:
timeinterval = [2017; 1 ; 17]; %start year, start month, number of months
timeinterval = [2017; 1 ; 18]; %start year, start month, number of months
% option to apply time correction
time_correction_str = 'n';
......
......@@ -30,7 +30,7 @@ station = 'NGK0';
%calculate IAF files for timeinterval:
timeinterval = [2018; 5; 1]; %start year, start month, number of months
timeinterval = [2018; 6; 1]; %start year, start month, number of months
% works for 2011 data onwards--
......@@ -60,7 +60,8 @@ 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
2018, 3, 16, 0, 2018, 4, 5, 0
2018, 6, 8, 0, 2018, 7, 2, 0
3000, 1, 1, 0, 3000, 1, 2, 0];%
end
if strcmp(station(1:3), 'NGK') % Hauptsystem in NGK0, use6 ist NGK1
......
......@@ -55,11 +55,11 @@ 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 = '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 = 'GAN0j', startime= jd2000(2014, 1, 1);
......@@ -70,7 +70,7 @@ clear all
% TIME SPAN,
%startime= jd2000(2014, 12, 03);
endtime = jd2000(2018, 05, 31);
endtime = jd2000(2018, 05, 25);
endtime = jd2000(2018, 12, 31);
......
......@@ -469,6 +469,7 @@ if strcmp(station, 'NGK0g')
6665 18929.5 'q' '2018-04-01'
6695 18929.5 'q' '2018-05-01'
6726 18929.5 'q' '2018-06-01'
6756 18929.5 'q' '2018-07-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 2.45166 'p' ''
......@@ -501,6 +502,7 @@ if strcmp(station, 'NGK0g')
6665 2.44700 'q' '2018-04-01'
6695 2.44900 'q' '2018-05-01'
6726 2.44900 'q' '2018-06-01'
6756 2.44900 'q' '2018-07-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 45426.17 'p' ''
......@@ -536,6 +538,7 @@ if strcmp(station, 'NGK0g')
6665 45425.90 'q' '2018-04-01'
6695 45425.80 'q' '2018-05-01'
6726 45425.50 'q' '2018-06-01'
6756 45425.50 'q' '2018-07-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......
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