Changeset 6020 in ntrip


Ignore:
Timestamp:
Aug 21, 2014, 9:19:24 AM (10 years ago)
Author:
mervart
Message:
 
Location:
trunk/BNC/src/PPP
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/PPP/pppClient.cpp

    r5914 r6020  
    206206                                  ColumnVector& xyzc, bool print) {
    207207
    208   t_lc::type tLC = (OPT->dualFreqRequired('G') ? t_lc::cIF : t_lc::c1);
     208  t_lc::type tLC = t_lc::dummy;
    209209
    210210  while (true) {
     
    213213    for (unsigned ii = 0; ii < obsVector.size(); ii++) {
    214214      const t_pppSatObs* satObs = obsVector.at(ii);
    215       if ( satObs->isValid() && satObs->prn().system() == 'G' &&
    216            (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) {
    217         ++iObs;   
    218         BB[iObs][0] = satObs->xc()[0];
    219         BB[iObs][1] = satObs->xc()[1];
    220         BB[iObs][2] = satObs->xc()[2];
    221         BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC);
     215      if (satObs->prn().system() == 'G') {
     216        if (tLC == t_lc::dummy) {
     217          tLC = satObs->isValid(t_lc::cIF) ? t_lc::cIF : t_lc::c1;
     218        }
     219        if ( satObs->isValid(tLC) && (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) {
     220          ++iObs;   
     221          BB[iObs][0] = satObs->xc()[0];
     222          BB[iObs][1] = satObs->xc()[1];
     223          BB[iObs][2] = satObs->xc()[2];
     224          BB[iObs][3] = satObs->obsValue(tLC) - satObs->cmpValueForBanc(tLC);
     225        }
    222226      }
    223227    }
     
    276280double t_pppClient::cmpOffGG(vector<t_pppSatObs*>& obsVector) {
    277281
    278   t_lc::type tLC   = (OPT->dualFreqRequired('R') ? t_lc::cIF : t_lc::c1);
     282  t_lc::type tLC   = t_lc::dummy;
    279283  double     offGG = 0.0;
    280284
     
    289293      for (unsigned ii = 0; ii < obsVector.size(); ii++) {
    290294        t_pppSatObs* satObs = obsVector.at(ii);
    291         if ( satObs->isValid() && satObs->prn().system() == 'R' &&
    292              (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle) ) {
    293           double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC);
    294           ++nObs;
    295           offGG += ll;
    296           if (fabs(ll) > fabs(maxRes)) {
    297             maxRes      = ll;
    298             maxResIndex = ii;
    299             maxResPrn   = satObs->prn();
     295        if (satObs->prn().system() == 'R') {
     296          if (tLC == t_lc::dummy) {
     297            tLC = satObs->isValid(t_lc::cIF) ? t_lc::cIF : t_lc::c1;
     298          }
     299          if (satObs->isValid(tLC) && (!satObs->modelSet() || satObs->eleSat() >= OPT->_minEle)) {
     300            double ll = satObs->obsValue(tLC) - satObs->cmpValue(tLC);
     301            ++nObs;
     302            offGG += ll;
     303            if (fabs(ll) > fabs(maxRes)) {
     304              maxRes      = ll;
     305              maxResIndex = ii;
     306              maxResPrn   = satObs->prn();
     307            }
    300308          }
    301309        }
  • trunk/BNC/src/PPP/pppSatObs.h

    r6018 r6020  
    1717  ~t_pppSatObs();
    1818  bool                isValid() const;
     19  bool                isValid(t_lc::type tLC) const;
    1920  const t_prn&        prn() const {return _prn;}
    2021  const ColumnVector& xc() const {return _xcSat;}
Note: See TracChangeset for help on using the changeset viewer.