quasi-real-time strongly improved (e.g., read-in and handling RT corr…

…ection streams); change to long stationnames and tropo files; improved CDDIS download

CODE/COMMON/DEF.m

@@ -29,13 +29,9 @@
         % maximum time difference for a correction from correction stream,
         % timely nearest correction from correction stream has to be under
         % this threshold
-        THRESHOLD_corr2brdc_clk_dt = 5;         % [s], for clock corrections
+        THRESHOLD_corr2brdc_clk_dt = 60;        % [s], for clock corrections
         THRESHOLD_corr2brdc_orb_dt = 120;       % [s], for orbit corrections
 
-        % time [s] which satellites is excluded if no broadcast ephemeris
-        % are found
-        EXCLUDE_TIME = 300;     % [s]
-
         % default observation type ranking
         RANKING_GPS = 'WCDPSLXYMNDIQF';
         RANKING_GLO = 'PCIQX';

CODE/COMMON/GeodeticToolbox/jd2gps_GT.m

@@ -28,4 +28,4 @@
 sow = (jd - (jdgps+nweek*7)) * 3600*24;
 rollover = fix(nweek/1024);  % rollover every 1024 weeks
 %gpsweek = mod(nweek,1024);
-gpsweek = nweek;
+gpsweek = nweek;

CODE/COMMON/setupRotation_LL2ECEF.m

@@ -3,6 +3,12 @@
 % Local Level origin in point with phi and lam
 % can be interpreted as the axes of the Local Level Frame expressed in the
 % ECEF frame
+% 
+% INPUT:
+%   phi         latitude [rad]
+%   lam         longitude [rad]
+% OUTPUT:
+%	R           3x3, rotation matrix from Local Level to ECEF frame
 %
 % This function belongs to raPPPid, Copyright (c) 2023, M.F. Glaner
 % *************************************************************************

CODE/FIXING/EW_fixing.m

@@ -74,8 +74,8 @@
 MW_gnss(:, remove) = [];
 
 MW_SD = MW_refSat - MW_gnss;                % collected MW LC single differenced to reference satellite
-%     std_MW = nanstd(MW_SD);               % stdev of collected MW LC, [cycles]
-mean_MW = nanmean(MW_SD);                   % mean of collected MW LC, [cycles]
+%     std_MW = std(MW_SD, 'omitnan');               % stdev of collected MW LC, [cycles]
+mean_MW = mean(MW_SD, 'omitnan');                   % mean of collected MW LC, [cycles]
 MW_round = round(mean_MW);                  % rounded mean of collected MW LC
 dist_round = abs(mean_MW - MW_round);   	% distance mean to rounded mean
 dist_MW = abs(mean_MW' - Epoch.WL_23(prns));	% distance to current EW fix

CODE/FIXING/HMW_fixing.m

@@ -69,8 +69,8 @@
 HMW_gnss(:, remove) = [];
 
 HMW_SD = HMW_refSat - HMW_gnss;             % collected HMW LC single differenced to reference satellite
-%     std_MW = nanstd(MW_SD);               % stdev of collected HMW LC, [cycles]
-mean_HMW = nanmean(HMW_SD);                 % mean of collected HMW LC, [cycles]
+%     std_MW = std(MW_SD, 'omitnan');               % stdev of collected HMW LC, [cycles]
+mean_HMW = mean(HMW_SD, 'omitnan');                 % mean of collected HMW LC, [cycles]
 HMW_round = round(mean_HMW);                % rounded mean of collected HMW LC
 dist_round = abs(mean_HMW - HMW_round);   	% distance mean to rounded mean
 dist_HMW = abs(mean_HMW' - fixed(prns));    % distance to current EW fix

CODE/FIXING/HMW_fixing_ZD.m

@@ -68,8 +68,8 @@
 prns(remove) = [];
 HMW_gnss(:, remove) = [];
 
-%     std_MW = nanstd(HMW_gnss);            % stdev of collected HMW LC, [cycles]
-mean_HMW = nanmean(HMW_gnss);               % mean of collected HMW LC, [cycles]
+%     std_MW = std(HMW_gnss, 'omitnan');            % stdev of collected HMW LC, [cycles]
+mean_HMW = mean(HMW_gnss, 'omitnan');               % mean of collected HMW LC, [cycles]
 HMW_round = round(mean_HMW);                % rounded mean of collected HMW LC
 dist_round = abs(mean_HMW - HMW_round);   	% distance mean to rounded mean
 dist_HMW = abs(mean_HMW' - fixed(prns));    % distance to current fix

CODE/FIXING/WL_fixing.m

@@ -98,8 +98,8 @@
 MW_gnss(:, remove) = [];
 
 MW_SD = MW_refSat - MW_gnss;                % collected MW LC single differenced to reference satellite
-%     std_MW = nanstd(MW_SD);               % stdev of collected MW LC, [cycles]
-mean_MW = nanmean(MW_SD);                   % mean of collected MW LC, [cycles]
+%     std_MW = std(MW_SD, 'omitnan');               % stdev of collected MW LC, [cycles]
+mean_MW = mean(MW_SD, 'omitnan');                   % mean of collected MW LC, [cycles]
 MW_round = round(mean_MW);                  % rounded mean of collected MW LC
 dist_round = abs(mean_MW - MW_round);   	% distance mean to rounded mean
 dist_MW = abs(mean_MW' - Epoch.WL_12(prns)); 	% distance to current WL fix

CODE/OBSERVATIONS/ObservationDownload/IGS_r3_stations.txt

@@ -2,7 +2,13 @@ BRUX00BEL
 AREG00PER
 MAW100ATA
 DJIG00DJI
-NYA200NOR
+ULAB00MNG
+NRC100CAN
+NNOR00AUS
+GAMB00PYF
+PNGM00PNG
+INEG00MEX
+HKSL00HKG
 
 
 % ABMF00GLP

CODE/OBSERVATIONS/apply_corr2brdc_biases.m

@@ -17,10 +17,14 @@
 % *************************************************************************
 
 
+age_biases = settings.ORBCLK.CorrectionStream_age(3);
+
 % CODE BIASES
 if settings.BIASES.code_corr2brdc_bool
     dt = Epoch.gps_time - input.ORBCLK.corr2brdc.t_code;
-    dt(dt < 0) = [];    % remove future data since to maintain real-time conditions
+    dt(dt < 0) = [];    % remove future data to maintain real-time conditions
+    dt(dt > age_biases) = [];   % remove corrections which are too old
+
     if ~isempty(dt)
         idx_corr2brdc = find(dt==min(dt));
         idx_corr2brdc = idx_corr2brdc(1);   % index of 1st timely nearest correction
@@ -40,7 +44,9 @@
 % PHASE BIASES
 if settings.BIASES.phase_corr2brdc_bool
     dt = Epoch.gps_time - input.ORBCLK.corr2brdc.t_code;
-    dt(dt < 0) = [];    % remove future data since to maintain real-time conditions
+    dt(dt < 0) = [];    % remove future data to maintain real-time conditions
+    dt(dt > age_biases) = [];   % remove corrections which are too old
+
     if ~isempty(dt)
         idx_corr2brdc = find(dt==min(dt));
         idx_corr2brdc = idx_corr2brdc(1);   % index of 1st timely nearest correction

CODE/OBSERVATIONS/checkMultipath.m

@@ -72,10 +72,10 @@
 % check if code time difference is above specified threshold. Thereby, subtract 
 % median (e.g., smartphone data) for each GNSS (e.g., different clock drifts)
 C1_diff_n_ = C1_diff_n;       % do keep dimension
-C1_diff_n_(Epoch.gps) = abs(C1_diff_n(Epoch.gps) - nanmedian(C1_diff_n(Epoch.gps)));   
-C1_diff_n_(Epoch.glo) = abs(C1_diff_n(Epoch.glo) - nanmedian(C1_diff_n(Epoch.glo)));  
-C1_diff_n_(Epoch.gal) = abs(C1_diff_n(Epoch.gal) - nanmedian(C1_diff_n(Epoch.gal)));  
-C1_diff_n_(Epoch.bds) = abs(C1_diff_n(Epoch.bds) - nanmedian(C1_diff_n(Epoch.bds))); 
+C1_diff_n_(Epoch.gps) = abs(C1_diff_n(Epoch.gps) - median(C1_diff_n(Epoch.gps), 'omitnan'));   
+C1_diff_n_(Epoch.glo) = abs(C1_diff_n(Epoch.glo) - median(C1_diff_n(Epoch.glo), 'omitnan'));  
+C1_diff_n_(Epoch.gal) = abs(C1_diff_n(Epoch.gal) - median(C1_diff_n(Epoch.gal), 'omitnan'));  
+C1_diff_n_(Epoch.bds) = abs(C1_diff_n(Epoch.bds) - median(C1_diff_n(Epoch.bds), 'omitnan')); 
 
 % check which code differences are above threshold and consider cooldown
 above_thresh = C1_diff_n_ > settings.OTHER.mp_thresh;

CODE/OBSERVATIONS/cycleSlip.m

@@ -50,9 +50,8 @@
     if settings.PROC.LLI
         Epoch = cycleSlip_LLI(Epoch, use_column, settings);
     end
-    
-    
+
+
     % save detected cycle slips in Epoch.sat_status
     Epoch.sat_status(any(Epoch.cs_found, 2)) = 3;
-
 end