Index: trunk/BNC/src/RTCM/rtcm_utils.cpp
===================================================================
--- trunk/BNC/src/RTCM/rtcm_utils.cpp	(revision 8234)
+++ trunk/BNC/src/RTCM/rtcm_utils.cpp	(revision 8483)
@@ -10,12 +10,12 @@
 
 void resolveEpoch (double secsHour,
-                   int  refWeek,   double  refSecs,  
+                   int  refWeek,   double  refSecs,
                    int& epochWeek, double& epochSecs) {
 
-  const double secsPerWeek = 604800.0;                            
+  const double secsPerWeek = 604800.0;
 
   epochWeek = refWeek;
   epochSecs = secsHour + 3600.0*(floor((refSecs-secsHour)/3600.0+0.5));
-  
+
   if (epochSecs<0          ) { epochWeek--; epochSecs+=secsPerWeek; };
   if (epochSecs>secsPerWeek) { epochWeek++; epochSecs-=secsPerWeek; };
@@ -29,11 +29,11 @@
            double& xSat, double& ySat, double& zSat, double& clkSat) {
 
-  const double omega_earth = 7292115.1467e-11; 
-  const double secsPerWeek = 604800.0;                            
+  const double omega_earth = 7292115.1467e-11;
+  const double secsPerWeek = 604800.0;
 
   // Initial values
   // --------------
   rho = 0.0;
-  ColumnVector xc(4);
+  ColumnVector xc(7);
   ColumnVector vv(3);
   eph->getCrd(bncTime(GPSWeek, GPSWeeks), xc, vv, false);
@@ -45,5 +45,5 @@
   ////cout << "----- cmpRho -----\n";
   ////eph->print(cout);
-  ////cout << "  pos " << setw(4)  << GPSWeek 
+  ////cout << "  pos " << setw(4)  << GPSWeek
   ////     << " "      << setw(14) << setprecision(6) << GPSWeeks
   ////     << " "      << setw(13) << setprecision(3) << xSat
@@ -57,10 +57,10 @@
   do {
     rhoLast = rho;
-    
+
     // Correction station position due to Earth Rotation
     // -------------------------------------------------
     double dPhi = omega_earth * rho / c_light;
-    double xRec = stax * cos(dPhi) - stay * sin(dPhi); 
-    double yRec = stay * cos(dPhi) + stax * sin(dPhi); 
+    double xRec = stax * cos(dPhi) - stay * sin(dPhi);
+    double yRec = stay * cos(dPhi) + stax * sin(dPhi);
     double zRec = staz;
 
@@ -81,5 +81,5 @@
       GPSWeek_tot  += 1;
     }
-      
+
     eph->getCrd(bncTime(GPSWeek_tot, GPSWeeks_tot), xc, vv, false);
     xSat   = xc(1);
@@ -94,5 +94,5 @@
     rho = sqrt(dx*dx + dy*dy + dz*dz);
 
-    ////cout << "  scrd "   << setw(4)  << GPSWeek_tot 
+    ////cout << "  scrd "   << setw(4)  << GPSWeek_tot
     ////	 << " "         << setw(15) << setprecision(8) << GPSWeeks_tot
     ////	 << " "         << setw(13) << setprecision(3) << xSat
@@ -108,9 +108,9 @@
     ////	 << " rho "     << setw(13) << setprecision(3) << rho
     ////	 << endl;
-    
+
 
     ////cout.setf(ios::fixed);
     ////
-    ////cout << "niter " << setw(3) << ++niter 
+    ////cout << "niter " << setw(3) << ++niter
     ////         << " " << setw(14) << setprecision(3) << rhoLast
     ////         << " " << setw(14) << setprecision(3) << rho