diff options
Diffstat (limited to 'Docs/source')
| -rw-r--r-- | Docs/source/running_cpp/parameters.rst | 7 | ||||
| -rw-r--r-- | Docs/source/visualization/yt.rst | 19 |
2 files changed, 26 insertions, 0 deletions
diff --git a/Docs/source/running_cpp/parameters.rst b/Docs/source/running_cpp/parameters.rst index 139811ffc..592ea1b0b 100644 --- a/Docs/source/running_cpp/parameters.rst +++ b/Docs/source/running_cpp/parameters.rst @@ -902,6 +902,13 @@ Diagnostics and output The time interval between the back-transformed reduced diagnostics (where this time interval is expressed in the laboratory frame). +* ``slice.particle_slice_width_lab`` (`float`, in meters) + Only used when ``warpx.do_boosted_frame_diagnostic`` is ``1`` and + ``slice.num_slice_snapshots_lab`` is non-zero. Particles are + copied from the full back-transformed diagnostic to the reduced + slice diagnostic if there are within the user-defined width from + the slice region defined by ``slice.dom_lo`` and ``slice.dom_hi``. + Checkpoints and restart ----------------------- WarpX supports checkpoints/restart via AMReX. diff --git a/Docs/source/visualization/yt.rst b/Docs/source/visualization/yt.rst index e2f8b377f..7c98e05aa 100644 --- a/Docs/source/visualization/yt.rst +++ b/Docs/source/visualization/yt.rst @@ -168,6 +168,25 @@ If the back-transformed diagnostics are written in the HDF5 format (This can be plt.figure() plt.plot(f2['Ez'][nx2//2,0,:]) +The back-transformed particle data on the full and reduced diagnostic can be visualized as follows + +:: + species='ions' + iteration = 1 + + snapshot = './lab_frame_data/snapshots/' + 'snapshot' + str(iteration).zfill(5) + xbo = get_particle_field(snapshot, species, 'x') # Read particle data + ybo = get_particle_field(snapshot, species, 'y') + zbo = get_particle_field(snapshot, species, 'z') + + snapshot = './lab_frame_data/slices/' + 'slice' + str(iteration).zfill(5) + xbo_slice = get_particle_field(snapshot, species, 'x') # Read particle data + ybo_slice = get_particle_field(snapshot, species, 'y') + zbo_slice = get_particle_field(snapshot, species, 'z') + plt.figure() + plt.plot(xbo, ybo, 'r.', markersize=1.) + plt.plot(xbo_slice, ybo_slice, 'bx', markersize=1.) + Further information ------------------- |