2 files changed, 9 insertions, 4 deletions
diff --git a/Source/Initialization/PlasmaInjector.H b/Source/Initialization/PlasmaInjector.H
index ed87a7982..f6b76e8b6 100644
--- a/Source/Initialization/PlasmaInjector.H
+++ b/Source/Initialization/PlasmaInjector.H
@@ -274,19 +274,21 @@ public:
     amrex::Real q_tot;
     long npart;
 
+    bool radially_weighted = true;
+
     std::string str_density_function;
     std::string str_momentum_function_ux;
     std::string str_momentum_function_uy;
     std::string str_momentum_function_uz;
 
-protected:
-
-    amrex::Real mass, charge;
-
     amrex::Real xmin, xmax;
     amrex::Real ymin, ymax;
     amrex::Real zmin, zmax;
 
+protected:
+
+    amrex::Real mass, charge;
+
     amrex::Real density;
 
     int species_id;
diff --git a/Source/Initialization/PlasmaInjector.cpp b/Source/Initialization/PlasmaInjector.cpp
index 3d3edf715..0da9318de 100644
--- a/Source/Initialization/PlasmaInjector.cpp
+++ b/Source/Initialization/PlasmaInjector.cpp
@@ -299,6 +299,9 @@ PlasmaInjector::PlasmaInjector(int ispecies, const std::string& name)
         StringParseAbortMessage("Injection style", part_pos_s);
     }
 
+    pp.query("radially_weighted", radially_weighted);
+    AMREX_ALWAYS_ASSERT_WITH_MESSAGE(radially_weighted, "ERROR: Only radially_weighted=true is supported");
+
     // parse plasma boundaries
     xmin = std::numeric_limits<amrex::Real>::lowest();
     ymin = std::numeric_limits<amrex::Real>::lowest();