diff --git a/src/main.cpp b/src/main.cpp
index 0a83b28d90aa0cb3a958169283c459a3020bb299..6325e25b9ce7f4d9c4df53411a87afe44b36db32 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -160,7 +160,7 @@ int main(int argc, char *argv[])
     double c = 0.;
     c = finite_volumes_diffusion.conservatif(unknowns);
 
-    /*
+    
     // Ecriture des valeurs initiales dans un fichier
     
     const Kokkos::View<const Rd*> xj = mesh_data.xj();
@@ -257,7 +257,7 @@ int main(int argc, char *argv[])
       }
     }
     diff.close();
-    */
+    
 
     while((t<tmax) and (iteration<itermax)) {
      
@@ -325,7 +325,7 @@ int main(int argc, char *argv[])
       ++iteration;
       std::cout << "temps t : " << t << std::endl;
 
-      /*
+      
       // ECRITURE DANS UN FICHIER 
       
       if ((std::fmod(t,0.01) < 0.0001) or (t == tmax)) {
@@ -544,7 +544,7 @@ int main(int argc, char *argv[])
       riffout.close();
       
       }
-      */
+      
       // ENTROPY TEST
       //finite_volumes_diffusion.entropie(unknowns);
       
diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp
index b1abbde581e6045076a371a33ef3ee6758d27b9e..98e5859ec03b192c467654bfaa7c23f6a58d5bd5 100644
--- a/src/scheme/FiniteVolumesEulerUnknowns.hpp
+++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp
@@ -279,7 +279,8 @@ public:
   	  m_kj[j]=0.;
   	} else {
 	  if (xj[j][0]<0.9){
-	    m_kj[j]=0.05;
+	    // m_kj[j]=0.05;
+	    m_kj[j] = 0.0028; // Re = 500
 	  } else {
 	    m_kj[j]=0. ;
 	  }