Changeset 8541 in ntrip for branches/BNC_2.12/src/ephemeris.cpp
- Timestamp:
- Nov 15, 2018, 11:14:49 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/BNC_2.12/src/ephemeris.cpp
r8497 r8541 65 65 << 900 << 1800 << 3600 << 7200 66 66 << 10800; 67 xc.ReSize( 7);67 xc.ReSize(6); 68 68 vv.ReSize(3); 69 69 if (position(tt.gpsw(), tt.gpssec(), xc.data(), vv.data()) != success) { … … 320 320 static const double gmGRS = 398.6005e12; 321 321 322 memset(xc, 0, 7*sizeof(double));322 memset(xc, 0, 6*sizeof(double)); 323 323 memset(vv, 0, 3*sizeof(double)); 324 324 … … 368 368 double tc = tt - _TOC; 369 369 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 370 371 xc[4] = _clock_bias;372 xc[5] = _clock_drift;373 xc[6] = _clock_driftrate;374 370 375 371 // Velocity … … 402 398 xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c; 403 399 400 xc[4] = _clock_drift + _clock_driftrate*tc; 401 xc[5] = _clock_driftrate; 402 404 403 return success; 405 404 } … … 618 617 static const double nominalStep = 10.0; 619 618 620 memset(xc, 0, 7*sizeof(double));619 memset(xc, 0, 6*sizeof(double)); 621 620 memset(vv, 0, 3*sizeof(double)); 622 621 … … 655 654 xc[3] = -_tau + _gamma * dtClk; 656 655 657 xc[4] = -_tau; 658 xc[5] = _gamma; 659 xc[6] = 0.0; 656 xc[4] = _gamma; 657 xc[5] = 0.0; 660 658 661 659 return success; … … 922 920 static const double gmWGS = 398.6004418e12; 923 921 924 memset(xc, 0, 7*sizeof(double));922 memset(xc, 0, 6*sizeof(double)); 925 923 memset(vv, 0, 3*sizeof(double)); 926 924 … … 970 968 double tc = tt - _TOC; 971 969 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 972 973 xc[4] = _clock_bias;974 xc[5] = _clock_drift;975 xc[6] = _clock_driftrate;976 970 977 971 // Velocity … … 1006 1000 //xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c; 1007 1001 1002 xc[4] = _clock_drift + _clock_driftrate*tc; 1003 xc[5] = _clock_driftrate; 1004 1008 1005 return success; 1009 1006 } … … 1317 1314 xc[3] = _agf0 + _agf1 * dt; 1318 1315 1319 xc[4] = _agf0; 1320 xc[5] = _agf1; 1321 xc[6] = 0.0; 1316 xc[4] = _agf1; 1317 xc[5] = 0.0; 1322 1318 1323 1319 return success; … … 1698 1694 xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc; 1699 1695 1700 xc[4] = _clock_bias;1701 xc[5] = _clock_drift;1702 xc[6] = _clock_driftrate;1703 1704 1696 // dotC = _clock_drift + _clock_driftrate*tc 1705 1697 // - 4.442807633e-10*_e*sqrt(a0)*cos(E) * dEdM * n; … … 1708 1700 // ----------------------- 1709 1701 // correspondent to BDS ICD and to SSR standard 1710 1702 xc[3] -= 4.442807633e-10 * _e * sqrt(a0) *sin(E); 1711 1703 // correspondent to IGS convention 1712 1704 // xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c; 1705 1706 xc[4] = _clock_drift + _clock_driftrate*tc; 1707 xc[5] = _clock_driftrate; 1713 1708 1714 1709 return success;
Note:
See TracChangeset
for help on using the changeset viewer.