diff options
Diffstat (limited to 'Examples/Modules/boosted_diags')
| -rwxr-xr-x | Examples/Modules/boosted_diags/analysis.py (renamed from Examples/Modules/boosted_diags/analysis_3Dbacktransformed_diag.py) | 24 | ||||
| -rw-r--r-- | Examples/Modules/boosted_diags/inputs_3d (renamed from Examples/Modules/boosted_diags/inputs_3d_slice) | 16 |
2 files changed, 0 insertions, 40 deletions
diff --git a/Examples/Modules/boosted_diags/analysis_3Dbacktransformed_diag.py b/Examples/Modules/boosted_diags/analysis.py index 6fa4e9c93..c6c089f98 100755 --- a/Examples/Modules/boosted_diags/analysis_3Dbacktransformed_diag.py +++ b/Examples/Modules/boosted_diags/analysis.py @@ -21,7 +21,6 @@ import sys import numpy as np import openpmd_api as io -import read_raw_data import yt yt.funcs.mylog.setLevel(0) @@ -35,22 +34,6 @@ filename = sys.argv[1] rtol = 1e-16 atol = 1e-16 -# Read data from legacy back-transformed diagnostics (entire domain) -snapshot = './lab_frame_data/snapshots/snapshot00003' -header = './lab_frame_data/snapshots/Header' -allrd, info = read_raw_data.read_lab_snapshot(snapshot, header) -Ez_legacy = allrd['Ez'] -print(f'Ez_legacy.shape = {Ez_legacy.shape}') -Ez_legacy_1D = np.squeeze(Ez_legacy[Ez_legacy.shape[0]//2,Ez_legacy.shape[1]//2,:]) - -# Read data from reduced back-transformed diagnostics (slice) -snapshot_slice = './lab_frame_data/slices/slice00003' -header_slice = './lab_frame_data/slices/Header' -allrd, info = read_raw_data.read_lab_snapshot(snapshot_slice, header_slice) -Ez_legacy_slice = allrd['Ez'] -print(f'Ez_legacy_slice.shape = {Ez_legacy_slice.shape}') -Ez_legacy_slice_1D = np.squeeze(Ez_legacy_slice[Ez_legacy_slice.shape[0]//2,1,:]) - # Read data from new back-transformed diagnostics (plotfile) ds_plotfile = yt.load(filename) data = ds_plotfile.covering_grid( @@ -69,12 +52,5 @@ series.flush() # Compare arrays to check consistency between new BTD formats (plotfile and openPMD) assert(np.allclose(Ez_plotfile, Ez_openpmd, rtol=rtol, atol=atol)) -# Check slicing -err = np.max(np.abs(Ez_legacy_slice_1D-Ez_legacy_1D)) / np.max(np.abs(Ez_legacy_1D)) -tol = 1e-16 -print(f'error = {err}') -print(f'tolerance = {tol}') -assert(err < tol) - test_name = os.path.split(os.getcwd())[1] checksumAPI.evaluate_checksum(test_name, filename) diff --git a/Examples/Modules/boosted_diags/inputs_3d_slice b/Examples/Modules/boosted_diags/inputs_3d index ff68b5909..ba98558be 100644 --- a/Examples/Modules/boosted_diags/inputs_3d_slice +++ b/Examples/Modules/boosted_diags/inputs_3d @@ -55,7 +55,6 @@ electrons.zmax = .003 electrons.profile = constant electrons.density = 3.5e24 electrons.do_continuous_injection = 1 -electrons.do_back_transformed_diagnostics = 1 ions.charge = q_e ions.mass = m_p @@ -71,7 +70,6 @@ ions.zmax = .003 ions.profile = constant ions.density = 3.5e24 ions.do_continuous_injection = 1 -ions.do_back_transformed_diagnostics = 1 beam.charge = -q_e beam.mass = m_e @@ -104,14 +102,6 @@ laser1.profile_t_peak = 40.e-15 # The time at which the laser reaches its pea laser1.profile_focal_distance = 0.5e-3 # Focal distance from the antenna (in meters) laser1.wavelength = 0.81e-6 # The wavelength of the laser (in meters) -slice.dom_lo = xmin 0.0 zmin -slice.dom_hi = xmax 0.0 zmax -slice.coarsening_ratio = 1 1 1 -slice.plot_int = -1 -slice.num_slice_snapshots_lab = 4 -slice.dt_slice_snapshots_lab = 3.3356409519815207e-12 -slice.particle_slice_width_lab = 2.e-6 - # Diagnostics diagnostics.diags_names = diag1 diag2 @@ -132,9 +122,3 @@ diag2.fields_to_plot = Ex Ey Ez Bx By Bz jx jy jz rho diag2.format = openpmd diag2.buffer_size = 32 diag2.openpmd_backend = h5 - -# old BTD diagnostics -warpx.do_back_transformed_diagnostics = 1 -warpx.num_snapshots_lab = 4 -warpx.dz_snapshots_lab = 0.001 -warpx.back_transformed_diag_fields= Ex Ey Ez By rho |