2 files changed, 4 insertions, 7 deletions
diff --git a/Source/FieldSolver/SpectralSolver/SpectralSolver.H b/Source/FieldSolver/SpectralSolver/SpectralSolver.H
index bdec4499d..91cfc84f1 100644
--- a/Source/FieldSolver/SpectralSolver/SpectralSolver.H
+++ b/Source/FieldSolver/SpectralSolver/SpectralSolver.H
@@ -36,7 +36,8 @@ class SpectralSolver
                         const amrex::RealVect dx, const amrex::Real dt,
                         const bool pml=false,
                         const bool periodic_single_box=false,
-                        const bool update_with_rho=false );
+                        const bool update_with_rho=false,
+                        const bool fft_do_time_averaging=false);
 
         /**
          * \brief Transform the component `i_comp` of MultiFab `mf`
@@ -83,8 +84,6 @@ class SpectralSolver
              algorithm->CurrentCorrection( field_data, current, rho );
         };
 
-        bool fft_do_time_averaging = false;
-
     private:
         void ReadParameters ();
 
diff --git a/Source/FieldSolver/SpectralSolver/SpectralSolver.cpp b/Source/FieldSolver/SpectralSolver/SpectralSolver.cpp
index 16894da78..7146c7d4c 100644
--- a/Source/FieldSolver/SpectralSolver/SpectralSolver.cpp
+++ b/Source/FieldSolver/SpectralSolver/SpectralSolver.cpp
@@ -39,7 +39,8 @@ SpectralSolver::SpectralSolver(
                 const amrex::Array<amrex::Real,3>& v_galilean,
                 const amrex::RealVect dx, const amrex::Real dt,
                 const bool pml, const bool periodic_single_box,
-                const bool update_with_rho ) {
+                const bool update_with_rho,
+                const bool fft_do_time_averaging) {
 
     // Initialize all structures using the same distribution mapping dm
 
@@ -51,9 +52,6 @@ SpectralSolver::SpectralSolver(
     // - Select the algorithm depending on the input parameters
     //   Initialize the corresponding coefficients over k space
 
-    amrex::ParmParse pp("psatd");
-    pp.query("do_time_averaging", fft_do_time_averaging);
-
     if (pml) {
         algorithm = std::unique_ptr<PMLPsatdAlgorithm>( new PMLPsatdAlgorithm(
             k_space, dm, norder_x, norder_y, norder_z, nodal, dt ) );