Index: branches/BNC_2.12/src/ephemeris.cpp
===================================================================
--- branches/BNC_2.12/src/ephemeris.cpp	(revision 9130)
+++ branches/BNC_2.12/src/ephemeris.cpp	(revision 9133)
@@ -404,8 +404,6 @@
   // Relativistic Correction
   // -----------------------
-  // correspondent to IGS convention and GPS ICD (and SSR standard)
-  // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
   xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
-  // correspondent to IGS-SSR
+
   xc[4] = _clock_drift + _clock_driftrate*tc;
   xc[5] = _clock_driftrate;
@@ -429,5 +427,5 @@
   QString fmt = version < 3.0 ? "   %1%2%3%4\n" : "    %1%2%3%4\n";
 
-   out << QString(fmt)
+  out << QString(fmt)
     .arg(_IODE,    19, 'e', 12)
     .arg(_Crs,     19, 'e', 12)
@@ -644,4 +642,5 @@
   acc[1] = _y_acceleration * 1.e3;
   acc[2] = _z_acceleration * 1.e3;
+
   for (int ii = 1; ii <= nSteps; ii++) {
     _xv = rungeKutta4(_tt.gpssec(), _xv, step, acc, glo_deriv);
@@ -1000,8 +999,5 @@
   // Relativistic Correction
   // -----------------------
-  // correspondent to Galileo ICD and to SSR standard
   xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
-  // correspondent to IGS convention
-  //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
 
   xc[4] = _clock_drift + _clock_driftrate*tc;
@@ -1662,13 +1658,10 @@
   xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
 
-  // dotC  = _clock_drift + _clock_driftrate*tc
-  //       - 4.442807633e-10*_e*sqrt(a0)*cos(E) * dEdM * n;
+//  dotC  = _clock_drift + _clock_driftrate*tc
+//        - 4.442807633e-10*_e * sqrt(a0) * cos(E) * dEdM * n;
 
   // Relativistic Correction
   // -----------------------
-  // correspondent to BDS ICD and to SSR standard
   xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
-  // correspondent to IGS convention
-  // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
 
   xc[4] = _clock_drift + _clock_driftrate*tc;
