diff --git a/src/main.cpp b/src/main.cpp
index 84cdf276d7ffdf0cc0e6ee0bd640a66f7f15f365..be703b8996425c34a09dcce7f9030495db8d3855 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -155,7 +155,7 @@ int main(int argc, char *argv[])
 
     while((t<tmax) and (iteration<itermax)) {
       //double dt = 0.4*acoustic_solver.acoustic_dt(Vj, cj);
-      double dt = 0.4*finite_volumes_diffusion.diffusion_dt(rhoj, kj);
+      double dt = 0.8*finite_volumes_diffusion.diffusion_dt(rhoj, kj);
 
       if (t+dt>tmax) {
 	dt=tmax-t;
diff --git a/src/mesh/Mesh.hpp b/src/mesh/Mesh.hpp
index 9de8bc118eae4fcf8fac8ad870e4a7f20c5930ee..a5ce53cca535cdffc3cab73058b848de2c8ca9ec 100644
--- a/src/mesh/Mesh.hpp
+++ b/src/mesh/Mesh.hpp
@@ -50,7 +50,7 @@ public:
 
   // pas constant
 
-  
+  /*
   Mesh(const Connectivity& connectivity)
     : m_connectivity(connectivity),
       m_xr("xr", connectivity.numberOfNodes())
@@ -60,11 +60,11 @@ public:
   	m_xr[r][0] = r*delta_x;
       });
   }
-    
+*/
 
   // pas non constant
 
-  /*
+    
   Mesh(const Connectivity& connectivity)
   : m_connectivity(connectivity),
     m_xr("xr", connectivity.numberOfNodes())
@@ -85,7 +85,7 @@ public:
   	}
       });
   }
-  */
+  
 
   ~Mesh()
   {
diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp
index c9aeb725d03d59a1d12bd965345b95f6b1067d1c..04558707e8ee74b1d813f511a459cb56ecac5db1 100644
--- a/src/scheme/FiniteVolumesDiffusion.hpp
+++ b/src/scheme/FiniteVolumesDiffusion.hpp
@@ -127,9 +127,14 @@ private:
 	}                
 	 
 	m_Gl(l) = (1./face_nb_cells(l))*Fl(l)*sum;
-
       });
 
+    // C'est juste un test, mais ca devrait marcher
+    m_Gl(0) = 0;
+    m_Gl(m_mesh.numberOfFaces()-1) = 0;
+    // m_Gl(0) = 0.5*m_Gl(0);
+    // m_Gl(m_mesh.numberOfFaces()-1) = 0.5*m_Gl(m_mesh.numberOfFaces()-1);
+
     return m_Gl ;
   }