diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
index ce30fc3520ff3d9239a8cbbec3f539fcdfaf6466..d56dc2a5632bfb6f40b9373b6e3ab063f95a30be 100644
--- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
@@ -352,7 +352,7 @@ class RoeViscousFormEulerianCompositeSolver_v2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               Rdxd MatriceProj(identity);
               MatriceProj -= tensorProduct(normal, normal);
@@ -459,7 +459,7 @@ class RoeViscousFormEulerianCompositeSolver_v2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
@@ -579,7 +579,7 @@ class RoeViscousFormEulerianCompositeSolver_v2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
index d22e940a9c3d53575747cf3d09bb1160377f209a..e707e3f8474387385fcb9a5d3a0b9b1f2c97054d 100644
--- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
+++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp
@@ -360,7 +360,7 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               Rdxd MatriceProj(identity);
               MatriceProj -= tensorProduct(normal, normal);
@@ -467,7 +467,7 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
@@ -587,7 +587,7 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
diff --git a/src/scheme/RusanovEulerianCompositeSolver.cpp b/src/scheme/RusanovEulerianCompositeSolver.cpp
index a13b34fb94582505b86683af9eed5efdca6f8af3..cda7e625cdbc5f83542124d51d1f3fbfd4ad784b 100644
--- a/src/scheme/RusanovEulerianCompositeSolver.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver.cpp
@@ -245,7 +245,7 @@ class RusanovEulerianCompositeSolver
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
@@ -481,7 +481,7 @@ class RusanovEulerianCompositeSolver
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               Rdxd MatriceProj(identity);
               MatriceProj -= tensorProduct(normal, normal);
@@ -588,7 +588,7 @@ class RusanovEulerianCompositeSolver
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
diff --git a/src/scheme/RusanovEulerianCompositeSolver_o2.cpp b/src/scheme/RusanovEulerianCompositeSolver_o2.cpp
index 867e6c4e5f7c2f9865820e954c26c219800f2cb9..7d53c2f884e925b3bba1c442c62ee3f13fd6533d 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_o2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_o2.cpp
@@ -253,7 +253,7 @@ class RusanovEulerianCompositeSolver_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
@@ -489,7 +489,7 @@ class RusanovEulerianCompositeSolver_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               Rdxd MatriceProj(identity);
               MatriceProj -= tensorProduct(normal, normal);
@@ -596,7 +596,7 @@ class RusanovEulerianCompositeSolver_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
index cde234af1fcfe35fa66938cbc30b5a37942f3392..f367cb161b634b9b2a4eff54fa0dcfc2387e2376 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp
@@ -246,7 +246,7 @@ class RusanovEulerianCompositeSolver_v2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
@@ -483,7 +483,7 @@ class RusanovEulerianCompositeSolver_v2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               Rdxd MatriceProj(identity);
               MatriceProj -= tensorProduct(normal, normal);
@@ -590,7 +590,7 @@ class RusanovEulerianCompositeSolver_v2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp
index 7386f4e84451e2d1b99d284df97eb5f9f44693d1..7a22ea9c73f324fd65fe5598d386a51fce0e5695 100644
--- a/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp
+++ b/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp
@@ -253,7 +253,7 @@ class RusanovEulerianCompositeSolver_v2_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;
 
@@ -488,7 +488,7 @@ class RusanovEulerianCompositeSolver_v2_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               Rdxd MatriceProj(identity);
               MatriceProj -= tensorProduct(normal, normal);
@@ -595,7 +595,7 @@ class RusanovEulerianCompositeSolver_v2_o2
 
               Rd vectorSym(zero);
               for (size_t dim = 0; dim < Dimension; ++dim)
-                vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim];
+                vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim];
 
               vectorSym -= dot(vectorSym, normal) * normal;