Changeset 7259 in ntrip


Ignore:
Timestamp:
Aug 28, 2015, 1:20:35 PM (9 years ago)
Author:
stuerze
Message:

small fix in the conversion from vtec to stec

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/pppModel.cpp

    r7256 r7259  
    428428    piercePoint(vTec->_layers[ii]._height, epoch, geocSta.data(), sphEle, sphAzi);
    429429    double vtec = vtecSingleLayerContribution(vTec->_layers[ii]);
    430     stec += vtec * sin(sphEle * _psiPP);
     430    stec += vtec * sin(sphEle + _psiPP);
    431431  }
    432432  return stec;
     
    461461    return 0.0;
    462462  }
     463
    463464  return vtec;
    464465}
     
    469470  double q = (t_CST::rgeoc + geocSta[2]) / (t_CST::rgeoc + layerHeight);
    470471
    471   _psiPP = M_PI / 2 - sphEle - asin(q * cos(sphEle));
     472  _psiPP = M_PI/2 - sphEle - asin(q * cos(sphEle));
    472473
    473474  _phiPP = asin(sin(geocSta[0]) * cos(_psiPP) + cos(geocSta[0]) * sin(_psiPP) * cos(sphAzi));
     
    480481  }
    481482
    482   _lonS = fmod((_lambdaPP + (epoch - 50400) * (M_PI / 43200)), 2*M_PI);
     483  _lonS = fmod((_lambdaPP + (epoch - 50400) * M_PI / 43200), 2*M_PI);
    483484
    484485  return;
Note: See TracChangeset for help on using the changeset viewer.