diff --git a/Solver/TESTCASES/AdvectionDiffusion.lua b/Solver/TESTCASES/AdvectionDiffusion.lua
index 847b2e9e2ea8be66e4756bf7ea9ed139cebd3f1f..90a85807b86d6b0a7c41c114e30a045a5922c036 100644
--- a/Solver/TESTCASES/AdvectionDiffusion.lua
+++ b/Solver/TESTCASES/AdvectionDiffusion.lua
@@ -1,6 +1,6 @@
 model = GModel  ()
 model:load ('square.geo')
-model:mesh(2)
+model:load ('square.msh')
 dg = dgSystemOfEquations (model)
 dg:setOrder(5)
 
@@ -40,7 +40,7 @@ dg:L2Projection(createFunction.lua(1,'initial_condition','XYZ'))
 dg:exportSolution('output/Advection_00000')
 
 -- main loop
-for i=1,10 do
+for i=1,10000 do
   norm = dg:RK44(0.001)
   if (i % 1 == 0) then 
     print('iter',i,norm)
diff --git a/Solver/TESTCASES/WavePulse.lua b/Solver/TESTCASES/WavePulse.lua
index a4d007405150394d830d82fe476e6b90feec1252..815eb769e9ebb4af100cb057a8d0772b5ca8f7d1 100644
--- a/Solver/TESTCASES/WavePulse.lua
+++ b/Solver/TESTCASES/WavePulse.lua
@@ -8,8 +8,8 @@ function initial_condition( _x , _f )
     X = xyz:get(i,0) - .5
     Y = xyz:get(i,1) - .5
     Z = xyz:get(i,2)
-    VALUE = math.exp(-40*(X*X+Y*Y+Z*Z));   
-    f:set(i,0,7) 
+    VALUE = math.exp(-40*(X*X+Y*Y+Z*Z)); 
+    f:set(i,0,VALUE) 
     f:set(i,1,0.0) 
     f:set(i,2,0.0) 
   end
@@ -21,8 +21,8 @@ end
 
 print'*** Loading the mesh and the model ***'
 myModel   = GModel  ()
-myModel:load('square.geo')
-myModel:mesh(2)
+myModel:load('box.geo')
+myModel:load('box.msh')
 
 print'*** Create a dg solver ***'
 DG = dgSystemOfEquations (myModel)