Changeset 2429 in ntrip


Ignore:
Timestamp:
Apr 15, 2010, 1:42:35 PM (8 years ago)
Author:
mervart
Message:

* empty log message *

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/RTCM3/ephemeris.cpp

    r2285 r2429  
    122122    tc += (GPSweek - _GPSweek) * secPerWeek;
    123123  }
    124   xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc
    125           - 4.442807633e-10 * _e * sqrt(a0) *sin(E);
     124  xc[3] = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
    126125
    127126  // Velocity
     
    148147
    149148  vv[2]  = sini    *doty  + yp*cosi      *doti;
     149
     150  // Relativistic Correction
     151  // -----------------------
     152  //  xc(4) -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
     153  xc[3] -= 2.0 * (xc[0]*vv[0] + xc[1]*vv[1] + xc[2]*vv[2]) / t_CST::c / t_CST::c;
    150154}
    151155
  • trunk/BNC/bncpppclient.cpp

    r2428 r2429  
    374374        in >> bType >> bValue;
    375375        if      (bType ==  CODETYPEGPS_L1_Z) {
    376           bb->p1 = bValue * 0.299792458;
     376          bb->p1 = bValue;
    377377        }
    378378        else if (bType ==  CODETYPEGPS_L1_CA) {
    379           bb->c1 = bValue * 0.299792458;
     379          bb->c1 = bValue;
    380380        }
    381381        else if (bType == CODETYPEGPS_L2_Z) {
    382           bb->p2 = bValue * 0.299792458;
     382          bb->p2 = bValue;
    383383        }
    384384      }
  • trunk/BNS/bns.cpp

    r2427 r2429  
    426426              xx(5) *= 1e-6;         // rel. corr.
    427427                                     // xx(6), xx(7), xx(8) ... PhaseCent - CoM
    428               xx(9)  /= 0.299792458; // xx(9)  ... P1-C1 DCB in ns
    429               xx(10) /= 0.299792458; // xx(10) ... P1-P2 DCB in ns
     428                                     // xx(9)  ... P1-C1 DCB in meters
     429                                     // xx(10) ... P1-P2 DCB in meters
    430430                                     // xx(11) ... dT
    431431              xx(12) *= 1e3;         // x-crd at time + dT
  • trunk/BNS/bnseph.cpp

    r2285 r2429  
    454454    tc += (GPSweek - _GPSweek) * secPerWeek;
    455455  }
    456   xc(4) = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc
    457           - 4.442807633e-10 * _e * sqrt(a0) *sin(E);
     456  xc(4) = _clock_bias + _clock_drift*tc + _clock_driftrate*tc*tc;
    458457
    459458  // Velocity
     
    480479
    481480  vv(3)  = sini    *doty  + yp*cosi      *doti;
     481
     482  // Relativistic Correction
     483  // -----------------------
     484  //  xc(4) -= 4.442807633e-10 * _e * sqrt(a0) *sin(E);
     485  const static double c_vel = 299792458.0;
     486  xc(4) -= 2.0 * DotProduct(xc.Rows(1,3),vv) / c_vel / c_vel;
    482487}
    483488
Note: See TracChangeset for help on using the changeset viewer.