1 files changed, 12 insertions, 0 deletions
diff --git a/Exec/plasma_acceleration/CustomDensityProb.cpp b/Exec/plasma_acceleration/CustomDensityProb.cpp
new file mode 100644
index 000000000..3efcb13c5
--- /dev/null
+++ b/Exec/plasma_acceleration/CustomDensityProb.cpp
@@ -0,0 +1,12 @@
+#include <PlasmaInjector.H>
+
+#include <iostream>
+
+using namespace amrex;
+
+///
+/// This "custom" density profile just does constant
+///
+Real CustomDensityProfile::getDensity(Real x, Real y, Real z) const {
+  return params[0];
+}