Index: trunk/BNC/src/PPP/pppSatObs.cpp
===================================================================
--- trunk/BNC/src/PPP/pppSatObs.cpp	(revision 9447)
+++ trunk/BNC/src/PPP/pppSatObs.cpp	(revision 9473)
@@ -17,4 +17,5 @@
 
 #include <iostream>
+#include <iomanip>
 #include <cmath>
 #include <newmatio.h>
@@ -304,9 +305,9 @@
 
   map<t_frequency::type, double>::const_iterator it;
-  for (it = codeCoeff.begin(); it != codeCoeff.end(); it++)   {//cout << "codeCoeff : " << t_frequency::toString(it->first).c_str() << ": " << it->second << endl;
+  for (it = codeCoeff.begin(); it != codeCoeff.end(); it++)   {
     retVal += it->second * it->second * OPT->_sigmaC1 * OPT->_sigmaC1;
   }
 
-  for (it = phaseCoeff.begin(); it != phaseCoeff.end(); it++) {//cout << "phaseCoeff: " << t_frequency::toString(it->first).c_str()  << ": " << it->second << endl;
+  for (it = phaseCoeff.begin(); it != phaseCoeff.end(); it++) {
     retVal += it->second * it->second * OPT->_sigmaL1 * OPT->_sigmaL1;
   }
@@ -314,10 +315,16 @@
   retVal = sqrt(retVal);
 
-  // De-Weight GLONASS and BDS
-  // --------------------------
-  if ((_prn.system() == 'R' && (t_lc::includesCode(tLC) ||  t_lc::includesPhase(tLC))) ||
-      (_prn.system() == 'C' && (t_lc::includesCode(tLC) ||  t_lc::includesPhase(tLC)))) {
+  // De-Weight GLONASS
+  // -----------------
+  if (_prn.system() == 'R' && (t_lc::includesCode(tLC) ||  t_lc::includesPhase(tLC))) {
     retVal *= 5.0;
   }
+
+  // De-Weight BDS
+  // -------------
+  if (_prn.system() == 'C' && (t_lc::includesCode(tLC) ||  t_lc::includesPhase(tLC))) {
+    retVal *= 2.0;
+  }
+
 
   // Elevation-Dependent Weighting
@@ -345,8 +352,8 @@
 
   map<t_frequency::type, double>::const_iterator it;
-  for (it = codeCoeff.begin(); it != codeCoeff.end(); it++)   {//cout << "codeCoeff: " << it->first << ": " << it->second << endl;
+  for (it = codeCoeff.begin(); it != codeCoeff.end(); it++)   {
     retVal += it->second * it->second * OPT->_maxResC1 * OPT->_maxResC1;
   }
-  for (it = phaseCoeff.begin(); it != phaseCoeff.end(); it++) {//cout << "phaseCoeff: " << it->first << ": " << it->second << endl;
+  for (it = phaseCoeff.begin(); it != phaseCoeff.end(); it++) {
     retVal += it->second * it->second * OPT->_maxResL1 * OPT->_maxResL1;
   }