Changeset 8483 in ntrip for trunk/BNC/src/RTCM/rtcm_utils.cpp

Timestamp:

Sep 17, 2018, 2:53:30 PM (6 years ago)

Author:

stuerze

Message:

SSR parameter clock rate, clock drift and URA are added within RTNET format

File:

• trunk/BNC/src/RTCM/rtcm_utils.cpp (modified) (7 diffs)

trunk/BNC/src/RTCM/rtcm_utils.cpp

@@ -10 +10 @@
 
 void resolveEpoch (double secsHour,
-                   int  refWeek,   double  refSecs,  
+                   int  refWeek,   double  refSecs,
                    int& epochWeek, double& epochSecs) {
 
-  const double secsPerWeek = 604800.0;                            
+  const double secsPerWeek = 604800.0;
 
   epochWeek = refWeek;
   epochSecs = secsHour + 3600.0*(floor((refSecs-secsHour)/3600.0+0.5));
-  
+
   if (epochSecs<0          ) { epochWeek--; epochSecs+=secsPerWeek; };
   if (epochSecs>secsPerWeek) { epochWeek++; epochSecs-=secsPerWeek; };
@@ -29 +29 @@
            double& xSat, double& ySat, double& zSat, double& clkSat) {
 
-  const double omega_earth = 7292115.1467e-11; 
-  const double secsPerWeek = 604800.0;                            
+  const double omega_earth = 7292115.1467e-11;
+  const double secsPerWeek = 604800.0;
 
   // Initial values
   // --------------
   rho = 0.0;
-  ColumnVector xc(4);
+  ColumnVector xc(7);
   ColumnVector vv(3);
   eph->getCrd(bncTime(GPSWeek, GPSWeeks), xc, vv, false);
@@ -45 +45 @@
   ////cout << "----- cmpRho -----\n";
   ////eph->print(cout);
-  ////cout << "  pos " << setw(4)  << GPSWeek 
+  ////cout << "  pos " << setw(4)  << GPSWeek
   ////     << " "      << setw(14) << setprecision(6) << GPSWeeks
   ////     << " "      << setw(13) << setprecision(3) << xSat
@@ -57 +57 @@
   do {
     rhoLast = rho;
-    
+
     // Correction station position due to Earth Rotation
     // -------------------------------------------------
     double dPhi = omega_earth * rho / c_light;
-    double xRec = stax * cos(dPhi) - stay * sin(dPhi); 
-    double yRec = stay * cos(dPhi) + stax * sin(dPhi); 
+    double xRec = stax * cos(dPhi) - stay * sin(dPhi);
+    double yRec = stay * cos(dPhi) + stax * sin(dPhi);
     double zRec = staz;
 
@@ -81 +81 @@
       GPSWeek_tot  += 1;
     }
-      
+
     eph->getCrd(bncTime(GPSWeek_tot, GPSWeeks_tot), xc, vv, false);
     xSat   = xc(1);
@@ -94 +94 @@
     rho = sqrt(dx*dx + dy*dy + dz*dz);
 
-    ////cout << "  scrd "   << setw(4)  << GPSWeek_tot 
+    ////cout << "  scrd "   << setw(4)  << GPSWeek_tot
     ////         << " "         << setw(15) << setprecision(8) << GPSWeeks_tot
     ////         << " "         << setw(13) << setprecision(3) << xSat
@@ -108 +108 @@
     ////         << " rho "     << setw(13) << setprecision(3) << rho
     ////         << endl;
-    
+
 
     ////cout.setf(ios::fixed);
     ////
-    ////cout << "niter " << setw(3) << ++niter 
+    ////cout << "niter " << setw(3) << ++niter
     ////         << " " << setw(14) << setprecision(3) << rhoLast
     ////         << " " << setw(14) << setprecision(3) << rho