From d59fa46d24417b67554132bc666e45886160bd09 Mon Sep 17 00:00:00 2001
From: Weiqun Zhang <weiqunzhang@lbl.gov>
Date: Wed, 24 Jul 2019 19:43:18 -0700
Subject: Reimplement AddPlasma. Commits related to AddPlasma in hackathon
 branch are squashed into one.

---
 Source/Initialization/PlasmaProfiles.cpp | 41 --------------------------------
 1 file changed, 41 deletions(-)
 delete mode 100644 Source/Initialization/PlasmaProfiles.cpp

(limited to 'Source/Initialization/PlasmaProfiles.cpp')

diff --git a/Source/Initialization/PlasmaProfiles.cpp b/Source/Initialization/PlasmaProfiles.cpp
deleted file mode 100644
index d9d207f7e..000000000
--- a/Source/Initialization/PlasmaProfiles.cpp
+++ /dev/null
@@ -1,41 +0,0 @@
-#include <PlasmaInjector.H>
-#include <cmath>
-#include <iostream>
-#include <WarpXConst.H>
-
-using namespace amrex;
-
-Real PredefinedDensityProfile::getDensity(Real x, Real y, Real z) const {
-    Real n;
-    if ( which_profile == predefined_profile_flag::parabolic_channel ) {
-        n = ParabolicChannel(x,y,z);
-    }
-    return n;
-}
-
-///
-/// plateau between linear upramp and downramp, and parab transverse profile
-///
-Real PredefinedDensityProfile::ParabolicChannel(Real x, Real y, Real z) const {
-    //  params = [z_start   ramp_up   plateau   ramp_down   rc       n0]
-    Real z_start   = params[0];
-    Real ramp_up   = params[1];
-    Real plateau   = params[2];
-    Real ramp_down = params[3];
-    Real rc        = params[4];
-    Real n0        = params[5];
-    Real n;
-    Real kp = PhysConst::q_e/PhysConst::c*sqrt( n0/(PhysConst::m_e*PhysConst::ep0) );
-
-    if        ((z-z_start)>=0               and (z-z_start)<ramp_up                  ) {
-        n = (z-z_start)/ramp_up;
-    } else if ((z-z_start)>=ramp_up         and (z-z_start)<ramp_up+plateau          ) {
-        n = 1;
-    } else if ((z-z_start)>=ramp_up+plateau and (z-z_start)<ramp_up+plateau+ramp_down) {
-        n = 1-((z-z_start)-ramp_up-plateau)/ramp_down;
-    } else {
-        n = 0;
-    }
-    n *= n0*(1+4*(x*x+y*y)/(kp*kp*std::pow(rc,4)));
-    return n;
-}
-- 
cgit v1.2.3