Changeset 6254 in ntrip for trunk


Ignore:
Timestamp:
Oct 23, 2014, 10:16:27 AM (10 years ago)
Author:
stuerze
Message:

modeled zenith tropo delay added in ppp logfile output

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/PPP/pppParlist.cpp

    r6107 r6254  
    1111 * Created:    29-Jul-2014
    1212 *
    13  * Changes:   
     13 * Changes:
    1414 *
    1515 * -----------------------------------------------------------------------*/
     
    103103}
    104104
    105 // 
    106 ////////////////////////////////////////////////////////////////////////////
    107 double t_pppParam::partial(const bncTime& /* epoTime */, const t_pppSatObs* obs, 
     105//
     106////////////////////////////////////////////////////////////////////////////
     107double t_pppParam::partial(const bncTime& /* epoTime */, const t_pppSatObs* obs,
    108108                        const t_lc::type& tLC) const {
    109109
     
    141141        obs->lcCoeff(tLC, codeCoeff, phaseCoeff);
    142142        if      (_tLC == t_lc::l1) {
    143           return obs->lambda(t_lc::l1) * phaseCoeff[t_lc::toFreq(obs->prn().system(),t_lc::l1)]; 
     143          return obs->lambda(t_lc::l1) * phaseCoeff[t_lc::toFreq(obs->prn().system(),t_lc::l1)];
    144144        }
    145145        else if (_tLC == t_lc::l2) {
    146           return obs->lambda(t_lc::l2) * phaseCoeff[t_lc::toFreq(obs->prn().system(),t_lc::l2)]; 
     146          return obs->lambda(t_lc::l2) * phaseCoeff[t_lc::toFreq(obs->prn().system(),t_lc::l2)];
    147147        }
    148148      }
     
    150150    return 0.0;
    151151  case trp:
    152     return 1.0 / sin(obs->eleSat()); 
     152    return 1.0 / sin(obs->eleSat());
    153153  }
    154154
     
    156156}
    157157
    158 // 
     158//
    159159////////////////////////////////////////////////////////////////////////////
    160160string t_pppParam::toString() const {
     
    163163  case crdX:
    164164    ss << "CRD_X";
    165     break; 
     165    break;
    166166  case crdY:
    167167    ss << "CRD_Y";
     
    199199}
    200200
    201 // 
     201//
    202202////////////////////////////////////////////////////////////////////////////
    203203t_irc t_pppParlist::set(const bncTime& epoTime, const std::vector<t_pppSatObs*>& obsVector) {
     
    211211    bool remove = false;
    212212
    213     if      (par->epoSpec()) { 
     213    if      (par->epoSpec()) {
    214214      remove = true;
    215215    }
     
    258258      }
    259259    }
    260   } 
     260  }
    261261
    262262  // Required Set of Parameters
     
    336336}
    337337
    338 // 
     338//
    339339////////////////////////////////////////////////////////////////////////////
    340340void t_pppParlist::printResult(const bncTime& epoTime, const SymmetricMatrix& QQ,
     
    342342
    343343  string epoTimeStr = string(epoTime);
     344  const t_pppStation* sta = PPP_CLIENT->staRover();
    344345
    345346  LOG << endl;
     
    361362    else {
    362363      int ind = par->indexNew();
     364      double apr = (par->type() == t_pppParam::trp) ?
     365        t_tropo::delay_saast(sta->xyzApr(), M_PI/2.0) :  par->x0();
    363366      LOG << epoTimeStr << ' ' << par->toString() << ' '
    364           << setw(10) << setprecision(4) << par->x0() << ' '
     367          << setw(10) << setprecision(4) << apr << ' '
    365368          << showpos << setw(10) << setprecision(4) << xx[ind] << noshowpos << " +- "
    366369          << setw(8)  << setprecision(4) << sqrt(QQ[ind][ind]);
     
    372375    }
    373376  }
    374  
     377
    375378  if (parX && parY && parZ) {
    376     const t_pppStation* sta = PPP_CLIENT->staRover();
    377 
    378     ColumnVector xyz(3);
     379
     380        ColumnVector xyz(3);
    379381    xyz[0] = xx[parX->indexNew()];
    380382    xyz[1] = xx[parY->indexNew()];
     
    392394        << " X = " << setprecision(4) << sta->xyzApr()[0] + xyz[0] << " +- "
    393395        << setprecision(4) << sqrt(QQxyz[0][0])
    394                                    
     396
    395397        << " Y = " << setprecision(4) << sta->xyzApr()[1] + xyz[1] << " +- "
    396398        << setprecision(4) << sqrt(QQxyz[1][1])
     
    407409        << " dU = " << setprecision(4) << neu[2] << " +- "
    408410        << setprecision(4) << sqrt(QQneu[2][2])
    409  
     411
    410412        << endl;
    411413  }
Note: See TracChangeset for help on using the changeset viewer.