diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp
index cf2d371586dd46a73a0f1028fb6478d9e8b95943..1db575470259d965c6362a9d35c345b29ddb15ba 100644
--- a/src/scheme/FiniteVolumesEulerUnknowns.hpp
+++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp
@@ -227,8 +227,8 @@ void initializeSod()
 
     Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){
 	//m_Ej[j] = m_ej[j]+0.5*(m_uj[j],m_uj[j]);
-	m_Ej[j] = 0.5*(m_uj[j],m_uj[j]);
-	//m_Ej[j] = 2.;
+	//m_Ej[j] = (std::cos(pi*xj[j][0]))*(std::cos(pi*xj[j][0]))*std::sin(pi*xj[j][0])*((pi*pi*pi)*0.5);
+	m_Ej[j] = 2.;
       });
 
     const Kokkos::View<const double*> Vj = m_mesh_data.Vj();