aboutsummaryrefslogtreecommitdiff
path: root/Examples/Tests/PML/analysis_pml_ckc.py
diff options
context:
space:
mode:
authorGravatar MaxThevenet <mthevenet@lbl.gov> 2018-12-11 10:37:21 -0800
committerGravatar MaxThevenet <mthevenet@lbl.gov> 2018-12-11 10:37:21 -0800
commitb53d4a7ea5c3963e16727aafdbda01771af04c0d (patch)
tree33ddccf2ceb4ae6670baa4ed31b9b97cb7035473 /Examples/Tests/PML/analysis_pml_ckc.py
parent4e4a84aba36b2d2168cf41077f596549180f0dd9 (diff)
downloadWarpX-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-xExamples/Tests/PML/analysis_pml_ckc.py34
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 )
-