diff --git a/src/mesh/Mesh.hpp b/src/mesh/Mesh.hpp
index 1061d0afd2fe1fe02edfe77358a032ee9f91ff67..b2a8c9b102df7fc9320134c9b0ce70ad89f33caa 100644
--- a/src/mesh/Mesh.hpp
+++ b/src/mesh/Mesh.hpp
@@ -57,7 +57,7 @@ public:
   // const double delta_x = 1./connectivity.numberOfCells();
   // Kokkos::parallel_for(connectivity.numberOfNodes(), KOKKOS_LAMBDA(const int& r){
   //	m_xr[r][0] = r*delta_x;
-  //  });
+  // });
   //}
 
   // pas non constant
@@ -79,13 +79,16 @@ public:
 
   // pas non constant avec fonction x^2
 
-   Mesh(const Connectivity& connectivity)
-   : m_connectivity(connectivity),
-     m_xr("xr", connectivity.numberOfNodes())
+  Mesh(const Connectivity& connectivity)
+    : m_connectivity(connectivity),
+      m_xr("xr", connectivity.numberOfNodes())
   {
-   const double delta_x = 1./connectivity.numberOfCells();
-  Kokkos::parallel_for(connectivity.numberOfNodes(), KOKKOS_LAMBDA(const int& r){
-  	m_xr[r][0] = (r*delta_x)*(r*delta_x);
+    const double delta_x = 1./connectivity.numberOfCells();
+    Kokkos::parallel_for(connectivity.numberOfNodes(), KOKKOS_LAMBDA(const int& r){
+  	m_xr[r][0] = r*delta_x;
+      });
+    Kokkos::parallel_for(connectivity.numberOfNodes(), KOKKOS_LAMBDA(const int& r){
+    	m_xr[r][0] = m_xr[r][0]*m_xr[r][0];
     });
   }
 
diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp
index adc46f0417c5a4f3dba6c2e4fc7c2b2eb8661b13..789dfa04eb99ca4739c28ac689cf31a938d4b058 100644
--- a/src/scheme/FiniteVolumesDiffusion.hpp
+++ b/src/scheme/FiniteVolumesDiffusion.hpp
@@ -100,9 +100,8 @@ private:
 	}
 
 	// k = x
-	                      
-	m_Fl(l) = ((sum2*0.5)/Vl(l))*sum;
-
+	// m_Fl(l) = ((sum2*0.5)/Vl(l))*sum;
+	m_Fl(l) = ((sum2/face_nb_cells(l))/Vl(l))*sum;
         // k = 2
 	 
 	//m_Fl(l)= (2./Vl(l))*sum;