Index: trunk/BNC/bncmodel.cpp
===================================================================
--- trunk/BNC/bncmodel.cpp	(revision 2669)
+++ trunk/BNC/bncmodel.cpp	(revision 2670)
@@ -724,4 +724,10 @@
                             epoData->satDataGlo) != 0);
 
+  // Remember the Epoch-specific Results for the computation of means
+  // ----------------------------------------------------------------
+  pppPos* newPos = new pppPos;
+  newPos->time   = epoData->tt;
+  _posAverage.push_back(newPos); 
+
   // Set Solution Vector
   // -------------------
@@ -745,9 +751,11 @@
     }
     else if (par->type == bncParam::TROPO) {
+      double aprTrp = delay_saast(M_PI/2.0);
       strB << "\n    trp     = " << par->prn.toAscii().data()
-           << setw(7) << setprecision(3) << delay_saast(M_PI/2.0) << " "
+           << setw(7) << setprecision(3) << aprTrp << " "
            << setw(6) << setprecision(3) << showpos << par->xx << noshowpos
            << " +- " << setw(6) << setprecision(3) 
            << sqrt(_QQ(par->index,par->index));
+      newPos->xnt[6] = aprTrp + par->xx;
     }
   }
@@ -755,5 +763,4 @@
   _log += strB.str().c_str();
   emit newMessage(_log, false);
-
 
   // Final Message (both log file and screen)
@@ -780,20 +787,17 @@
     xyzRef[2] = settings.value("pppRefCrdZ").toDouble();
 
-    pppPos* newPos = new pppPos;
-    newPos->time   = epoData->tt;
-    newPos->xn[0] = x() - xyzRef[0];
-    newPos->xn[1] = y() - xyzRef[1];
-    newPos->xn[2] = z() - xyzRef[2];
+    newPos->xnt[0] = x() - xyzRef[0];
+    newPos->xnt[1] = y() - xyzRef[1];
+    newPos->xnt[2] = z() - xyzRef[2];
 
     double ellRef[3];
     xyz2ell(xyzRef, ellRef);
-    xyz2neu(ellRef, newPos->xn, &newPos->xn[3]);
+    xyz2neu(ellRef, newPos->xnt, &newPos->xnt[3]);
 
     strC << "  NEU "
-         << setw(8) << setprecision(3) << newPos->xn[3] << " "
-         << setw(8) << setprecision(3) << newPos->xn[4] << " "
-         << setw(8) << setprecision(3) << newPos->xn[5];
-
-    _posAverage.push_back(newPos); // remember for the computation of mean
+         << setw(8) << setprecision(3) << newPos->xnt[3] << " "
+         << setw(8) << setprecision(3) << newPos->xnt[4] << " "
+         << setw(8) << setprecision(3) << newPos->xnt[5];
+
   }
 
@@ -814,5 +818,5 @@
     // Compute the Mean
     // ----------------
-    ColumnVector mean(6); mean = 0.0;
+    ColumnVector mean(7); mean = 0.0;
 
     QMutableVectorIterator<pppPos*> it(_posAverage);
@@ -824,6 +828,6 @@
       }
       else {
-        for (int ii = 0; ii < 6; ++ii) {
-          mean[ii] += pp->xn[ii];
+        for (int ii = 0; ii < 7; ++ii) {
+          mean[ii] += pp->xnt[ii];
         }
       }
@@ -838,13 +842,13 @@
       // Compute the Deviation
       // ---------------------
-      ColumnVector std(6); std = 0.0;
+      ColumnVector std(7); std = 0.0;
       QVectorIterator<pppPos*> it2(_posAverage);
       while (it2.hasNext()) {
         pppPos* pp = it2.next();
-        for (int ii = 0; ii < 6; ++ii) {
-          std[ii] += (pp->xn[ii] - mean[ii]) * (pp->xn[ii] - mean[ii]);
-        }
-      }
-      for (int ii = 0; ii < 6; ++ii) {
+        for (int ii = 0; ii < 7; ++ii) {
+          std[ii] += (pp->xnt[ii] - mean[ii]) * (pp->xnt[ii] - mean[ii]);
+        }
+      }
+      for (int ii = 0; ii < 7; ++ii) {
         std[ii] = sqrt(std[ii] / nn);
       }
@@ -869,7 +873,15 @@
            << setw(6)  << setprecision(3) << std[4]   << " "
            << setw(14) << setprecision(3) << mean[5]  << " +- "
-           << setw(6)  << setprecision(3) << std[5]   << endl;
+           << setw(6)  << setprecision(3) << std[5];
 
       emit newMessage(QByteArray(strE.str().c_str()), true);
+
+      ostringstream strF; strF.setf(ios::fixed);
+      strF << _staID.data() << "  AVE-TRP " 
+           << epoData->tt.timestr(1) << " "
+           << setw(13) << setprecision(3) << mean[6]  << " +- "
+           << setw(6)  << setprecision(3) << std[6]   << endl;
+
+      emit newMessage(QByteArray(strF.str().c_str()), true);
     }
   }
Index: trunk/BNC/bncmodel.h
===================================================================
--- trunk/BNC/bncmodel.h	(revision 2669)
+++ trunk/BNC/bncmodel.h	(revision 2670)
@@ -99,6 +99,11 @@
   class pppPos {
    public:
+    pppPos() {
+      for (int ii = 0; ii < 7; ++ii) {
+        xnt[ii] = 0.0;
+      }
+    }
     bncTime time;
-    double  xn[6];
+    double  xnt[7];
   };