Changeset 9132 in ntrip


Ignore:
Timestamp:
Sep 23, 2020, 2:49:25 PM (5 weeks ago)
Author:
stuerze
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/ephemeris.cpp

    r9126 r9132  
    401401  // Relativistic Correction
    402402  // -----------------------
    403   // correspondent to IGS convention and GPS ICD (and SSR standard)
    404   // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    405403  xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
    406   // correspondent to IGS-SSR
     404
    407405  xc[4] = _clock_drift + _clock_driftrate*tc;
    408406  xc[5] = _clock_driftrate;
     
    641639  acc[1] = _y_acceleration * 1.e3;
    642640  acc[2] = _z_acceleration * 1.e3;
    643  
     641
    644642  for (int ii = 1; ii <= nSteps; ii++) {
    645643    _xv = rungeKutta4(_tt.gpssec(), _xv, step, acc, glo_deriv);
     
    998996  // Relativistic Correction
    999997  // -----------------------
    1000   // correspondent to Galileo ICD and to SSR standard
    1001998  xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
    1002   // correspondent to IGS convention
    1003   //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    1004999
    10051000  xc[4] = _clock_drift + _clock_driftrate*tc;
     
    16601655  xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
    16611656
    1662   // dotC  = _clock_drift + _clock_driftrate*tc
    1663   //       - 4.442807633e-10*_e*sqrt(a0)*cos(E) * dEdM * n;
     1657// dotC  = _clock_drift + _clock_driftrate*tc
     1658//        - 4.442807633e-10*_e * sqrt(a0) * cos(E) * dEdM * n;
    16641659
    16651660  // Relativistic Correction
    16661661  // -----------------------
    1667   // correspondent to BDS ICD and to SSR standard
    16681662  xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
    1669   // correspondent to IGS convention
    1670   // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    16711663
    16721664  xc[4] = _clock_drift + _clock_driftrate*tc;
Note: See TracChangeset for help on using the changeset viewer.