diff options
| author |  RevathiJambunathan <revanathan@gmail.com>
| 2019-12-13 13:42:00 -0800 |
|---|---|---|
| committer | RevathiJambunathan <revanathan@gmail.com>
| 2019-12-13 13:42:00 -0800 |
| commit | db6a476daa9d8547b96436d5cc9b5ef5ab212a43 (patch) | |
| tree | fbc37e424e0737af5083746fab959f09f1786131 /Source/Initialization/WarpXInitData.cpp | |
| parent | d430c8236fbca43475cb371205c3b4b030970fab (diff) | |
| download | WarpX-db6a476daa9d8547b96436d5cc9b5ef5ab212a43.tar.gz WarpX-db6a476daa9d8547b96436d5cc9b5ef5ab212a43.tar.zst WarpX-db6a476daa9d8547b96436d5cc9b5ef5ab212a43.zip | |
for 2D, the function is (x,y=0.0,z)
Diffstat (limited to 'Source/Initialization/WarpXInitData.cpp')
| -rw-r--r-- | Source/Initialization/WarpXInitData.cpp | 25 |
1 files changed, 16 insertions, 9 deletions
diff --git a/Source/Initialization/WarpXInitData.cpp b/Source/Initialization/WarpXInitData.cpp index ce11fa0fd..488c351dc 100644 --- a/Source/Initialization/WarpXInitData.cpp +++ b/Source/Initialization/WarpXInitData.cpp @@ -596,13 +596,14 @@ WarpX::InitializeExternalFieldsOnGridUsingParser ( // Shift required in the x-, y-, or z- position // depending on the index type of the multifab Real fac_x = (1.0 - mfx_type[0]) * dx_lev[0]*0.5; - Real fac_y = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; - Real x = i*dx_lev[0] + real_box.lo(0) + fac_x; - Real y = j*dx_lev[1] + real_box.lo(1) + fac_y; #if (AMREX_SPACEDIM==2) - Real z = 0.0; + Real y = 0.0; + Real fac_z = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; + Real z = j*dx_lev[1] + real_box.lo(1) + fac_z; #else + Real fac_y = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; + Real y = j*dx_lev[1] + real_box.lo(1) + fac_y; Real fac_z = (1.0 - mfx_type[2]) * dx_lev[2]*0.5; Real z = k*dx_lev[2] + real_box.lo(2) + fac_z; #endif @@ -611,13 +612,15 @@ WarpX::InitializeExternalFieldsOnGridUsingParser ( }, [=] AMREX_GPU_DEVICE (int i, int j, int k) { Real fac_x = (1.0 - mfy_type[0]) * dx_lev[0]*0.5; - Real fac_y = (1.0 - mfy_type[1]) * dx_lev[1]*0.5; Real x = i*dx_lev[0] + real_box.lo(0) + fac_x; - Real y = j*dx_lev[1] + real_box.lo(1) + fac_y; #if (AMREX_SPACEDIM==2) - Real z = 0.0; + Real y = 0.0; + Real fac_z = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; + Real z = j*dx_lev[1] + real_box.lo(1) + fac_z; #elif (AMREX_SPACEDIM==3) - Real fac_z = (1.0 - mfy_type[2]) * dx_lev[2]*0.5; + Real fac_y = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; + Real y = j*dx_lev[1] + real_box.lo(1) + fac_y; + Real fac_z = (1.0 - mfx_type[2]) * dx_lev[2]*0.5; Real z = k*dx_lev[2] + real_box.lo(2) + fac_z; #endif // Initialize the y-component of the field. @@ -629,8 +632,12 @@ WarpX::InitializeExternalFieldsOnGridUsingParser ( Real x = i*dx_lev[0] + real_box.lo(0) + fac_x; Real y = j*dx_lev[1] + real_box.lo(1) + fac_y; #if (AMREX_SPACEDIM==2) - Real z = 0.0; + Real y = 0.0; + Real fac_z = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; + Real z = j*dx_lev[1] + real_box.lo(1) + fac_z; #elif (AMREX_SPACEDIM==3) + Real fac_y = (1.0 - mfx_type[1]) * dx_lev[1]*0.5; + Real y = j*dx_lev[1] + real_box.lo(1) + fac_y; Real fac_z = (1.0 - mfz_type[2]) * dx_lev[2]*0.5; Real z = k*dx_lev[2] + real_box.lo(2) + fac_z; #endif |