Index: trunk/BNC/src/bnctime.h
===================================================================
--- trunk/BNC/src/bnctime.h	(revision 4583)
+++ trunk/BNC/src/bnctime.h	(revision 4584)
@@ -21,4 +21,5 @@
   unsigned int gpsw()   const;
   double       gpssec() const;
+  double       mjddec() const {return (_mjd + _sec / 86400.0);}
   void         civil_date (unsigned int& year, unsigned int& month,
                            unsigned int& day) const;
Index: trunk/BNC/src/rinex/availplot.cpp
===================================================================
--- trunk/BNC/src/rinex/availplot.cpp	(revision 4583)
+++ trunk/BNC/src/rinex/availplot.cpp	(revision 4584)
@@ -7,10 +7,9 @@
 
 #include "availplot.h"
+#include "reqcanalyze.h"
 
 t_availPlot::t_availPlot(QWidget* parent, 
-                         QMap<QString, QVector<int> >* prnAvail) 
+                        QMap<QString, t_availData>* availDataMap) 
 : QwtPlot(parent) {
-
-  setTitle("Availability Plot");  
 
   setCanvasBackground(QColor(Qt::white));
@@ -37,10 +36,11 @@
   // ------
   int iC = 0;
-  QMapIterator<QString, QVector<int> > it(*prnAvail);
+  QMapIterator<QString, t_availData > it(*availDataMap);
   while (it.hasNext()) {
     it.next();
     ++iC;
-    const QString&      prn    = it.key();
-    const QVector<int>& epochs = it.value();
+    const QString&         prn       = it.key();
+    const t_availData&     availData = it.value();
+    const QVector<double>& epochs    = availData._epoL1;
 
     double xData[epochs.size()];
Index: trunk/BNC/src/rinex/availplot.h
===================================================================
--- trunk/BNC/src/rinex/availplot.h	(revision 4583)
+++ trunk/BNC/src/rinex/availplot.h	(revision 4584)
@@ -5,9 +5,11 @@
 #include <qwt_plot.h>
 
+class t_availData;
+
 class t_availPlot: public QwtPlot {
  Q_OBJECT
 
 public:
-  t_availPlot(QWidget* parent, QMap<QString, QVector<int> >* prnAvail);
+  t_availPlot(QWidget* parent, QMap<QString, t_availData>* availDataMap);
 
 private:
Index: trunk/BNC/src/rinex/reqcanalyze.cpp
===================================================================
--- trunk/BNC/src/rinex/reqcanalyze.cpp	(revision 4583)
+++ trunk/BNC/src/rinex/reqcanalyze.cpp	(revision 4584)
@@ -84,11 +84,6 @@
                                     const QByteArray&, double)));
 
-  connect(this, SIGNAL(dspAvailPlot(const QString&, 
-                                    const QByteArray&,
-                                    QMap<QString, QVector<int> >*)),
-          this, SLOT(slotDspAvailPlot(const QString&, 
-                                      const QByteArray&,
-                                      QMap<QString, QVector<int> >*)));
-
+  connect(this, SIGNAL(dspAvailPlot(const QString&, const QByteArray&)),
+          this, SLOT(slotDspAvailPlot(const QString&, const QByteArray&)));
 }
 
@@ -214,5 +209,6 @@
   }
 
-  _satStat.clear();
+  _allObsMap.clear();
+  _availDataMap.clear();
 
   // A priori Coordinates
@@ -239,7 +235,5 @@
                                      .arg(obs.satNum, 2, 10, QChar('0'));
   
-        t_satStat& satStat = _satStat[prn];
-  
-        satStat.addObs(obs);
+        _allObsMap[prn].addObs(obs);
       }
   
@@ -262,15 +256,11 @@
   QVector<t_polarPoint*>*       dataSNR1 = new QVector<t_polarPoint*>;
   QVector<t_polarPoint*>*       dataSNR2 = new QVector<t_polarPoint*>;
-  QMap<QString, QVector<int> >* availL1  = new QMap<QString, QVector<int> >;
-
-
-  QMapIterator<QString, t_satStat> it(_satStat);
+
+  QMutableMapIterator<QString, t_allObs> it(_allObsMap);
   while (it.hasNext()) {
     it.next();
-    QString          prn     = it.key();
-    const t_satStat& satStat = it.value();
-    QVector<int>&    dataL1  = (*availL1)[prn];
-    preparePlotData(prn, satStat, xyz, obsFile->interval(), 
-                    dataMP1, dataMP2, dataSNR1, dataSNR2, dataL1);
+    QString    prn     = it.key();
+    preparePlotData(prn, xyz, obsFile->interval(), 
+                    dataMP1, dataMP2, dataSNR1, dataSNR2);
   }
 
@@ -281,5 +271,5 @@
                   "", 9.0);
 
-  emit dspAvailPlot(obsFile->fileName(), "Availability L1", availL1);
+  emit dspAvailPlot(obsFile->fileName(), "Availability L1");
 
   if (_log) {
@@ -290,7 +280,7 @@
 //  
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::t_satStat::addObs(const t_obs& obs) { 
-
-  t_anaObs* newObs = new t_anaObs(obs.GPSWeek, obs.GPSWeeks);
+void t_reqcAnalyze::t_allObs::addObs(const t_obs& obs) { 
+
+  t_oneObs* newObs = new t_oneObs(obs.GPSWeek, obs.GPSWeeks);
   bool      okFlag = false;
 
@@ -364,5 +354,5 @@
   // ------------------------
   if (okFlag) {
-    anaObs << newObs;
+    _oneObsVec << newObs;
   }
   else {
@@ -373,18 +363,17 @@
 //  
 ////////////////////////////////////////////////////////////////////////////
-void t_reqcAnalyze::preparePlotData(const QString& prn, 
-                                    const t_satStat& satStat,
-                                    const ColumnVector& xyz,
+void t_reqcAnalyze::preparePlotData(const QString& prn, const ColumnVector& xyz,
                                     double obsInterval,
                                     QVector<t_polarPoint*>* dataMP1, 
                                     QVector<t_polarPoint*>* dataMP2,
                                     QVector<t_polarPoint*>* dataSNR1, 
-                                    QVector<t_polarPoint*>* dataSNR2,
-                                    QVector<int>&           dataL1) {
+                                    QVector<t_polarPoint*>* dataSNR2) {
 
   const int chunkStep = int( 30.0 / obsInterval); // chunk step (30 sec)  
   const int numEpo    = int(600.0 / obsInterval); // # epochs in one chunk (10 min)
 
-  for (int chunkStart = 0; chunkStart + numEpo < satStat.anaObs.size();
+  t_allObs& allObs = _allObsMap[prn];
+
+  for (int chunkStart = 0; chunkStart + numEpo < allObs._oneObsVec.size();
        chunkStart += chunkStep) {
 
@@ -401,18 +390,19 @@
     for (int ii = 0; ii < numEpo; ii++) {
       int iEpo = chunkStart + ii;
-      const t_anaObs* anaObs = satStat.anaObs[iEpo];
+      const t_oneObs* oneObs = allObs._oneObsVec[iEpo];
 
       if (ii == 0) {
-        firstEpoch.set(anaObs->_GPSWeek, anaObs->_GPSWeeks);
-      }
-
-      mean1 += anaObs->_MP1;
-      mean2 += anaObs->_MP2;
-
-      if ( anaObs->_SNR1 > 0 && (SNR1 == 0 || SNR1 > anaObs->_SNR1) ) {
-        SNR1 = anaObs->_SNR1;
-      }
-      if ( anaObs->_SNR2 > 0 && (SNR2 == 0 || SNR2 > anaObs->_SNR2) ) {
-        SNR2 = anaObs->_SNR2;
+        bncTime epoTime(oneObs->_GPSWeek, oneObs->_GPSWeeks);
+        _availDataMap[prn]._epoL1 << epoTime.mjddec();
+      }
+
+      mean1 += oneObs->_MP1;
+      mean2 += oneObs->_MP2;
+
+      if ( oneObs->_SNR1 > 0 && (SNR1 == 0 || SNR1 > oneObs->_SNR1) ) {
+        SNR1 = oneObs->_SNR1;
+      }
+      if ( oneObs->_SNR2 > 0 && (SNR2 == 0 || SNR2 > oneObs->_SNR2) ) {
+        SNR2 = oneObs->_SNR2;
       }
   
@@ -420,6 +410,6 @@
       // ----------
       if (ii > 0) {
-        double diff1 = anaObs->_MP1 - satStat.anaObs[iEpo-1]->_MP1;
-        double diff2 = anaObs->_MP2 - satStat.anaObs[iEpo-1]->_MP2;
+        double diff1 = oneObs->_MP1 - allObs._oneObsVec[iEpo-1]->_MP1;
+        double diff2 = oneObs->_MP2 - allObs._oneObsVec[iEpo-1]->_MP2;
         if (fabs(diff1) > SLIPTRESH || fabs(diff2) > SLIPTRESH) {
           slipFlag = true;
@@ -428,6 +418,4 @@
       }
     }
-
-    dataL1 << int(firstEpoch.gpssec());
 
     if (slipFlag) {
@@ -444,7 +432,7 @@
     for (int ii = 0; ii < numEpo; ii++) {
       int iEpo = chunkStart + ii;
-      const t_anaObs* anaObs = satStat.anaObs[iEpo];
-      double diff1 = anaObs->_MP1 - mean1;
-      double diff2 = anaObs->_MP2 - mean2;
+      const t_oneObs* oneObs = allObs._oneObsVec[iEpo];
+      double diff1 = oneObs->_MP1 - mean1;
+      double diff2 = oneObs->_MP2 - mean2;
       stddev1 += diff1 * diff1;
       stddev2 += diff2 * diff2;
@@ -453,5 +441,5 @@
     double MP2 = sqrt(stddev2 / (numEpo-1));
 
-    const t_anaObs* anaObs0 = satStat.anaObs[chunkStart];
+    const t_oneObs* oneObs0 = allObs._oneObsVec[chunkStart];
 
     // Compute the Azimuth and Zenith Distance
@@ -470,5 +458,5 @@
       if (eph) {
         double xSat, ySat, zSat, clkSat;
-        eph->position(anaObs0->_GPSWeek, anaObs0->_GPSWeeks, 
+        eph->position(oneObs0->_GPSWeek, oneObs0->_GPSWeeks, 
                       xSat, ySat, zSat, clkSat);
       
@@ -509,10 +497,10 @@
 ////////////////////////////////////////////////////////////////////////////
 void t_reqcAnalyze::slotDspAvailPlot(const QString& fileName, 
-                                     const QByteArray& title,
-                                     QMap<QString, QVector<int> >* prnAvail){
+                                     const QByteArray& title) {
 
   if (dynamic_cast<bncApp*>(qApp)->GUIenabled()) {
 
-    t_availPlot* plot = new t_availPlot(0, prnAvail);
+    t_availPlot* plot = new t_availPlot(0, &_availDataMap);
+    plot->setTitle(title);
 
     QVector<QWidget*> plots;
Index: trunk/BNC/src/rinex/reqcanalyze.h
===================================================================
--- trunk/BNC/src/rinex/reqcanalyze.h	(revision 4583)
+++ trunk/BNC/src/rinex/reqcanalyze.h	(revision 4584)
@@ -34,4 +34,9 @@
 class t_polarPoint;
 
+class t_availData {
+ public:
+  QVector<double> _epoL1;
+};
+
 class t_reqcAnalyze : public QThread {
 Q_OBJECT
@@ -52,6 +57,5 @@
                   const QByteArray&, double);
 
-  void dspAvailPlot(const QString&, const QByteArray&,
-                    QMap<QString, QVector<int> >*);
+  void dspAvailPlot(const QString&, const QByteArray&);
    
  private slots:
@@ -63,6 +67,5 @@
                       const QByteArray& scaleTitle, double maxValue);
 
-  void slotDspAvailPlot(const QString& fileName, const QByteArray& title,
-                        QMap<QString, QVector<int> >* prnAvail);
+  void slotDspAvailPlot(const QString& fileName, const QByteArray& title);
 
  public:
@@ -70,7 +73,7 @@
  
  private:
-  class t_anaObs {
+  class t_oneObs {
    public:
-    t_anaObs(int GPSWeek, double GPSWeeks) {
+    t_oneObs(int GPSWeek, double GPSWeeks) {
       _GPSWeek  = GPSWeek;
       _GPSWeeks = GPSWeeks;
@@ -92,34 +95,34 @@
   };
 
-  class t_satStat {
+  class t_allObs {
    public:
-    t_satStat() {}
-    ~t_satStat() {
-      for (int ii = 0; ii < anaObs.size(); ii++) {
-        delete anaObs[ii];
+    t_allObs() {}
+    ~t_allObs() {
+      for (int ii = 0; ii < _oneObsVec.size(); ii++) {
+        delete _oneObsVec[ii];
       }
     }
     void addObs(const t_obs& obs);
-    QVector<t_anaObs*> anaObs;
+    QVector<t_oneObs*> _oneObsVec;
   };
 
   void analyzeFile(t_rnxObsFile* obsFile);
-  void preparePlotData(const QString& prn, const t_satStat& satStat,
-                       const ColumnVector& xyz, double obsInterval,
+  void preparePlotData(const QString& prn, const ColumnVector& xyz, 
+                       double obsInterval,
                        QVector<t_polarPoint*>* dataMP1, 
                        QVector<t_polarPoint*>* dataMP2,
                        QVector<t_polarPoint*>* dataSNR1, 
-                       QVector<t_polarPoint*>* dataSNR2,
-                       QVector<int>&           dataL1);
+                       QVector<t_polarPoint*>* dataSNR2);
 
-  QString                  _logFileName;
-  QFile*                   _logFile;
-  QTextStream*             _log;
-  QStringList              _obsFileNames;
-  QVector<t_rnxObsFile*>   _rnxObsFiles;
-  QStringList              _navFileNames;
-  QVector<t_eph*>          _ephs;
-  QMap<QString, t_satStat> _satStat;
-  t_rnxObsFile::t_rnxEpo*  _currEpo;
+  QString                    _logFileName;
+  QFile*                     _logFile;
+  QTextStream*               _log;
+  QStringList                _obsFileNames;
+  QVector<t_rnxObsFile*>     _rnxObsFiles;
+  QStringList                _navFileNames;
+  QVector<t_eph*>            _ephs;
+  t_rnxObsFile::t_rnxEpo*    _currEpo;
+  QMap<QString, t_allObs>    _allObsMap;
+  QMap<QString, t_availData> _availDataMap;
 };