Merge pull request #417 from gtrichardson/QED_Hybrid_algorithm

QED Corrections to Maxwell's Equations (Hybrid)

1 files changed, 4 insertions, 0 deletions

diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp
index d93fab7df..3fe61274a 100644
--- a/Source/WarpX.cpp
+++ b/Source/WarpX.cpp
@@ -55,6 +55,7 @@ int WarpX::do_moving_window = 0;
 int WarpX::moving_window_dir = -1;
 Real WarpX::moving_window_v = std::numeric_limits<amrex::Real>::max();
 
+Real WarpX::quantum_xi = PhysConst::xi;
 Real WarpX::gamma_boost = 1.;
 Real WarpX::beta_boost = 0.;
 Vector<int> WarpX::boost_direction = {0,0,0};
@@ -341,6 +342,7 @@ WarpX::ReadParameters ()
         pp.query("verbose", verbose);
         pp.query("regrid_int", regrid_int);
         pp.query("do_subcycling", do_subcycling);
+        pp.query("use_hybrid_QED", use_hybrid_QED);
         pp.query("exchange_all_guard_cells", exchange_all_guard_cells);
         pp.query("override_sync_int", override_sync_int);
 
@@ -456,6 +458,8 @@ WarpX::ReadParameters ()
         pp.query("n_current_deposition_buffer", n_current_deposition_buffer);
         pp.query("sort_int", sort_int);
 
+        pp.query("quantum_xi", quantum_xi);
+
         pp.query("do_pml", do_pml);
         pp.query("pml_ncell", pml_ncell);
         pp.query("pml_delta", pml_delta);