diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
index 2d4c3b74a82f39f3e1db0a29d8a68d3162f34d1c..2aab892d8df4533b7d2aef934add050cfe926dce 100644
--- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
@@ -1641,17 +1641,17 @@ class RoeViscousFormEulerianCompositeSolver_v2
     }   // dim 3
 
     // Pour les assemblages
-    double theta = .5;
-    double eta   = .2;
+    double theta = 2. / 3.;   //.5;
+    double eta   = 1. / 6.;   //.2;
     if constexpr (Dimension == 2) {
       eta = 0;
-    } else {
-      theta = 1. / 3.;
-      eta   = 1. / 3.;
-
-      // theta = .5;
-      // eta   = 0;
     }
+    // else {
+    // theta = 1. / 3.;
+    // eta   = 1. / 3.;
+    // theta = .5;
+    // eta   = 0;
+    //}
     //
     parallel_for(
       p_mesh->numberOfCells(), PUGS_LAMBDA(CellId j) {
diff --git a/src/scheme/RusanovEulerianCompositeSolver.cpp b/src/scheme/RusanovEulerianCompositeSolver.cpp
index 4192bacc4a677921c7151e801f17d39d02b65032..d8db36aa23080b4c29b2ea3e60910964a8f14520 100644
--- a/src/scheme/RusanovEulerianCompositeSolver.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver.cpp
@@ -1218,17 +1218,17 @@ class RusanovEulerianCompositeSolver
     }   // dim 3
 
     // Pour les assemblages
-    double theta = .5;
-    double eta   = .2;
+    double theta = 2. / 3.;   //.5;
+    double eta   = 1. / 6.;   //.2;
     if constexpr (Dimension == 2) {
       eta = 0;
-    } else {
-      theta = 1. / 3.;
-      eta   = 1. / 3.;
-
-      // theta = .5;
-      // eta   = 0;
     }
+    // else {
+    // theta = 1. / 3.;
+    // eta   = 1. / 3.;
+    // theta = .5;
+    // eta   = 0;
+    //}
     //
     parallel_for(
       p_mesh->numberOfCells(), PUGS_LAMBDA(CellId j) {
diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
index 1a003ade61c35f99fc8bbb6b251b478999f14dc1..d490202f5f433aaa13867f97d30134f38d2be5fd 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
@@ -1116,17 +1116,17 @@ class RusanovEulerianCompositeSolver_v2
     }   // dim 3
 
     // Pour les assemblages
-    double theta = .5;
-    double eta   = .2;
+    double theta = 2. / 3.;   //.5;
+    double eta   = 1. / 6.;   //.2;
     if constexpr (Dimension == 2) {
       eta = 0;
-    } else {
-      theta = 1. / 3.;
-      eta   = 1. / 3.;
-
-      // theta = .5;
-      // eta   = 0;
     }
+    // else{
+    // theta = 1. / 3.;
+    // eta   = 1. / 3.;
+    // theta = .5;
+    // eta   = 0;
+    //}
     //
     parallel_for(
       p_mesh->numberOfCells(), PUGS_LAMBDA(CellId j) {