Changeset 9133 in ntrip


Ignore:
Timestamp:
Sep 23, 2020, 2:56:18 PM (4 weeks ago)
Author:
stuerze
Message:

cleanup of the relativistic effects

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/BNC_2.12/src/ephemeris.cpp

    r9127 r9133  
    404404  // Relativistic Correction
    405405  // -----------------------
    406   // correspondent to IGS convention and GPS ICD (and SSR standard)
    407   // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    408406  xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
    409   // correspondent to IGS-SSR
     407
    410408  xc[4] = _clock_drift + _clock_driftrate*tc;
    411409  xc[5] = _clock_driftrate;
     
    429427  QString fmt = version < 3.0 ? "   %1%2%3%4\n" : "    %1%2%3%4\n";
    430428
    431    out << QString(fmt)
     429  out << QString(fmt)
    432430    .arg(_IODE,    19, 'e', 12)
    433431    .arg(_Crs,     19, 'e', 12)
     
    644642  acc[1] = _y_acceleration * 1.e3;
    645643  acc[2] = _z_acceleration * 1.e3;
     644
    646645  for (int ii = 1; ii <= nSteps; ii++) {
    647646    _xv = rungeKutta4(_tt.gpssec(), _xv, step, acc, glo_deriv);
     
    1000999  // Relativistic Correction
    10011000  // -----------------------
    1002   // correspondent to Galileo ICD and to SSR standard
    10031001  xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
    1004   // correspondent to IGS convention
    1005   //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    10061002
    10071003  xc[4] = _clock_drift + _clock_driftrate*tc;
     
    16621658  xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
    16631659
    1664   // dotC  = _clock_drift + _clock_driftrate*tc
    1665   //       - 4.442807633e-10*_e*sqrt(a0)*cos(E) * dEdM * n;
     1660// dotC  = _clock_drift + _clock_driftrate*tc
     1661//        - 4.442807633e-10*_e * sqrt(a0) * cos(E) * dEdM * n;
    16661662
    16671663  // Relativistic Correction
    16681664  // -----------------------
    1669   // correspondent to BDS ICD and to SSR standard
    16701665  xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
    1671   // correspondent to IGS convention
    1672   // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    16731666
    16741667  xc[4] = _clock_drift + _clock_driftrate*tc;
Note: See TracChangeset for help on using the changeset viewer.