diff --git a/dart/pyVII/vCoupler.py b/dart/pyVII/vCoupler.py
index 449f2a2ce2c669f6c4e5744220c3d2a5eae57e03..e0d5813bfb4f6b51c38e086d2567b7cf5db612e1 100644
--- a/dart/pyVII/vCoupler.py
+++ b/dart/pyVII/vCoupler.py
@@ -89,7 +89,7 @@ class Coupler:
       error = abs((self.vSolver.Cdt - cdPrev)/self.vSolver.Cdt) if self.vSolver.Cdt != 0 else 1
       
       if error  <= self.couplTol and exitCode==0:
-
+        print('{:>5s}| {:>10s} {:>10s} | {:>8s} {:>6s} {:>8s}'.format('It', 'Cl', 'Cd', 'xtrT', 'xtrB', 'Error'))
         print(ccolors.ANSI_GREEN,'{:>4.0f}| {:>10.5f} {:>10.5f} | {:>8.2f} {:>6.2f} {:>8.2f}\n'.format(couplIter,
         self.iSolver.Cl, self.vSolver.Cdt, xtr[0], xtr[1], np.log10(error)),ccolors.ANSI_RESET)
         return 0
diff --git a/dart/src/wViscSolver.cpp b/dart/src/wViscSolver.cpp
index 586cdacfd4e989d1ff2f9c82a660464871ca8cf0..84f5388d6e5d7bc8a42e77456a3bb991317f9eb0 100644
--- a/dart/src/wViscSolver.cpp
+++ b/dart/src/wViscSolver.cpp
@@ -220,6 +220,11 @@ int ViscSolver::Run(unsigned int couplIter)
 
                 pointExitCode = tSolver->Integration(iPoint, Sections[iSec][iRegion]);
 
+                if (Sections[iSec][iRegion]->Cf[iPoint] < 0)
+                {
+                    std::cout << "Negative Cf at point " << iPoint << " : " << Sections[iSec][iRegion]->Cf[iPoint] << std::endl;
+                }
+
                 /* if (iRegion == 1 && iPoint == 265)
                 {
                 std::cout << " iSec " << iSec
diff --git a/dart/tests/bli_lowLift.py b/dart/tests/bli_lowLift.py
index 78b802f30cc977d041c3eeb115481fb5b8ead3bf..eddabfa0ad56bdb97d12f89d00dda8498c442ff6 100644
--- a/dart/tests/bli_lowLift.py
+++ b/dart/tests/bli_lowLift.py
@@ -57,9 +57,9 @@ def main():
 
     # define flow variables
     Re = 1e7
-    alpha = 12*math.pi/180
+    alpha = 2*math.pi/180
     #alphaSeq = np.arange(-5, 5.5, 0.5)
-    M_inf = 0.
+    M_inf = 0.715
 
     meshFactor = 1
     CFL0 = 1
@@ -76,8 +76,8 @@ def main():
     # mesh the geometry
     print(ccolors.ANSI_BLUE + 'PyMeshing...' + ccolors.ANSI_RESET)
     tms['msh'].start()
-    pars = {'xLgt' : 5, 'yLgt' : 5, 'msF' : 1, 'msTe' : 0.01, 'msLe' : 0.001}
-    msh, gmshWriter = floD.mesh(dim, 'models/n0012.geo', pars, ['field', 'airfoil', 'downstream'])
+    pars = {'xLgt' : 5, 'yLgt' : 5, 'msF' : 1, 'msTe' : 0.01, 'msLe' : 0.005}
+    msh, gmshWriter = floD.mesh(dim, 'models/rae2822.geo', pars, ['field', 'airfoil', 'downstream'])
     tms['msh'].stop()
 
     # set the problem and add medium, initial/boundary conditions