Comment out PrintToFile Debugging

Author: ax3l

src/particles/distribution/Thermal.H

@@ -89,7 +89,7 @@ namespace distribution
             amrex::ParticleReal r_scale = matched_scale_radius(data);
             amrex::ParticleReal rmin = rin*r_scale;
             amrex::ParticleReal rmax = rout*r_scale;
-            amrex::PrintToFile("equilibrium_params.out") << r_scale << " " << data.Cintensity << "\n";
+            // amrex::PrintToFile("equilibrium_params.out") << r_scale << " " << data.Cintensity << "\n";
 
             // Scale the parameters p1 and p2
             amrex::ParticleReal rt2pi = sqrt(2.0_prt*pi);
@@ -156,10 +156,12 @@ namespace distribution
                 data.cdf1[j+1] = data.f1;
                 data.cdf2[j+1] = data.f2;
                 // write cumulative density to file for debugging
+                /* comment in for debugging
                 amrex::PrintToFile("density1.out") << reval << " " << data.f1 << "\n";
                 amrex::PrintToFile("density2.out") << reval << " " << data.f2 << "\n";
                 amrex::PrintToFile("phi1.out") << reval << " " << data.phi1 << "\n";
                 amrex::PrintToFile("phi2.out") << reval << " " << data.phi2 << "\n";
+                */
             }
 
         }
@@ -213,7 +215,7 @@ namespace distribution
             amrex::ParticleReal Pdensity1 = data.p1*exp(-potential/T1);
             amrex::ParticleReal Pdensity2 = data.p2*exp(-potential/T2);
             amrex::ParticleReal Pdensity_tot = (1.0_prt-w)*Pdensity1 + w*Pdensity2;
-            amrex::PrintToFile("Pdensity.out") << reval << " " << Pdensity_tot << "\n";
+            // amrex::PrintToFile("Pdensity.out") << reval << " " << Pdensity_tot << "\n";
 
             // Apply map to update f1 and f2:
             data.phi1 = phi1;