Index: /trunk/BNC/src/pppModel.cpp
===================================================================
--- /trunk/BNC/src/pppModel.cpp	(revision 7258)
+++ /trunk/BNC/src/pppModel.cpp	(revision 7259)
@@ -428,5 +428,5 @@
     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,4 +461,5 @@
     return 0.0;
   }
+
   return vtec;
 }
@@ -469,5 +470,5 @@
   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,5 +481,5 @@
   }
 
-  _lonS = fmod((_lambdaPP + (epoch - 50400) * (M_PI / 43200)), 2*M_PI);
+  _lonS = fmod((_lambdaPP + (epoch - 50400) * M_PI / 43200), 2*M_PI);
 
   return;