Changeset 10034 in ntrip for trunk/BNC/src/PPP/pppSatObs.cpp
- Timestamp:
- Apr 21, 2023, 11:48:24 AM (12 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/PPP/pppSatObs.cpp
r10014 r10034 44 44 _valid = true; 45 45 _reference = false; 46 _stecRefSat = 0.0;47 46 _stecSat = 0.0; 48 47 _signalPriorities = QString::fromStdString(OPT->_signalPriorities); 49 48 if (!_signalPriorities.size()) { 50 if (OPT->_obsModelType == OPT->DCMcodeBias || 51 OPT->_obsModelType == OPT->DCMphaseBias) { 52 // at the moment only one code or phase bias per system (G,R,E,C)/modulation considered 53 _signalPriorities = "G:12&CW R:12&CP E:1&CX E:5&QX C:26&I"; 54 } 55 else { 56 _signalPriorities = "G:12&CWPSLX R:12&CP E:1&CBX E:5&QIX C:26&IQX"; 57 } 49 _signalPriorities = "G:12&CWPSLX R:12&CP E:1&CBX E:5&QIX C:26&IQX"; 58 50 } 59 51 … … 109 101 if (obs->_rnxType2ch == obsType.toStdString() && 110 102 obs->_codeValid && obs->_code && 111 obs->_phaseValid && obs->_phase) { 103 obs->_phaseValid && obs->_phase && 104 obs->_lockTimeValid && obs->_lockTime > 5.0) { 112 105 _obs[iFreq] = new t_frqObs(*obs); //cout << "================> newObs: " << obs->_rnxType2ch <<endl; 113 106 } … … 261 254 // Pseudo observations 262 255 if (tLC == t_lc::GIM) { 263 if (_stec RefSat == 0.0 || _stecSat == 0.0) {256 if (_stecSat == 0.0) { 264 257 if (valid) *valid = false; 265 258 return 0.0; 266 259 } 267 260 else { 268 return _stec RefSat;261 return _stecSat; 269 262 } 270 263 } … … 339 332 340 333 if (tLC == t_lc::GIM) { 341 retVal = OPT->_sigmaGIM * OPT->_sigmaGIM + OPT->_sigmaGIM * OPT->_sigmaGIM;334 retVal = OPT->_sigmaGIM * OPT->_sigmaGIM; 342 335 } 343 336 … … 587 580 double f1GPS = t_CST::freq(t_frequency::G1, 0); 588 581 for (unsigned iFreq = 1; iFreq < t_frequency::max; iFreq++) { 589 if (OPT->_pseudoObsIono) { // DCMcodeBias, DCMphaseBias582 if (OPT->_pseudoObsIono) { 590 583 // For scaling the slant ionospheric delays the trick is to be consistent with units! 591 584 // The conversion of TECU into meters requires the frequency of the signal. … … 737 730 // 738 731 //////////////////////////////////////////////////////////////////////////// 739 void t_pppSatObs::setPseudoObsIono(t_frequency::type freq, double stecRefSat) { 732 bool t_pppSatObs::setPseudoObsIono(t_frequency::type freq) { 733 bool pseudoObsIono = false; 740 734 _stecSat = _model._ionoCodeDelay[freq]; 741 _stecRefSat = stecRefSat; 742 } 735 if (_stecSat) { 736 pseudoObsIono = true; 737 } 738 return pseudoObsIono; 739 }
Note:
See TracChangeset
for help on using the changeset viewer.