diff --git a/flow/viscous/solver.py b/flow/viscous/solver.py
index c6b5a9ccc8d15cec12657224b46bc9eb2817b8d9..002f585a8e23c3b41b8c469701864ca7ac97db05 100644
--- a/flow/viscous/solver.py
+++ b/flow/viscous/solver.py
@@ -20,6 +20,14 @@ class Solver():
     def run(self):
         ''' Solve BLE
         '''
+        # define tangential velocity in 2D
+        vt = self.v[:,0] + self.v[:,1]
+        # define stagnation point 
+        stag = np.argmin(abs(self.v[:,1])) #define the stagnation point 
+        print 'Stag is' 
+        print stag
+        print len(self.boundary.nodes)
+
         print '--- I am a fake BL solver and I am setting dummy blowing velocities! ---'
         # 6th order fitting of delta_star for naca 0012 at alpha = 0, from XFoil
         a = 0.000054925976379640116
@@ -32,4 +40,4 @@ class Solver():
         for i in range(0, self.nE):
             x = 0.5*(self.boundary.groups[0].tag.elems[i].nodes[0].pos.x[0]+self.boundary.groups[0].tag.elems[i].nodes[1].pos.x[0])
             # CAUTION: a positive velocity will result in a suction!
-            self.u[i] = -(b + 2*c*x + 3*d*x*x* + 4*e*x*x*x + 5*f*x*x*x*x + 6*g*x*x*x*x*x) # ue = d(delta_star) / dx
+            self.u[i] = -(b + 2*c*x + 3*d*x*x + 4*e*x*x*x + 5*f*x*x*x*x + 6*g*x*x*x*x*x) # ue = d(delta_star) / dx