aboutsummaryrefslogtreecommitdiff
path: root/Source/Utils/WarpXMovingWindow.cpp
diff options
context:
space:
mode:
authorGravatar Remi Lehe <remi.lehe@normalesup.org> 2019-10-28 16:12:50 -0700
committerGravatar GitHub <noreply@github.com> 2019-10-28 16:12:50 -0700
commitf326192ee5140278c6061b1171b98c6a0f932320 (patch)
treeb9e687ace23d47a301e2ac07123bd41b912590ad /Source/Utils/WarpXMovingWindow.cpp
parent2eac8727874e11c3c8163b3d9ed1ec8854717753 (diff)
parentca6e70897035190dfdd2e706163b73bd383775b3 (diff)
downloadWarpX-f326192ee5140278c6061b1171b98c6a0f932320.tar.gz
WarpX-f326192ee5140278c6061b1171b98c6a0f932320.tar.zst
WarpX-f326192ee5140278c6061b1171b98c6a0f932320.zip
Merge pull request #476 from MaxThevenet/initial_fields
Add uniform initial field to the grid
Diffstat (limited to 'Source/Utils/WarpXMovingWindow.cpp')
-rw-r--r--Source/Utils/WarpXMovingWindow.cpp13
1 files changed, 7 insertions, 6 deletions
diff --git a/Source/Utils/WarpXMovingWindow.cpp b/Source/Utils/WarpXMovingWindow.cpp
index f84701a02..c577da7f3 100644
--- a/Source/Utils/WarpXMovingWindow.cpp
+++ b/Source/Utils/WarpXMovingWindow.cpp
@@ -99,8 +99,8 @@ WarpX::MoveWindow (bool move_j)
for (int dim = 0; dim < 3; ++dim) {
// Fine grid
- shiftMF(*Bfield_fp[lev][dim], geom[lev], num_shift, dir);
- shiftMF(*Efield_fp[lev][dim], geom[lev], num_shift, dir);
+ shiftMF(*Bfield_fp[lev][dim], geom[lev], num_shift, dir, B_external_grid[dim]);
+ shiftMF(*Efield_fp[lev][dim], geom[lev], num_shift, dir, E_external_grid[dim]);
if (move_j) {
shiftMF(*current_fp[lev][dim], geom[lev], num_shift, dir);
}
@@ -113,8 +113,8 @@ WarpX::MoveWindow (bool move_j)
if (lev > 0) {
// Coarse grid
- shiftMF(*Bfield_cp[lev][dim], geom[lev-1], num_shift_crse, dir);
- shiftMF(*Efield_cp[lev][dim], geom[lev-1], num_shift_crse, dir);
+ shiftMF(*Bfield_cp[lev][dim], geom[lev-1], num_shift_crse, dir, B_external_grid[dim]);
+ shiftMF(*Efield_cp[lev][dim], geom[lev-1], num_shift_crse, dir, E_external_grid[dim]);
shiftMF(*Bfield_aux[lev][dim], geom[lev], num_shift, dir);
shiftMF(*Efield_aux[lev][dim], geom[lev], num_shift, dir);
if (move_j) {
@@ -203,7 +203,8 @@ WarpX::MoveWindow (bool move_j)
}
void
-WarpX::shiftMF (MultiFab& mf, const Geometry& geom, int num_shift, int dir)
+WarpX::shiftMF (MultiFab& mf, const Geometry& geom, int num_shift, int dir,
+ amrex::Real external_field)
{
BL_PROFILE("WarpX::shiftMF()");
const BoxArray& ba = mf.boxArray();
@@ -257,7 +258,7 @@ WarpX::shiftMF (MultiFab& mf, const Geometry& geom, int num_shift, int dir)
if (outbox.ok()) {
AMREX_PARALLEL_FOR_4D ( outbox, nc, i, j, k, n,
{
- srcfab(i,j,k,n) = 0.0;
+ srcfab(i,j,k,n) = external_field;
});
}