diff --git a/src/main.cpp b/src/main.cpp
index a0acdd5f753a05c9c65ab6bbfcfeefb8d894709f..e939a34cd9b15ecbdd43edf7f68fb81cfae78adb 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -225,6 +225,16 @@ int main(int argc, char *argv[])
     //method_cost_map["AcousticSolverWithMesh"] = timer.seconds();
     method_cost_map["FiniteVolumesDiffusionWithMesh"] = timer.seconds();
     
+    { // gnuplot output for density
+     const Kokkos::View<const Rd*> xj   = mesh_data.xj();
+     const Kokkos::View<const Rd*> uj = unknowns.uj();
+     std::ofstream fout("resultat rho");
+     fout.precision(15);
+     for (size_t j=0; j<mesh.numberOfCells(); ++j) {
+       fout << xj[j][0] << ' ' << rhoj[j] << '\n';
+     }
+     }
+    
      { // gnuplot output for vitesse
      const Kokkos::View<const Rd*> xj   = mesh_data.xj();
      const Kokkos::View<const Rd*> uj = unknowns.uj();