diff options
Diffstat (limited to 'Source/Evolve')
| -rw-r--r-- | Source/Evolve/WarpXComputeDt.cpp | 9 | ||||
| -rw-r--r-- | Source/Evolve/WarpXEvolve.cpp | 1 |
2 files changed, 7 insertions, 3 deletions
diff --git a/Source/Evolve/WarpXComputeDt.cpp b/Source/Evolve/WarpXComputeDt.cpp index 392ef81b5..08d12dde6 100644 --- a/Source/Evolve/WarpXComputeDt.cpp +++ b/Source/Evolve/WarpXComputeDt.cpp @@ -38,7 +38,9 @@ WarpX::ComputeDt () if (maxwell_solver_id == MaxwellSolverAlgo::PSATD) { // Computation of dt for spectral algorithm // (determined by the minimum cell size in all directions) -#if (AMREX_SPACEDIM == 2) +#if (AMREX_SPACEDIM == 1) + deltat = cfl * dx[0] / PhysConst::c; +#elif (AMREX_SPACEDIM == 2) deltat = cfl * std::min(dx[0], dx[1]) / PhysConst::c; #else deltat = cfl * std::min(dx[0], std::min(dx[1], dx[2])) / PhysConst::c; @@ -84,10 +86,13 @@ WarpX::PrintDtDxDyDz () for (int lev=0; lev <= max_level; lev++) { const amrex::Real* dx_lev = geom[lev].CellSize(); amrex::Print() << "Level " << lev << ": dt = " << dt[lev] +#if (defined WARPX_DIM_1D_Z) + << " ; dz = " << dx_lev[0] << '\n'; +#elif (defined WARPX_DIM_XZ) || (defined WARPX_DIM_RZ) << " ; dx = " << dx_lev[0] -#if (defined WARPX_DIM_XZ) || (defined WARPX_DIM_RZ) << " ; dz = " << dx_lev[1] << '\n'; #elif (defined WARPX_DIM_3D) + << " ; dx = " << dx_lev[0] << " ; dy = " << dx_lev[1] << " ; dz = " << dx_lev[2] << '\n'; #endif diff --git a/Source/Evolve/WarpXEvolve.cpp b/Source/Evolve/WarpXEvolve.cpp index d04a54ae3..770887cfe 100644 --- a/Source/Evolve/WarpXEvolve.cpp +++ b/Source/Evolve/WarpXEvolve.cpp @@ -253,7 +253,6 @@ WarpX::Evolve (int numsteps) bool move_j = is_synchronized; // If is_synchronized we need to shift j too so that next step we can evolve E by dt/2. // We might need to move j because we are going to make a plotfile. - int num_moved = MoveWindow(step+1, move_j); mypc->ContinuousFluxInjection(dt[0]); |