diff options
Diffstat (limited to 'Source/WarpX.cpp')
| -rw-r--r-- | Source/WarpX.cpp | 7 |
1 files changed, 7 insertions, 0 deletions
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index 585369d70..c30eb79d1 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -891,6 +891,13 @@ WarpX::ReadParameters () AMREX_ALWAYS_ASSERT_WITH_MESSAGE(Geom(0).isPeriodic(0) == 0, "The problem must not be periodic in the radial direction"); + // Ensure code aborts if PEC is specified at r=0 for RZ + if (Geom(0).ProbLo(0) == 0){ + AMREX_ALWAYS_ASSERT_WITH_MESSAGE( + WarpX::field_boundary_lo[0] == FieldBoundaryType::None, + "Error : Field boundary at r=0 must be ``none``. \n"); + } + if (maxwell_solver_id == MaxwellSolverAlgo::PSATD) { // Force do_nodal=true (that is, not staggered) and // use same shape factors in all directions, for gathering |