Changeset 10023 in ntrip for trunk/BNC/src/PPP/pppParlist.cpp


Ignore:
Timestamp:
Apr 3, 2023, 4:02:31 PM (13 months ago)
Author:
stuerze
Message:

change from receiver clock offsets to system clocks

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/BNC/src/PPP/pppParlist.cpp

    r10018 r10023  
    6262     _noise   = OPT->_noiseCrd[2];
    6363     break;
    64    case clkR:
     64   case rClkG:
     65     _epoSpec = true;
     66     _sigma0  = OPT->_aprSigClk;
     67     break;
     68   case rClkR:
     69     _epoSpec = true;
     70     _sigma0  = OPT->_aprSigClk;
     71     break;
     72   case rClkE:
     73     _epoSpec = true;
     74     _sigma0  = OPT->_aprSigClk;
     75     break;
     76   case rClkC:
    6577     _epoSpec = true;
    6678     _sigma0  = OPT->_aprSigClk;
    … …  
    7486         const t_pppSatObs* obs = obsVector->at(ii);
    7587         if (obs->prn() == _prn) {
    76            double offGlo = 0;
    77            if (_prn.system() == 'R' && tLC != t_lc::MW) {
    78              offGlo = PPP_CLIENT->offGlo();
    79            }
    80            double offGal = 0;
    81            if (_prn.system() == 'E' && tLC != t_lc::MW) {
    82              offGal = PPP_CLIENT->offGal();
    83            }
    84            double offBds = 0;
    85            if (_prn.system() == 'C' && tLC != t_lc::MW) {
    86              offBds = PPP_CLIENT->offBds();
    87            }
    88            _x0 = floor((obs->obsValue(tLC) - offGlo - offGal - offBds - obs->cmpValue(tLC)) / obs->lambda(tLC) + 0.5);
     88           _x0 = floor((obs->obsValue(tLC) - obs->cmpValue(tLC)) / obs->lambda(tLC) + 0.5);
    8989           break;
    9090         }
    9191       }
    9292     }
    93      break;
    94    case offGlo:
    95      _epoSpec = true;
    96      _sigma0  = OPT->_aprSigClk;
    97      _x0      = PPP_CLIENT->offGlo();
    98      break;
    99    case offGal:
    100      _epoSpec = true;
    101      _sigma0  = OPT->_aprSigClk;
    102      _x0      = PPP_CLIENT->offGal();
    103      break;
    104    case offBds:
    105      _epoSpec = true;
    106      _sigma0  = OPT->_aprSigClk;
    107      _x0      = PPP_CLIENT->offBds();
    10893     break;
    10994   case trp:
    … …  
    191176    if (tLC == t_lc::GIM) {return 0.0;}
    192177    return (sta->xyzApr()[2] - obs->xc()[2]) / rhoV.NormFrobenius();
    193   case clkR:
     178  case rClkG:
    194179    if (tLC == t_lc::GIM) {return 0.0;}
    195     return 1.0;
    196   case offGlo:
     180    return (obs->prn().system() == 'G') ? 1.0 : 0.0;
     181  case rClkR:
    197182    if (tLC == t_lc::GIM) {return 0.0;}
    198183    return (obs->prn().system() == 'R') ? 1.0 : 0.0;
    199   case offGal:
     184  case rClkE:
    200185    if (tLC == t_lc::GIM) {return 0.0;}
    201186    return (obs->prn().system() == 'E') ? 1.0 : 0.0;
    202   case offBds:
     187  case rClkC:
    203188    if (tLC == t_lc::GIM) {return 0.0;}
    204189    return (obs->prn().system() == 'C') ? 1.0 : 0.0;
    … …  
    325310    ss << "CRD_Z";
    326311    break;
    327   case clkR:
    328     ss << "REC_CLK     ";
    329     break;
    330   case offGlo:
    331     ss << "OFF_GLO     ";
    332     break;
    333   case offGal:
    334     ss << "OFF_GAL     ";
    335     break;
    336   case offBds:
    337     ss << "OFF_BDS     ";
     312  case rClkG:
     313    ss << "REC_CLK  G  ";
     314    break;
     315  case rClkR:
     316    ss << "REC_CLK  R  ";
     317    break;
     318  case rClkE:
     319    ss << "REC_CLK  E  ";
     320    break;
     321  case rClkC:
     322    ss << "REC_CLK  C  ";
    338323    break;
    339324  case trp:
    … …  
    688673  }
    689674
    690   // Receiver Clock
    691   // --------------
    692   required.push_back(new t_pppParam(t_pppParam::clkR, t_prn(), t_lc::dummy));
    693 
    694   // GLONASS Clock Offset
    695   // --------------------
    696   if ( _usedSystems.contains('R')  &&
    697       (_usedSystems.contains('G') || _usedSystems.contains('E') || _usedSystems.contains('C'))) {
    698     required.push_back(new t_pppParam(t_pppParam::offGlo, t_prn(), t_lc::dummy));
    699   }
    700 
    701   // Galileo Clock Offset
    702   // --------------------
    703   if (_usedSystems.contains('E') && _usedSystems.contains('G')) {
    704     required.push_back(new t_pppParam(t_pppParam::offGal, t_prn(), t_lc::dummy));
    705   }
    706 
    707   // BDS Clock Offset
    708   // ----------------
    709   if (_usedSystems.contains('C')  && (_usedSystems.contains('G') || _usedSystems.contains('E'))) {
    710     required.push_back(new t_pppParam(t_pppParam::offBds, t_prn(), t_lc::dummy));
    711   }
     675  // Receiver Clocks
     676  // ---------------
     677   if (_usedSystems.contains('G')) {
     678     required.push_back(new t_pppParam(t_pppParam::rClkG, t_prn(), t_lc::dummy));
     679   }
     680
     681   if (_usedSystems.contains('R')) {
     682     required.push_back(new t_pppParam(t_pppParam::rClkR, t_prn(), t_lc::dummy));
     683   }
     684
     685   if (_usedSystems.contains('E')) {
     686     required.push_back(new t_pppParam(t_pppParam::rClkE, t_prn(), t_lc::dummy));
     687   }
     688
     689   if (_usedSystems.contains('C')) {
     690     required.push_back(new t_pppParam(t_pppParam::rClkC, t_prn(), t_lc::dummy));
     691   }
    712692
    713693  // Troposphere
Note: See TracChangeset for help on using the changeset viewer.