- Timestamp:
- Sep 23, 2020, 2:49:25 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/ephemeris.cpp
r9126 r9132 401 401 // Relativistic Correction 402 402 // ----------------------- 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;405 403 xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E); 406 // correspondent to IGS-SSR 404 407 405 xc[4] = _clock_drift + _clock_driftrate*tc; 408 406 xc[5] = _clock_driftrate; … … 641 639 acc[1] = _y_acceleration * 1.e3; 642 640 acc[2] = _z_acceleration * 1.e3; 643 641 644 642 for (int ii = 1; ii <= nSteps; ii++) { 645 643 _xv = rungeKutta4(_tt.gpssec(), _xv, step, acc, glo_deriv); … … 998 996 // Relativistic Correction 999 997 // ----------------------- 1000 // correspondent to Galileo ICD and to SSR standard1001 998 xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E); 1002 // correspondent to IGS convention1003 //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;1004 999 1005 1000 xc[4] = _clock_drift + _clock_driftrate*tc; … … 1660 1655 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 1661 1656 1662 //dotC = _clock_drift + _clock_driftrate*tc1663 // - 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; 1664 1659 1665 1660 // Relativistic Correction 1666 1661 // ----------------------- 1667 // correspondent to BDS ICD and to SSR standard1668 1662 xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E); 1669 // correspondent to IGS convention1670 // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;1671 1663 1672 1664 xc[4] = _clock_drift + _clock_driftrate*tc;
Note:
See TracChangeset
for help on using the changeset viewer.