Changeset 4591 in ntrip for trunk/BNC/src/rinex/reqcanalyze.cpp


Ignore:
Timestamp:
Aug 30, 2012, 12:14:37 PM (12 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/rinex/reqcanalyze.cpp

    r4590 r4591  
    375375  t_allObs& allObs = _allObsMap[prn];
    376376
    377   bncTime currTime;
    378   bncTime prevTime;
    379 
     377  // Loop over all Chunks of Data
     378  // ----------------------------
    380379  for (int chunkStart = 0; chunkStart + numEpo < allObs._oneObsVec.size();
    381380       chunkStart += chunkStep) {
    382381
    383 
    384     // Compute Mean
    385     // ------------
     382    // Chunk-Speicific Variables
     383    // -------------------------
     384    bncTime currTime;
     385    bncTime prevTime;
    386386    bncTime chunkStartTime;
    387387    bool    availL1  = false;
    388388    bool    availL2  = false;
    389     bool    missL1   = false;
    390     bool    missL2   = false;
     389    bool    gapL1    = false;
     390    bool    gapL2    = false;
    391391    bool    slipFlag = false;
    392392    double  meanMP1  = 0.0;
     
    397397    double  zenDeg   = 0.0;
    398398
     399    // Loop over all Epochs within one Chunk of Data
     400    // ---------------------------------------------
    399401    for (int ii = 0; ii < numEpo; ii++) {
    400402      int iEpo = chunkStart + ii;
     
    433435      // Check Interval
    434436      // --------------
    435       double dt = 0.0;
    436437      if (prevTime.valid()) {
    437         dt = currTime - prevTime;
    438       }
    439       if (dt != obsInterval) {
    440         missL1 = true;
    441         missL2 = true;
     438        double dt = currTime - prevTime;
     439        if (dt != obsInterval) {
     440          gapL1 = true;
     441          gapL2 = true;
     442        }
    442443      }
    443444      prevTime = currTime;
     
    449450      }
    450451      else {
    451         missL1 = true;
     452        gapL1 = true;
    452453      }
    453454      if (oneObs->_hasL2) {
     
    455456      }
    456457      else {
    457         missL2 = true;
     458        gapL2 = true;
    458459      }
    459460
     
    483484    // Availability Plot Data
    484485    // ----------------------
     486    double mjd = chunkStartTime.mjddec();
    485487    if (availL1) {
    486       _availDataMap[prn]._epoL1 << chunkStartTime.mjddec();
     488      if      (slipFlag) {
     489        _availDataMap[prn]._L1slip << mjd;
     490      }
     491      else if (gapL1) {
     492        _availDataMap[prn]._L1gap << mjd;
     493      }
     494      else {
     495        _availDataMap[prn]._L1ok << mjd;
     496      }
     497    }
     498    if (availL2) {
     499      if      (slipFlag) {
     500        _availDataMap[prn]._L2slip << mjd;
     501      }
     502      else if (gapL2) {
     503        _availDataMap[prn]._L2gap << mjd;
     504      }
     505      else {
     506        _availDataMap[prn]._L2ok << mjd;
     507      }
    487508    }
    488509
Note: See TracChangeset for help on using the changeset viewer.