diff options
| author |  MaxThevenet <mthevenet@lbl.gov>
| 2020-02-29 09:15:57 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2020-02-29 09:15:57 -0800 |
| commit | 04c5d80a62d82ed29c09bfd2fbed0522ad613f72 (patch) | |
| tree | 390c71fba1fc44ca45c99ec60d589a53c68729d4 /Examples/Modules | |
| parent | ce895b10fd90e584ea62256a2f4834bdce2715d9 (diff) | |
| download | WarpX-04c5d80a62d82ed29c09bfd2fbed0522ad613f72.tar.gz WarpX-04c5d80a62d82ed29c09bfd2fbed0522ad613f72.tar.zst WarpX-04c5d80a62d82ed29c09bfd2fbed0522ad613f72.zip | |
replace tau by optical_depth_QSR and optical_depth_BW (#742)
* replace tau by optical_depth_QSR and optical_depth_BW
* update python analysis scripts accordingly
* fix bug for BW
Diffstat (limited to 'Examples/Modules')
3 files changed, 4 insertions, 4 deletions
diff --git a/Examples/Modules/qed/breit_wheeler/analysis_2d_tau_init.py b/Examples/Modules/qed/breit_wheeler/analysis_2d_tau_init.py index ce21aa059..b5e470e61 100755 --- a/Examples/Modules/qed/breit_wheeler/analysis_2d_tau_init.py +++ b/Examples/Modules/qed/breit_wheeler/analysis_2d_tau_init.py @@ -22,7 +22,7 @@ def check(): data_set = yt.load(filename) all_data = data_set.all_data() - res_tau = all_data["photons", 'particle_tau'] + res_tau = all_data["photons", 'particle_optical_depth_BW'] loc, scale = st.expon.fit(res_tau) diff --git a/Examples/Modules/qed/breit_wheeler/analysis_3d_optical_depth_evolution.py b/Examples/Modules/qed/breit_wheeler/analysis_3d_optical_depth_evolution.py index 4ed0e6953..e5540ffa8 100755 --- a/Examples/Modules/qed/breit_wheeler/analysis_3d_optical_depth_evolution.py +++ b/Examples/Modules/qed/breit_wheeler/analysis_3d_optical_depth_evolution.py @@ -104,7 +104,7 @@ def check(): all_data_end = data_set_end.all_data() tau_end_avg = np.array([ - np.average(all_data_end[name, 'particle_tau']) + np.average(all_data_end[name, 'particle_optical_depth_BW']) for name in spec_names]) dNBW_dt_sim = (tau_begin_avg - tau_end_avg)/sim_time diff --git a/Examples/Modules/qed/quantum_synchrotron/analysis_2d_tau_init.py b/Examples/Modules/qed/quantum_synchrotron/analysis_2d_tau_init.py index 2cceed2e4..bdbcf78f3 100755 --- a/Examples/Modules/qed/quantum_synchrotron/analysis_2d_tau_init.py +++ b/Examples/Modules/qed/quantum_synchrotron/analysis_2d_tau_init.py @@ -23,8 +23,8 @@ def check(): data_set = yt.load(filename) all_data = data_set.all_data() - res_ele_tau = all_data["electrons", 'particle_tau'] - res_pos_tau = all_data["positrons", 'particle_tau'] + res_ele_tau = all_data["electrons", 'particle_optical_depth_QSR'] + res_pos_tau = all_data["positrons", 'particle_optical_depth_QSR'] loc_ele, scale_ele = st.expon.fit(res_ele_tau) loc_pos, scale_pos = st.expon.fit(res_pos_tau) |