Ignore:
Timestamp:
Feb 12, 2013, 11:29:38 AM (12 years ago)
Author:
weber
Message:

Clock correction due to heavy scale in GDA94 transformation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/upload/bncrtnetuploadcaster.cpp

    r4808 r4908  
    574574                                    double& dc) {
    575575
    576   // Clock Correction
    577   // ----------------
    578   dc = 0.0; // TODO
    579 
    580576  // Current epoch minus 2000.0 in years
    581577  // ------------------------------------
     
    595591
    596592  double sc = 1.0 + _sc * 1e-9 + dt * _scr * 1e-9;
     593
     594  // Clock correction due to heavy scale in GDA94 transformation
     595  // -----------------------------------------------------------
     596  if (_crdTrafo != "GDA94") {
     597    dc = 0.0;
     598  }
     599  else {
     600
     601    // Spefify mean of Australia (here: GDA94 coordinates of site ALIC)
     602    // ----------------------------------------------------------------
     603    ColumnVector meanSta(3);
     604    meanSta(1) = -4052051.767;
     605    meanSta(2) =  4212836.197;
     606    meanSta(3) = -2545106.022;
     607
     608    // Correction proportional to topocentric distance to satellites
     609    // -------------------------------------------------------------
     610    double rho = (xyz - meanSta).norm_Frobenius();
     611    dc = rho * (sc - 1.0) / t_CST::c;
     612  }
    597613
    598614  Matrix rMat(3,3);
Note: See TracChangeset for help on using the changeset viewer.