 Timestamp:
 Sep 23, 2020, 2:56:18 PM (8 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/BNC_2.12/src/ephemeris.cpp
r9127 r9133 404 404 // Relativistic Correction 405 405 //  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;408 406 xc[3] = 4.442807633e10 * _e * sqrt(a0) *sin(E); 409 // correspondent to IGSSSR 407 410 408 xc[4] = _clock_drift + _clock_driftrate*tc; 411 409 xc[5] = _clock_driftrate; … … 429 427 QString fmt = version < 3.0 ? " %1%2%3%4\n" : " %1%2%3%4\n"; 430 428 431 429 out << QString(fmt) 432 430 .arg(_IODE, 19, 'e', 12) 433 431 .arg(_Crs, 19, 'e', 12) … … 644 642 acc[1] = _y_acceleration * 1.e3; 645 643 acc[2] = _z_acceleration * 1.e3; 644 646 645 for (int ii = 1; ii <= nSteps; ii++) { 647 646 _xv = rungeKutta4(_tt.gpssec(), _xv, step, acc, glo_deriv); … … 1000 999 // Relativistic Correction 1001 1000 //  1002 // correspondent to Galileo ICD and to SSR standard1003 1001 xc[3] = 4.442807633e10 * _e * sqrt(a0) *sin(E); 1004 // correspondent to IGS convention1005 //xc[3] = 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;1006 1002 1007 1003 xc[4] = _clock_drift + _clock_driftrate*tc; … … 1662 1658 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 1663 1659 1664 //dotC = _clock_drift + _clock_driftrate*tc1665 //  4.442807633e10*_e*sqrt(a0)*cos(E) * dEdM * n;1660 // dotC = _clock_drift + _clock_driftrate*tc 1661 //  4.442807633e10*_e * sqrt(a0) * cos(E) * dEdM * n; 1666 1662 1667 1663 // Relativistic Correction 1668 1664 //  1669 // correspondent to BDS ICD and to SSR standard1670 1665 xc[3] = 4.442807633e10 * _e * sqrt(a0) *sin(E); 1671 // correspondent to IGS convention1672 // xc[3] = 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;1673 1666 1674 1667 xc[4] = _clock_drift + _clock_driftrate*tc;
Note: See TracChangeset
for help on using the changeset viewer.