changed saving method

diff --git a/src/domains.cpp b/src/domains.cpp
index 0f46e6285f7914b5f943a14d62fb8559834203bf..6cbb6dfa4100a9462cbbdc9f9fb7418fc9e8a8a6 100644 (file)
--- a/src/domains.cpp
+++ b/src/domains.cpp
@@ -179,24 +179,27 @@ bool check_domsizes(std::vector< std::vector< int > > cd_domsizes, int cd_domain
 }
 
 void print_domains(std::vector< std::vector< int > > pd_domains, std::vector< int > index, std::string fnNdx_) {
-    std::freopen(fnNdx_.c_str(), "w+", stdout);
+    FILE               *fpNdx_;
+    fpNdx_ = std::fopen(fnNdx_.c_str(), "w+");
     int write_count;
     for (int i = 0; i < pd_domains.size(); i++) {
-        std::cout << "[domain_" << i + 1 << "]\n";
+        std::fprintf(fpNdx_, "[domain_%d]\n", i + 1);
         write_count = 0;
         for (int j = 0; j < pd_domains[i].size(); j++) {
             write_count++;
             if (write_count > 20) {
                 write_count -= 20;
-                std::cout << "\n";
+                std::fprintf(fpNdx_, "\n");
             }
-            std::cout << index[pd_domains[i][j]] + 1 << " ";
+            std::fprintf(fpNdx_, "%5d ", index[pd_domains[i][j]] + 1);
         }
-        std::cout << "\n\n";
+        std::fprintf(fpNdx_,"\n\n");
     }
-    std::fclose(stdout);
+    std::fprintf(fpNdx_,"\n");
+    std::fclose(fpNdx_);
 }
 
+
 /*! \brief
  * Class used to compute free volume in a simulations box.
  *