Changeset 5750 in ntrip for trunk/BNC/src/bncmodel.cpp
- Timestamp:
- Aug 3, 2014, 11:30:39 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/bncmodel.cpp
r5656 r5750 342 342 satData->rho = (satData->xx - xRec).norm_Frobenius(); 343 343 344 double tropDelay = delay_saast( satData->eleSat) +344 double tropDelay = delay_saast(xRec, satData->eleSat) + 345 345 trp() / sin(satData->eleSat); 346 346 … … 385 385 // Tropospheric Model (Saastamoinen) 386 386 //////////////////////////////////////////////////////////////////////////// 387 double bncModel::delay_saast( double Ele) {387 double bncModel::delay_saast(const ColumnVector& xyz, double Ele) { 388 388 389 389 Tracer tracer("bncModel::delay_saast"); 390 390 391 double xyz[3];392 xyz[0] = x();393 xyz[1] = y();394 xyz[2] = z();395 391 double ell[3]; 396 xyz2ell(xyz , ell);392 xyz2ell(xyz.data(), ell); 397 393 double height = ell[2]; 398 394 … … 649 645 } 650 646 else if (par->type == bncParam::TROPO) { 651 double aprTrp = delay_saast(M_PI/2.0); 647 ColumnVector xyz(3); xyz(1) = x(); xyz(2) = y(); xyz(3) = z(); 648 double aprTrp = delay_saast(xyz, M_PI/2.0); 652 649 strB << "\n trp = " << par->prn.toAscii().data() 653 650 << setw(7) << setprecision(3) << aprTrp << " "
Note:
See TracChangeset
for help on using the changeset viewer.