Skip to content

GitLab

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

Daily work including wng5

parent bf903500
Hide whitespace changes
Inline Side-by-side
IAGA2002_Kyoto_transfer_Katrin.m
View file @ 6879234a
......@@ -21,17 +21,18 @@
% attachment: same file
clear all
startdate = jd2000(2018, 07, 10); % HERE START DATE FOR non-online
enddate = jd2000(2018, 07, 10); % HERE END DATE FOR non-online
startdate = jd2000(2018, 07, 18); % HERE START DATE FOR non-online
enddate = jd2000(2018, 07, 18); % HERE END DATE FOR non-online
station = 'VNA1';
%station = 'WNG1';
%station = 'BFO1';
%station = 'TDC1';
%station = 'NGK0';
%station = 'GAN0';
%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'};
......
calculate_baselines.m
View file @ 6879234a
......@@ -1195,8 +1195,46 @@ if exist(DI_file, 'file') == 2;
if pillar_nr == 1
% Sockel aus Altwingst.
error('Altwingst Mire und Pfeilerdiff noch hier kopieren')
deltaI = 0.0000 + 0.0000;
deltaD = 0.0000 + 0.0000;
%instrument difference and pillar difference
% Neue, zeitabhaengige Pfeilerdiff.
if timeDI(1) > jd2000(2019, 1, 1)
error('Keine Pfeilerdiff nach 1.1.2019')
elseif timeDI(1) > jd2000(2018, 1, 1)
deltappm = 0.0 + 0.0 % IST BEI WNG5 PER DEFINITION 0 nT
% instrument difference and pillar difference
else
error('Kein WNG deltappm')
end % Zeit
if mire_nr == 1
mire_azim = 12.93259 + ((0.2/60)/0.9) % addiere
% 0.2', damit Ergebnis konsistent ist zu N+NO, JM 17.7.2017
% in gon % from pillar 1 = NE-Pillar mire 1 = N-mire %11deg 38.36'
% =11deg 38' 22'' = 11.63933deg = 12.93259 gon
elseif mire_nr == 2
mire_azim = 14.87389 % in gon % from pillar 1 =NE-Pillar, mire 2 =NE-mire =13deg 23.19'
% =13deg 23' 19'' = 13.3865deg = 14,87389 gon % corrected
% on 20170714
elseif mire_nr == 3
mire_azim = 343.01741 % in gon %from pillar 1 =NE-pillar, mire 3 =W-mire = 308deg 42.94'
% = 308deg 42' 56'' = 308.71567 deg = 343.01741 gon
elseif mire_nr == 4 % mittelwert mire N und NE
mire_azim = 13.90324 % 12.93259 + 14.87389
elseif mire_nr == 5 %mittelwert aus W und N und NE
mire_azim = 123.60796 % 12.93259 + 14.87389 + 343.01741
else
error('This mire doesnt exist')
end % Mire
elseif pillar_nr == 3
......@@ -1246,6 +1284,8 @@ if exist(DI_file, 'file') == 2;
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')
......@@ -1332,6 +1372,26 @@ if strcmp(station(1:4), 'WNG0') || strcmp(station(1:4), 'WNG1')...
else
error('This mire doesnt exist')
end % Mire
elseif pillar_nr == 3 % das ist Neuwingst Absolutmessung mit Altwingst Variometer
% Hauptpfeiler in Neuwingst
display('Neuwingst Hauptpfeiler')
deltaI = 0.0000 + 0.0000;
deltaD = 0.0000 + 0.0000;
deltappm = 0.0 % SOLLTE PER DEFINITION RICHTIG SEIN!!!!
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
else
error('This pillar doesnt exist')
......
cdf_supplementing_from_SUP_GFZ.m 100644 → 100755
View file @ 6879234a
......@@ -39,7 +39,8 @@ clear all
close all
pack
% supply
station = 'NGK0g';
%station = 'NGK0g';
station = 'WNG1g';
%with
supplementer1 = 'file'; % supplementer ist stationsname z. Bsp. 'THL4d'
% oder 'file', dann wird die station aus den
......
check_IAF_vs_absolutes.m
View file @ 6879234a
......@@ -3,15 +3,15 @@ clear all
close all
%station = 'TDC1d';
station = 'NGK0g';
%station = 'GAN0j';
%station = 'NGK0g';
station = 'GAN0j';
%station = 'WNG1g';
%station = 'VNA1s';
%station = 'BFO1s';
%calculate IAF files for timeinterval:
timeinterval = [2017; 1 ; 18]; %start year, start month, number of months
timeinterval = [2017; 1 ; 16]; %start year, start month, number of months
% option to apply time correction
time_correction_str = 'n';
......
def2iaf.m
View file @ 6879234a
......@@ -20,9 +20,9 @@ clear all
IAFversion = '2.10';
%station = 'GDH6d';
%station = 'WNG1';
% station = 'GAN0';
station = 'GAN0';
%station = 'TDC1';
station = 'NGK0';
%station = 'NGK0';
%station = 'SHE0';
%station = 'VNA1';
%station = 'BFO1';
......@@ -30,7 +30,7 @@ station = 'NGK0';
%calculate IAF files for timeinterval:
timeinterval = [2018; 6; 1]; %start year, start month, number of months
timeinterval = [2018; 1; 4]; %start year, start month, number of months
% works for 2011 data onwards--
......
plot_baseline_universal.m
View file @ 6879234a
......@@ -55,14 +55,14 @@ 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 = '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);
% station = 'NGK0g', startime= jd2000(2014, 1, 1);
% station = 'GAN0j', startime= jd2000(2014, 1, 1);
station = 'GAN0j', startime= jd2000(2014, 1, 1);
% station = 'BFO1s', startime= jd2000(2014, 1, 1);
% station = 'VSS1g', startime= jd2000(2016, 1, 1);
% station = 'SMA0g', startime= jd2000(2016, 1, 1);
......
plot_baseline_universal_adopted.m
View file @ 6879234a
......@@ -150,6 +150,9 @@ if strcmp(station, 'GAN0j')
6420 37668.3 'd' ''
6575 37667.5 'd' '2018-01-01'
6606 37667.3 'd' '2018-02-01'
6634 37666.8 'd' '2018-03-01'
6665 37665.8 'd' '2018-04-01'
6695 37664.8 'd' '2018-05-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 -4.50600 'p' 'very prelim'
......@@ -172,8 +175,10 @@ if strcmp(station, 'GAN0j')
6450 -4.94600 'd' ''
6500 -4.94450 'd' ''
6575 -4.94450 'd' '2018-01-01'
6606 -4.94460 'd' '2018-02-01'
6606 -4.94300 'd' '2018-02-01'
6634 -4.94490 'd' '2018-03-01'
6665 -4.94890 'd' '2018-04-01'
6695 -4.95100 'd' '2018-05-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 -13329.50 'p' 'very prelim'
......@@ -196,7 +201,10 @@ if strcmp(station, 'GAN0j')
6480 -13326.10 'p' ''
6550 -13326.10 'p' ''
6575 -13325.00 'p' '2018-01-01 very prelim'
6606 -13325.00 'p' '2018-02-01'
6606 -13323.00 'p' '2018-02-01'
6634 -13323.00 'p' '2018-03-01'
6665 -13326.00 'p' '2018-04-01'
6695 -13327.00 'p' '2018-05-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......@@ -283,6 +291,7 @@ if strcmp(station, 'WNG1g')
6210 18064.00 'd' '2017-01-01'
6575 18064.00 'd' '2018-01-01'
6665 18065.30 'q' '2018-04-01'
6695 18065.30 'q' '2018-05-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
D0alines = {
5113 0.684000 'd' ''
......@@ -297,7 +306,8 @@ if strcmp(station, 'WNG1g')
6460 0.631000 'd' ''
6483 0.630000 'd' ''
6575 0.630000 'd' '2018-01-01'
6665 0.63000 'q' '2018-04-01'
6665 0.630000 'q' '2018-04-01'
6695 0.628000 'q' '2018-05-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
Z0alines = {
5113 45847.00 'd' ''
......@@ -321,7 +331,8 @@ if strcmp(station, 'WNG1g')
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'
6665 45846.20 'q' '2018-04-01'
6695 45845.50 'q' '2018-05-01'
}; % REMEMBER for QD: default is the first day of the next month!!!
end
......
read_temperature_10days.m 100644 → 100755
View file @ 6879234a
......@@ -9,10 +9,10 @@ close all
wievieleTage = 20;
%station_array = ['WNG1'; 'WNG0'; 'NGK1'; 'NGK0'; 'VNA1'; 'TDC1'; 'BFO1';'SMA0'];
station_array = ['WNG1'; 'WNG0'; 'WNG5';'NGK1'; 'NGK0'; 'VNA1';'VSS0'; 'TDC1'; 'BFO1';'SMA0';'GAN0'];
%station_array = ['WNG1'; 'WNG0';];
%station_array = ['NGK6'; 'TTB0';]
station_array = ['TDC1';]
%station_array = ['TDC1';]
%station_array = ['SMA0';]
[year1, month1, day1] = datevec(date)
......
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