diff --git a/Solver/TESTCASES/Stommel.lua b/Solver/TESTCASES/Stommel.lua
index 97539ed9af6b8540e2f772a5f0df541605e0ddc0..741048dbf76e8db8b9ca8ac02a1364b658de95af 100644
--- a/Solver/TESTCASES/Stommel.lua
+++ b/Solver/TESTCASES/Stommel.lua
@@ -1,6 +1,6 @@
 model = GModel()
 model:load ('stommel_square.msh')
-order=3
+order=1
 dimension = 2
 claw = dgConservationLawShallowWater2d()
 claw:addBoundaryCondition('Wall',claw:newBoundaryWall())
@@ -8,17 +8,11 @@ groups = dgGroupCollection(model, dimension, order)
 -- groups:split...
 groups:buildGroupsOfInterfaces()
 solution = dgDofContainer(groups, claw:getNbFields())
-a = fullMatrix(3,1);
-a:set(0,0, 1.);
-a:set(1,0, 2.);
-a:set(2,0, 3.);
-f = functionConstant(a):getName()
-solution:L2Projection(f)
 solution:exportMsh('output/init')
-
+rk=dgRungeKutta()
 
 for i=1,100000 do
-  norm = dg:RK44(150*(3/(2.*order+1)))
+  norm = rk:iterate33(claw,150*(3/(2.*order+1)/2),solution)
   if ( i%100 ==0 ) then
     print ('iter ', i, norm)
   end