diff options
author | 2018-12-11 10:37:21 -0800 | |
---|---|---|
committer | 2018-12-11 10:37:21 -0800 | |
commit | b53d4a7ea5c3963e16727aafdbda01771af04c0d (patch) | |
tree | 33ddccf2ceb4ae6670baa4ed31b9b97cb7035473 /Examples/Tests/PML/analysis_pml_ckc.py | |
parent | 4e4a84aba36b2d2168cf41077f596549180f0dd9 (diff) | |
download | WarpX-b53d4a7ea5c3963e16727aafdbda01771af04c0d.tar.gz WarpX-b53d4a7ea5c3963e16727aafdbda01771af04c0d.tar.zst WarpX-b53d4a7ea5c3963e16727aafdbda01771af04c0d.zip |
fix conflicts for merge revert
Diffstat (limited to 'Examples/Tests/PML/analysis_pml_ckc.py')
-rwxr-xr-x | Examples/Tests/PML/analysis_pml_ckc.py | 34 |
1 files changed, 0 insertions, 34 deletions
diff --git a/Examples/Tests/PML/analysis_pml_ckc.py b/Examples/Tests/PML/analysis_pml_ckc.py deleted file mode 100755 index 90d13cecc..000000000 --- a/Examples/Tests/PML/analysis_pml_ckc.py +++ /dev/null @@ -1,34 +0,0 @@ -#! /usr/bin/env python - -import sys -import yt ; yt.funcs.mylog.setLevel(0) -import numpy as np -import scipy.constants as scc - -filename = sys.argv[1] - -############################ -### INITIAL LASER ENERGY ### -############################ -energy_start = 9.1301289517e-08 - -########################## -### FINAL LASER ENERGY ### -########################## -ds = yt.load( filename ) -all_data_level_0 = ds.covering_grid(level=0,left_edge=ds.domain_left_edge, dims=ds.domain_dimensions) -Bx = all_data_level_0['boxlib', 'Bx'].v.squeeze() -By = all_data_level_0['boxlib', 'By'].v.squeeze() -Bz = all_data_level_0['boxlib', 'Bz'].v.squeeze() -Ex = all_data_level_0['boxlib', 'Ex'].v.squeeze() -Ey = all_data_level_0['boxlib', 'Ey'].v.squeeze() -Ez = all_data_level_0['boxlib', 'Ez'].v.squeeze() -energyE = np.sum(scc.epsilon_0/2*(Ex**2+Ey**2+Ez**2)) -energyB = np.sum(1./scc.mu_0/2*(Bx**2+By**2+Bz**2)) -energy_end = energyE + energyB - -Reflectivity = energy_end/energy_start -Reflectivity_theory = 1.8015e-06 - -assert( abs(Reflectivity-Reflectivity_theory) < 5./100 * Reflectivity_theory ) - |