Index: /trunk/BNC/bncmodel.cpp
===================================================================
--- /trunk/BNC/bncmodel.cpp	(revision 2065)
+++ /trunk/BNC/bncmodel.cpp	(revision 2066)
@@ -185,6 +185,4 @@
   double tropDelay = delay_saast(satData->eleSat);
 
-  cout << "tropDelay " << tropDelay << endl;
-
   return satData->rho + _xcBanc(4) - satData->clk + tropDelay;
 }
@@ -200,4 +198,7 @@
   double hh =  50.0 * exp(-6.396e-4 * height);
   double ee =  hh / 100.0 * exp(-37.2465 + 0.213166*TT - 0.000256908*TT*TT);
+
+  cout << "Ele: " << 180.0 * Ele / M_PI << " " << height << " " 
+       << pp << " " << TT << " " << hh << " " << ee << endl;
 
   double h_km = height / 1000.0;
@@ -230,4 +231,6 @@
   unsigned nObs = epoData->size();
 
+  // Create First-Design Matrix
+  // --------------------------
   _AA.ReSize(nObs, nPar);  // variance-covariance matrix
   _ll.ReSize(nObs);        // tems observed-computed
@@ -251,4 +254,6 @@
   }
 
+  // Compute Least-Squares Solution
+  // ------------------------------
   _QQ.ReSize(nPar);
   _QQ << _AA.t() * _AA;
@@ -256,6 +261,14 @@
   _dx = _QQ * _AA.t() * _ll;
 
+  // Compute Residuals
+  // -----------------
+  ColumnVector vv = _AA * _dx - _ll;
+
+  cout << setprecision(3) << vv.t(); 
+  cout.flush();
+
+  // Set Solution Vector
+  // -------------------
   _xx.ReSize(nPar);
-
   unsigned iPar = 0;
   QListIterator<bncParam*> itPar(_params);