Ignore:
Timestamp:
Mar 20, 2013, 11:30:12 AM (12 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/combination/bnccomb.cpp

    r4978 r4992  
    678678    t_corr* corr = it.value();
    679679
    680     double dT = 60.0;
    681 
    682680    ColumnVector xc(4);
    683681    ColumnVector vv(3);
     
    686684    bncPPPclient::applyCorr(_resTime, corr, xc, vv);
    687685   
    688     // Relativistic Correction
    689     // -----------------------
    690     double relCorr = - 2.0 * DotProduct(xc.Rows(1,3),vv) / t_CST::c / t_CST::c;
    691     xc(4) -= relCorr;
    692    
    693     // Code Biases
    694     // -----------
    695     double dcbP1C1 = 0.0;
    696     double dcbP1P2 = 0.0;
    697    
    698686    // Correction Phase Center --> CoM
    699687    // -------------------------------
     
    701689    if (_antex) {
    702690      double Mjd = _resTime.mjd() + _resTime.daysec()/86400.0;
    703       if (_antex->satCoMcorrection(corr->prn, Mjd, xc.Rows(1,3), dx) == success) {
    704         xc(1) -= dx(1);
    705         xc(2) -= dx(2);
    706         xc(3) -= dx(3);
    707       }
    708       else {
     691      if (_antex->satCoMcorrection(corr->prn, Mjd, xc.Rows(1,3), dx) != success) {
     692        dx = 0;
    709693        cout << "antenna not found" << endl;
    710694      }
    711695    }
    712696   
    713     out << 'P' << corr->prn.toAscii().data()
    714         << setw(14) << setprecision(6) << xc(1) / 1000.0
    715         << setw(14) << setprecision(6) << xc(2) / 1000.0
    716         << setw(14) << setprecision(6) << xc(3) / 1000.0
    717         << setw(14) << setprecision(6) << xc(4) * 1e6
    718         << setw(14) << setprecision(6) << relCorr * 1e6
    719         << setw(8)  << setprecision(3) << dx(1)
    720         << setw(8)  << setprecision(3) << dx(2)
    721         << setw(8)  << setprecision(3) << dx(3)
    722         << setw(8)  << setprecision(3) << dcbP1C1
    723         << setw(8)  << setprecision(3) << dcbP1P2
    724         << setw(6)  << setprecision(1) << dT
    725         << setw(14) << setprecision(6) << (xc(1) + vv(1) * dT) / 1000.0
    726         << setw(14) << setprecision(6) << (xc(2) + vv(2) * dT) / 1000.0
    727         << setw(14) << setprecision(6) << (xc(3) + vv(3) * dT) / 1000.0 << endl;
    728 
     697    out << corr->prn.toAscii().data()
     698        << " APC 3 "
     699        << setw(15) << setprecision(4) << xc(1)
     700        << setw(15) << setprecision(4) << xc(2)
     701        << setw(15) << setprecision(4) << xc(3)
     702        << " Clk 1 "
     703        << setw(15) << setprecision(4) << xc(4) * t_CST::c
     704        << " Vel 3 "
     705        << setw(15) << setprecision(4) << vv(1)
     706        << setw(15) << setprecision(4) << vv(2)
     707        << setw(15) << setprecision(4) << vv(3)
     708        << " CoM 3 "
     709        << setw(15) << setprecision(4) << xc(1) - dx(1)
     710        << setw(15) << setprecision(4) << xc(2) - dx(2)
     711        << setw(15) << setprecision(4) << xc(3) - dx(2)
     712        << endl;
     713 
    729714    QString line;
    730715    int messageType = COTYPE_GPSCOMBINED;
Note: See TracChangeset for help on using the changeset viewer.