diff --git a/appl/coupling-ff-pm/common/outputhelper.hh b/appl/coupling-ff-pm/common/outputhelper.hh
index d414edab42881ce24db48ee09ef5c8183937d868..51baab53613b8095a1bbc63674fe5931324f7892 100644
--- a/appl/coupling-ff-pm/common/outputhelper.hh
+++ b/appl/coupling-ff-pm/common/outputhelper.hh
@@ -52,7 +52,7 @@ std::tuple<double,double,double> writeVelocitiesOnInterfaceToFile( const std::st
         const int prec = ofs.precision();
         ofs << std::setprecision(std::numeric_limits<double>::digits10 + 1) << v << "\n";
         ofs.precision( prec );
-        ofs << v / 1e3 << "\n";
+//        ofs << v / 1e3 << "\n";
       }
     }
   }
diff --git a/appl/coupling-ff-pm/fvca-monolithic-reversed/fvca-monolithic.cc b/appl/coupling-ff-pm/fvca-monolithic-reversed/fvca-monolithic.cc
index 7bddbfc30e84f40e67719f413e69aad753ce2e9f..48ae76dca2222a90376ed484ce00ca15f7733e63 100644
--- a/appl/coupling-ff-pm/fvca-monolithic-reversed/fvca-monolithic.cc
+++ b/appl/coupling-ff-pm/fvca-monolithic-reversed/fvca-monolithic.cc
@@ -151,9 +151,9 @@ int main(int argc, char** argv) try
     auto darcyGridVariables = std::make_shared<DarcyGridVariables>(darcyProblem, darcyFvGridGeometry);
     darcyGridVariables->init(sol[darcyIdx]);
 
-//    couplingManager->setGridVariables(std::make_tuple(stokesGridVariables->cellCenterGridVariablesPtr(),
-//                                                      stokesGridVariables->faceGridVariablesPtr(),
-//                                                      darcyGridVariables));
+    couplingManager->setGridVariables(std::make_tuple(stokesGridVariables->cellCenterGridVariablesPtr(),
+                                                      stokesGridVariables->faceGridVariablesPtr(),
+                                                      darcyGridVariables));
 
     // intialize the vtk output module
     StaggeredVtkOutputModule<StokesGridVariables, decltype(stokesSol)> stokesVtkWriter(*stokesGridVariables, stokesSol, stokesProblem->name());
@@ -249,13 +249,26 @@ int main(int argc, char** argv) try
                                                                                 *darcyProblem,
                                                                                 *darcyGridVariables,
                                                                                 sol[darcyIdx] );
-      const int prec = std::cout.precision();
-      std::cout << "Velocity statistics:" << std::endl
+      const auto prec = std::cout.precision();
+      std::cout << "Velocity statistics (Darcy):" << std::endl
                 << std::setprecision(std::numeric_limits<double>::digits10 + 1)
                 << "  min: " << min << std::endl
                 << "  max: " << max << std::endl
                 << "  sum: " << sum << std::endl;
       std::cout.precision( prec );
+
+      {
+        const std::string filenameDarcy="darcy-flow-statistics.txt";
+        std::ofstream ofs( filenameDarcy+".txt", std::ofstream::out | std::ofstream::trunc);
+        const auto prec = ofs.precision();
+        ofs << "Velocity statistics (Darcy):" << std::endl
+                  << std::setprecision(std::numeric_limits<double>::digits10 + 1)
+                  << "  min: " << min << std::endl
+                  << "  max: " << max << std::endl
+                  << "  sum: " << sum << std::endl;
+        ofs.precision( prec );
+        ofs.close();
+      }
     }
 
     //TODO make freeflow
@@ -268,13 +281,25 @@ int main(int argc, char** argv) try
                                                                        *couplingManager,
                                                                        *stokesProblem,
                                                                        sol[stokesFaceIdx] );
-      const int prec = std::cout.precision();
-      std::cout << "Velocity statistics:" << std::endl
+      const auto prec = std::cout.precision();
+      std::cout << "Velocity statistics (free flow):" << std::endl
                 << std::setprecision(std::numeric_limits<double>::digits10 + 1)
                 << "  min: " << min << std::endl
                 << "  max: " << max << std::endl
                 << "  sum: " << sum << std::endl;
       std::cout.precision( prec );
+      {
+        const std::string filenameFlow="free-flow-statistics.txt";
+        std::ofstream ofs( filenameFlow+".txt", std::ofstream::out | std::ofstream::trunc);
+        const auto prec = ofs.precision();
+        ofs << "Velocity statistics (free flow):" << std::endl
+            << std::setprecision(std::numeric_limits<double>::digits10 + 1)
+            << "  min: " << min << std::endl
+            << "  max: " << max << std::endl
+            << "  sum: " << sum << std::endl;
+        ofs.precision( prec );
+        ofs.close();
+      }
     }
 
     ////////////////////////////////////////////////////////////