diff --git a/Solver/TESTCASES/CylinderEddies.lua b/Solver/TESTCASES/CylinderEddies.lua
index 0a7d9cf3c5fd8cee34c834c6c4ee47b9d1dc0465..92987c61d022ac23b26d132629cfa003ab93186a 100644
--- a/Solver/TESTCASES/CylinderEddies.lua
+++ b/Solver/TESTCASES/CylinderEddies.lua
@@ -1,7 +1,7 @@
 MACH = .05;
 RHO  = 1.0;
 PRES = 1./(MACH*RHO*RHO*1.4*1.4) 
-V = 1.0 
+V = 1 
 SOUND = V/MACH
 
 --[[ 
@@ -26,7 +26,7 @@ FS = functionLua(4, 'free_stream', {'XYZ'}):getName()
 
 -- diffusivity
 mu=fullMatrix(1,1);
-mu:set(0,0,0.01)
+mu:set(0,0,0.025)
 kappa=fullMatrix(1,1);
 kappa:set(0,0,0.01)
 
@@ -62,7 +62,7 @@ CFL = 20.1;
 dt = CFL * DG:computeInvSpectralRadius();
 print('DT = ',dt)
 T = 0;
-for i=1,1000000 do
+for i=1,2 do
     dt = CFL * DG:computeInvSpectralRadius();    
     norm = DG:RK44(dt)
     T = T + dt
diff --git a/Solver/TESTCASES/cyl.geo b/Solver/TESTCASES/cyl.geo
index 5779ec056692987e17dcc56ab56dfa994f7268d8..8b8595b31e20050136db0975d27b9f94b7b9c466 100644
--- a/Solver/TESTCASES/cyl.geo
+++ b/Solver/TESTCASES/cyl.geo
@@ -73,3 +73,4 @@ Physical Surface("BoundaryLayer") = {100};
 Physical Line("Cylinder") = {4, 1, 2, 3};
 Physical Line("Box") = {5, 6, 7, 8};
 Mesh.CharacteristicLengthExtendFromBoundary=1;
+Recombine Surface {4, 1, 2, 3};