Changeset 3314 in ntrip for trunk/BNC


Ignore:
Timestamp:
Jun 21, 2011, 7:27:53 PM (13 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/bncmodel.cpp

    r3313 r3314  
    6060const double   MINELE_GAL       = 10.0 * M_PI / 180.0;
    6161const double   MAXRES_CODE_GPS  = 10.0;
    62 const double   MAXRES_PHASE_GPS = 0.10;
    63 const double   MAXRES_PHASE_GLO = 0.05;
     62const double   MAXRES_PHASE_GPS = 0.03;
     63const double   MAXRES_PHASE_GLO = 0.03;
    6464const double   MAXRES_CODE_GAL  = 10.0;
    65 const double   MAXRES_PHASE_GAL = 0.10;
     65const double   MAXRES_PHASE_GAL = 0.03;
    6666
    6767// Constructor
     
    12161216  Tracer tracer("bncModel::addObs");
    12171217
     1218  const double ELEWGHT = 60.0;
     1219  double ellWgtCoef = 1.0;
     1220  double eleD = satData->eleSat * 180.0 / M_PI;
     1221  if (eleD < ELEWGHT) {
     1222    ellWgtCoef = 2.5 - 1.5 / (ELEWGHT - 10.0) * (eleD - 10.0);
     1223  }
     1224
    12181225  // Phase Observations
    12191226  // ------------------
     
    12211228    ++iObs;
    12221229    ll(iObs)      = satData->L3 - cmpValue(satData, true);
    1223     PP(iObs,iObs) = 1.0 / (_sigL3 * _sigL3);
     1230    PP(iObs,iObs) = 1.0 / (_sigL3 * _sigL3) / (ellWgtCoef * ellWgtCoef);
    12241231    for (int iPar = 1; iPar <= _params.size(); iPar++) {
    12251232      if (_params[iPar-1]->type == bncParam::AMB_L3 &&
     
    12371244    ++iObs;
    12381245    ll(iObs)      = satData->P3 - cmpValue(satData, false);
    1239     PP(iObs,iObs) = 1.0 / (_sigP3 * _sigP3);
     1246    PP(iObs,iObs) = 1.0 / (_sigP3 * _sigP3) / (ellWgtCoef * ellWgtCoef);
    12401247    for (int iPar = 1; iPar <= _params.size(); iPar++) {
    12411248      AA(iObs, iPar) = _params[iPar-1]->partial(satData, false);
Note: See TracChangeset for help on using the changeset viewer.