Changeset 7259 in ntrip

Timestamp:

Aug 28, 2015, 1:20:35 PM (9 years ago)

Author:

stuerze

Message:

small fix in the conversion from vtec to stec

File:

• trunk/BNC/src/pppModel.cpp (modified) (4 diffs)

trunk/BNC/src/pppModel.cpp

@@ -428 +428 @@
     piercePoint(vTec->_layers[ii]._height, epoch, geocSta.data(), sphEle, sphAzi);
     double vtec = vtecSingleLayerContribution(vTec->_layers[ii]);
-    stec += vtec * sin(sphEle * _psiPP);
+    stec += vtec * sin(sphEle + _psiPP);
   }
   return stec;
@@ -461 +461 @@
     return 0.0;
   }
+
   return vtec;
 }
@@ -469 +470 @@
   double q = (t_CST::rgeoc + geocSta[2]) / (t_CST::rgeoc + layerHeight);
 
-  _psiPP = M_PI / 2 - sphEle - asin(q * cos(sphEle));
+  _psiPP = M_PI/2 - sphEle - asin(q * cos(sphEle));
 
   _phiPP = asin(sin(geocSta[0]) * cos(_psiPP) + cos(geocSta[0]) * sin(_psiPP) * cos(sphAzi));
@@ -480 +481 @@
   }
 
-  _lonS = fmod((_lambdaPP + (epoch - 50400) * (M_PI / 43200)), 2*M_PI);
+  _lonS = fmod((_lambdaPP + (epoch - 50400) * M_PI / 43200), 2*M_PI);
 
   return;