obs251

IAGA2002_Kyoto_transfer_Katrin.m

@@ -21,8 +21,8 @@
 % attachment: same file
 
 clear all
-startdate = jd2000(2018, 05, 03);    % HERE START DATE FOR non-online
-enddate   = jd2000(2018, 05, 03);    % HERE END   DATE FOR non-online
+startdate = jd2000(2018, 05, 04);    % HERE START DATE FOR non-online
+enddate   = jd2000(2018, 05, 06);    % HERE END   DATE FOR non-online
 
 station = 'VNA1'; 
 %station = 'WNG1';

calculate_baselines.m

@@ -1763,6 +1763,10 @@ if exist(DI_file, 'file') == 2;
         H10 = 18868.09;    
         Z10 = 45418.48;
     end
+    if strcmp(station(1:4), 'WNG5') % passt gut
+        H10 = 18868.09;    
+        Z10 = 45418.48;
+    end
     
     if strcmp(station(1:4), 'NGK3') %  passt gut
         H10 = 18879.14;         

plot_baseline_universal.m

@@ -55,10 +55,10 @@ clear all
 % remember to check scale values!
 % TDC works only after 2011 01 01
 % station = 'TDC1d';  startime= jd2000(2011, 1, 1);
-% station = 'WNG1g',  startime= jd2000(2014, 1, 1);
+ station = 'WNG1g',  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 = '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);

ppmsec2cdf.m

@@ -936,8 +936,10 @@ for k = 1:size(Cfiles,1)
         F0(:,1) = 32;                % preliminary, only valid for 2010
         elseif idate < jd2000(2017, 05, 01)
         F0(:,1) = 32+156;            % preliminary  
+        elseif idate >= jd2000(2017, 05, 01)
+            F0(:,1) = 11;    
         else
-            F0(:,1) = 0;
+            error('keine Sockeldifferenz')
         end
     end
     jd2date(idate)

read_temperature.m

 
 clear all
 
-station = 'TDC1d';
-startdate = jd2000(2016, 1, 1);
-enddate   = jd2000(2016, 12, 31);
+station = 'WNG1g';
+startdate = jd2000(2018, 1, 1);
+enddate   = jd2000(2018, 4, 31);
 
 opsys = computer;
 if strcmp(opsys(1:2), 'GL')  == 1