Changeset 8619 in ntrip for trunk/BNC/src/PPP/pppSatObs.cpp
- Timestamp:
- Feb 28, 2019, 2:49:35 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/PPP/pppSatObs.cpp
r8495 r8619 62 62 // ------------------------------------------ 63 63 const string preferredAttrib = "CWPXI_"; 64 //const string preferredAttrib = "G:12&PWCSLXYN G:5&IQX R:12&PC R:3&IQX E:16&BCX E:578&IQX J:1&SLXCZ J:26&SLX J:5&IQX C:IQX I:ABCX S:1&C S:5&IQX"; 64 65 65 66 for (unsigned iFreq = 1; iFreq < t_frequency::max; iFreq++) { … … 106 107 // Compute Satellite Coordinates at Time of Transmission 107 108 // ----------------------------------------------------- 108 _xcSat.ReSize( 7); _xcSat = 0.0;109 _xcSat.ReSize(6); _xcSat = 0.0; 109 110 _vvSat.ReSize(3); _vvSat = 0.0; 110 111 bool totOK = false; 111 ColumnVector satPosOld( 7); satPosOld = 0.0;112 ColumnVector satPosOld(6); satPosOld = 0.0; 112 113 t_lc::type tLC = isValid(t_lc::cIF) ? t_lc::cIF : t_lc::c1; 113 114 double prange = obsValue(tLC); … … 328 329 _model._rho = rhoV.norm_Frobenius(); 329 330 331 ColumnVector vSat = _vvSat; 332 330 333 ColumnVector neu(3); 331 334 xyz2neu(station->ellApr().data(), rhoV.data(), neu.data()); … … 371 374 // ------------------ 372 375 _model._tropo = t_tropo::delay_saast(station->xyzApr(), _model._eleSat); 373 374 // Phase Wind-Up375 // -------------376 _model._windUp = station->windUp(_time, _prn, rSat);377 376 378 377 // Code Biases … … 393 392 // Phase Biases 394 393 // ----------- 395 // TODO: consideration of fix indicators , yaw angleand jump counter394 // TODO: consideration of fix indicators and jump counter 396 395 const t_satPhaseBias* satPhaseBias = PPP_CLIENT->obsPool()->satPhaseBias(_prn); 396 double yaw = 0.0; 397 bool ssr = false; 397 398 if (satPhaseBias) { 399 yaw = satPhaseBias->_yaw; 400 ssr = true; 398 401 for (unsigned ii = 0; ii < satPhaseBias->_bias.size(); ii++) { 399 402 const t_frqPhaseBias& bias = satPhaseBias->_bias[ii]; … … 407 410 } 408 411 412 // Phase Wind-Up 413 // ------------- 414 _model._windUp = station->windUp(_time, _prn, rSat, ssr, yaw, vSat) ; 415 416 409 417 // Tidal Correction 410 418 // ---------------- … … 428 436 } 429 437 } 438 439 // Relativistic effect due to earth gravity 440 // ---------------------------------------- 441 // TODO 430 442 431 443 // Ocean Loading
Note:
See TracChangeset
for help on using the changeset viewer.