Renamed "particle_vel" to "particle_mom" (#3672)

* Renamed particle_vel to particle_mom

* Changed vel to u

* Updated gitignore with vscode-specific folder

* Updated docs for single particle

* Revert "Updated gitignore with vscode-specific folder"

This reverts commit cb44ef8d66c123daa93d6a13d5a1d88416e866bb.

* Fixed typo

* [pre-commit.ci] auto fixes from pre-commit.com hooks

for more information, see https://pre-commit.ci

* Fixed other typo

* Add species indication in warning

---------

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>

27 files changed, 84 insertions, 84 deletions

diff --git a/Examples/Tests/AcceleratorLattice/analysis.py b/Examples/Tests/AcceleratorLattice/analysis.py
index d2fd7f6ff..35da52443 100755
--- a/Examples/Tests/AcceleratorLattice/analysis.py
+++ b/Examples/Tests/AcceleratorLattice/analysis.py
@@ -70,7 +70,7 @@ read_lattice('lattice', z_location)
 
 # Fetch the initial position of the particle
 x0 = [float(x) for x in ds.parameters.get('electron.single_particle_pos').split()]
-ux0 = [float(x)*c for x in ds.parameters.get('electron.single_particle_vel').split()]
+ux0 = [float(x)*c for x in ds.parameters.get('electron.single_particle_u').split()]
 
 xx = x0[0]
 zz = x0[2]
diff --git a/Examples/Tests/AcceleratorLattice/inputs_quad_3d b/Examples/Tests/AcceleratorLattice/inputs_quad_3d
index 0a3b433f6..613ba7c9d 100644
--- a/Examples/Tests/AcceleratorLattice/inputs_quad_3d
+++ b/Examples/Tests/AcceleratorLattice/inputs_quad_3d
@@ -20,7 +20,7 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "SingleParticle"
 electron.single_particle_pos = 0.05 0.0 0.0
-electron.single_particle_vel = 0.0 0.0 0.1  # gamma*beta
+electron.single_particle_u = 0.0 0.0 0.1  # gamma*beta
 electron.single_particle_weight = 1.0
 
 lattice.elements = lattice1 lattice2
diff --git a/Examples/Tests/AcceleratorLattice/inputs_quad_boosted_3d b/Examples/Tests/AcceleratorLattice/inputs_quad_boosted_3d
index ac1ea3b92..668ec73d2 100644
--- a/Examples/Tests/AcceleratorLattice/inputs_quad_boosted_3d
+++ b/Examples/Tests/AcceleratorLattice/inputs_quad_boosted_3d
@@ -20,7 +20,7 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "SingleParticle"
 electron.single_particle_pos = 0.05 0.0 0.0
-electron.single_particle_vel = 0.0 0.0 2.0  # gamma*beta
+electron.single_particle_u = 0.0 0.0 2.0  # gamma*beta
 electron.single_particle_weight = 1.0
 
 lattice.elements = drift1 quad1 drift2 quad2
diff --git a/Examples/Tests/AcceleratorLattice/inputs_quad_moving_3d b/Examples/Tests/AcceleratorLattice/inputs_quad_moving_3d
index 98d3d40ad..e5ec0a595 100644
--- a/Examples/Tests/AcceleratorLattice/inputs_quad_moving_3d
+++ b/Examples/Tests/AcceleratorLattice/inputs_quad_moving_3d
@@ -24,7 +24,7 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "SingleParticle"
 electron.single_particle_pos = 0.05 0.0 0.0
-electron.single_particle_vel = 0.0 0.0 0.1  # gamma*beta
+electron.single_particle_u = 0.0 0.0 0.1  # gamma*beta
 electron.single_particle_weight = 1.0
 
 lattice.elements = drift1 quad1 drift2 quad2
diff --git a/Examples/Tests/boundaries/inputs_3d b/Examples/Tests/boundaries/inputs_3d
index 71399fa5e..e3223af25 100644
--- a/Examples/Tests/boundaries/inputs_3d
+++ b/Examples/Tests/boundaries/inputs_3d
@@ -37,9 +37,9 @@ reflecting_particles.injection_style = "MultipleParticles"
 reflecting_particles.multiple_particles_pos_x = -0.9 +0.91
 reflecting_particles.multiple_particles_pos_y = 0. 0.
 reflecting_particles.multiple_particles_pos_z = 0. 0.
-reflecting_particles.multiple_particles_vel_x = -0.9 +0.91
-reflecting_particles.multiple_particles_vel_y = 0. 0.
-reflecting_particles.multiple_particles_vel_z = 0. 0.
+reflecting_particles.multiple_particles_ux = -0.9 +0.91
+reflecting_particles.multiple_particles_uy = 0. 0.
+reflecting_particles.multiple_particles_uz = 0. 0.
 reflecting_particles.multiple_particles_weight = 1. 1.
 
 absorbing_particles.charge = 0.
@@ -48,9 +48,9 @@ absorbing_particles.injection_style = "MultipleParticles"
 absorbing_particles.multiple_particles_pos_x = 0. 0. 0.
 absorbing_particles.multiple_particles_pos_y = -0.92 +0.93 0.
 absorbing_particles.multiple_particles_pos_z = 0. 0. 0.
-absorbing_particles.multiple_particles_vel_x = 0. 0. 0.
-absorbing_particles.multiple_particles_vel_y = -0.92 +0.93 0.
-absorbing_particles.multiple_particles_vel_z = 0. 0. 0.
+absorbing_particles.multiple_particles_ux = 0. 0. 0.
+absorbing_particles.multiple_particles_uy = -0.92 +0.93 0.
+absorbing_particles.multiple_particles_uz = 0. 0. 0.
 absorbing_particles.multiple_particles_weight = 1. 1. 1.
 
 periodic_particles.charge = 0.
@@ -59,9 +59,9 @@ periodic_particles.injection_style = "MultipleParticles"
 periodic_particles.multiple_particles_pos_x = 0. 0.
 periodic_particles.multiple_particles_pos_y = 0. 0.
 periodic_particles.multiple_particles_pos_z = -0.94 +0.95
-periodic_particles.multiple_particles_vel_x = 0. 0.
-periodic_particles.multiple_particles_vel_y = 0. 0.
-periodic_particles.multiple_particles_vel_z = -0.94 +0.95
+periodic_particles.multiple_particles_ux = 0. 0.
+periodic_particles.multiple_particles_uy = 0. 0.
+periodic_particles.multiple_particles_uz = -0.94 +0.95
 periodic_particles.multiple_particles_weight = 1. 1.
 
 # Diagnostics
diff --git a/Examples/Tests/ion_stopping/inputs_3d b/Examples/Tests/ion_stopping/inputs_3d
index 9352946b8..291e1ca0a 100644
--- a/Examples/Tests/ion_stopping/inputs_3d
+++ b/Examples/Tests/ion_stopping/inputs_3d
@@ -22,9 +22,9 @@ ions1.injection_style = "MultipleParticles"
 ions1.multiple_particles_pos_x = 0. 0. 0. 0.
 ions1.multiple_particles_pos_y = 0. 0. 0. 0.
 ions1.multiple_particles_pos_z = 0. 0. 0. 0.
-ions1.multiple_particles_vel_x = max_ion_velocity12 0. 0. max_ion_velocity12/2.
-ions1.multiple_particles_vel_y = 0. max_ion_velocity12 0. max_ion_velocity12/2.
-ions1.multiple_particles_vel_z = 0. 0. -max_ion_velocity12 -max_ion_velocity12/2.
+ions1.multiple_particles_ux = max_ion_velocity12 0. 0. max_ion_velocity12/2.
+ions1.multiple_particles_uy = 0. max_ion_velocity12 0. max_ion_velocity12/2.
+ions1.multiple_particles_uz = 0. 0. -max_ion_velocity12 -max_ion_velocity12/2.
 ions1.multiple_particles_weight = 1. 1. 1. 1.
 ions1.do_not_deposit = 1
 
@@ -33,9 +33,9 @@ ions2.injection_style = "MultipleParticles"
 ions2.multiple_particles_pos_x = -0.1 +0.1 -0.1 +0.1
 ions2.multiple_particles_pos_y = -0.1 -0.1 +0.1 +0.1
 ions2.multiple_particles_pos_z = 0. 0. 0. 0.
-ions2.multiple_particles_vel_x = max_ion_velocity12 max_ion_velocity12 max_ion_velocity12 max_ion_velocity12
-ions2.multiple_particles_vel_y = 0. 0. 0. 0.
-ions2.multiple_particles_vel_z = 0. 0. 0. 0.
+ions2.multiple_particles_ux = max_ion_velocity12 max_ion_velocity12 max_ion_velocity12 max_ion_velocity12
+ions2.multiple_particles_uy = 0. 0. 0. 0.
+ions2.multiple_particles_uz = 0. 0. 0. 0.
 ions2.multiple_particles_weight = 1. 1. 1. 1.
 ions2.do_not_deposit = 1
 
@@ -45,9 +45,9 @@ ions3.injection_style = "MultipleParticles"
 ions3.multiple_particles_pos_x = 0. 0. 0. 0.
 ions3.multiple_particles_pos_y = 0. 0. 0. 0.
 ions3.multiple_particles_pos_z = 0. 0. 0. 0.
-ions3.multiple_particles_vel_x = max_ion_velocity34 0. 0. max_ion_velocity34/2.
-ions3.multiple_particles_vel_y = 0. max_ion_velocity34 0. max_ion_velocity34/2.
-ions3.multiple_particles_vel_z = 0. 0. -max_ion_velocity34 -max_ion_velocity34/2.
+ions3.multiple_particles_ux = max_ion_velocity34 0. 0. max_ion_velocity34/2.
+ions3.multiple_particles_uy = 0. max_ion_velocity34 0. max_ion_velocity34/2.
+ions3.multiple_particles_uz = 0. 0. -max_ion_velocity34 -max_ion_velocity34/2.
 ions3.multiple_particles_weight = 1. 1. 1. 1.
 ions3.do_not_deposit = 1
 
@@ -57,9 +57,9 @@ ions4.injection_style = "MultipleParticles"
 ions4.multiple_particles_pos_x = -0.1 +0.1 -0.1 +0.1
 ions4.multiple_particles_pos_y = -0.1 -0.1 +0.1 +0.1
 ions4.multiple_particles_pos_z = 0. 0. 0. 0.
-ions4.multiple_particles_vel_x = max_ion_velocity34 max_ion_velocity34 max_ion_velocity34 max_ion_velocity34
-ions4.multiple_particles_vel_y = 0. 0. 0. 0.
-ions4.multiple_particles_vel_z = 0. 0. 0. 0.
+ions4.multiple_particles_ux = max_ion_velocity34 max_ion_velocity34 max_ion_velocity34 max_ion_velocity34
+ions4.multiple_particles_uy = 0. 0. 0. 0.
+ions4.multiple_particles_uz = 0. 0. 0. 0.
 ions4.multiple_particles_weight = 1. 1. 1. 1.
 ions4.do_not_deposit = 1
 
diff --git a/Examples/Tests/larmor/inputs_2d_mr b/Examples/Tests/larmor/inputs_2d_mr
index baf10e085..5d7af3d67 100644
--- a/Examples/Tests/larmor/inputs_2d_mr
+++ b/Examples/Tests/larmor/inputs_2d_mr
@@ -48,13 +48,13 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "SingleParticle"
 electron.single_particle_pos = 0.0  0.0  -1.25
-electron.single_particle_vel = -0.45825756949558416  0.0  0.0   # gamma*beta
+electron.single_particle_u = -0.45825756949558416  0.0  0.0   # gamma*beta
 
 positron.charge = q_e
 positron.mass = m_e
 positron.injection_style = "SingleParticle"
 positron.single_particle_pos = 0.0  0.0  -1.25
-positron.single_particle_vel = 0.45825756949558416  0.0  0.0   # gamma*beta
+positron.single_particle_u = 0.45825756949558416  0.0  0.0   # gamma*beta
 
 electron.single_particle_weight = 1.0e12
 positron.single_particle_weight = 1.0e12
diff --git a/Examples/Tests/particle_pusher/inputs_3d b/Examples/Tests/particle_pusher/inputs_3d
index 8f336fe66..31ab6800d 100644
--- a/Examples/Tests/particle_pusher/inputs_3d
+++ b/Examples/Tests/particle_pusher/inputs_3d
@@ -33,7 +33,7 @@ positron.charge = 1.0
 positron.mass = 1.0
 positron.injection_style = "SingleParticle"
 positron.single_particle_pos = 0.0  0.0  0.0
-positron.single_particle_vel = 0.0  19.974984355438178  0.0
+positron.single_particle_u = 0.0  19.974984355438178  0.0
 positron.single_particle_weight = 0.0
 
 # External fields
diff --git a/Examples/Tests/particles_in_pml/inputs_2d b/Examples/Tests/particles_in_pml/inputs_2d
index 14fb4f423..7cddc12b4 100644
--- a/Examples/Tests/particles_in_pml/inputs_2d
+++ b/Examples/Tests/particles_in_pml/inputs_2d
@@ -37,14 +37,14 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "singleparticle"
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 2. 0. 0.
+electron.single_particle_u = 2. 0. 0.
 electron.single_particle_weight = 1.
 
 proton.charge = q_e
 proton.mass = m_p    # Very heavy ; should not move
 proton.injection_style = "singleparticle"
 proton.single_particle_pos = 0. 0. 0.
-proton.single_particle_vel = -2. 0. 0.
+proton.single_particle_u = -2. 0. 0.
 proton.single_particle_weight = 1.
 
 # Order of particle shape factors
diff --git a/Examples/Tests/particles_in_pml/inputs_3d b/Examples/Tests/particles_in_pml/inputs_3d
index 6b92ad630..ad9e61c29 100644
--- a/Examples/Tests/particles_in_pml/inputs_3d
+++ b/Examples/Tests/particles_in_pml/inputs_3d
@@ -37,14 +37,14 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "singleparticle"
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 2. 0. 0.
+electron.single_particle_u = 2. 0. 0.
 electron.single_particle_weight = 1.
 
 proton.charge = q_e
 proton.mass = m_p    # Very heavy ; should not move
 proton.injection_style = "singleparticle"
 proton.single_particle_pos = 0. 0. 0.
-proton.single_particle_vel = -2. 0. 0.
+proton.single_particle_u = -2. 0. 0.
 proton.single_particle_weight = 1.
 
 # Order of particle shape factors
diff --git a/Examples/Tests/particles_in_pml/inputs_mr_2d b/Examples/Tests/particles_in_pml/inputs_mr_2d
index d56a3dc95..a6db42a38 100644
--- a/Examples/Tests/particles_in_pml/inputs_mr_2d
+++ b/Examples/Tests/particles_in_pml/inputs_mr_2d
@@ -39,14 +39,14 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "singleparticle"
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 2. 0. 0.
+electron.single_particle_u = 2. 0. 0.
 electron.single_particle_weight = 1.
 
 proton.charge = q_e
 proton.mass = m_p    # Very heavy ; should not move
 proton.injection_style = "singleparticle"
 proton.single_particle_pos = 0. 0. 0.
-proton.single_particle_vel = -2. 0. 0.
+proton.single_particle_u = -2. 0. 0.
 proton.single_particle_weight = 1.
 
 # Order of particle shape factors
diff --git a/Examples/Tests/particles_in_pml/inputs_mr_3d b/Examples/Tests/particles_in_pml/inputs_mr_3d
index b4b8ceb2d..ef2c4418c 100644
--- a/Examples/Tests/particles_in_pml/inputs_mr_3d
+++ b/Examples/Tests/particles_in_pml/inputs_mr_3d
@@ -39,14 +39,14 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "singleparticle"
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 2. 0. 0.
+electron.single_particle_u = 2. 0. 0.
 electron.single_particle_weight = 1.
 
 proton.charge = q_e
 proton.mass = m_p    # Very heavy ; should not move
 proton.injection_style = "singleparticle"
 proton.single_particle_pos = 0. 0. 0.
-proton.single_particle_vel = -2. 0. 0.
+proton.single_particle_u = -2. 0. 0.
 proton.single_particle_weight = 1.
 
 # Order of particle shape factors
diff --git a/Examples/Tests/pec/inputs_particle_PEC_3d b/Examples/Tests/pec/inputs_particle_PEC_3d
index 278eb4d47..ded8582b0 100644
--- a/Examples/Tests/pec/inputs_particle_PEC_3d
+++ b/Examples/Tests/pec/inputs_particle_PEC_3d
@@ -37,14 +37,14 @@ electron.charge = -q_e
 electron.mass = m_p
 electron.injection_style = "singleparticle"
 electron.single_particle_pos = 31.998e-6 0. 0.
-electron.single_particle_vel = 0. 0. 0.
+electron.single_particle_u = 0. 0. 0.
 electron.single_particle_weight = 1.0
 
 proton.charge = q_e
 proton.mass = m_p    # Very heavy ; should not move
 proton.injection_style = "singleparticle"
 proton.single_particle_pos = 31.998e-6 0.e-6 0.
-proton.single_particle_vel = 0. -2. 0.
+proton.single_particle_u = 0. -2. 0.
 proton.single_particle_weight = 1.0
 
 # Particle shape factor in each direction
diff --git a/Examples/Tests/photon_pusher/analysis_photon_pusher.py b/Examples/Tests/photon_pusher/analysis_photon_pusher.py
index 92d1e95f7..7518bd5ad 100755
--- a/Examples/Tests/photon_pusher/analysis_photon_pusher.py
+++ b/Examples/Tests/photon_pusher/analysis_photon_pusher.py
@@ -149,7 +149,7 @@ def generate():
             f.write('{}.injection_style = "SingleParticle"\n'.format(name))
             f.write("{}.single_particle_pos = {} {} {}\n".
                 format(name, init_pos[0], init_pos[1], init_pos[2]))
-            f.write("{}.single_particle_vel = {} {} {}\n".
+            f.write("{}.single_particle_u = {} {} {}\n".
                 format(name, velx, vely, velz))
             f.write("{}.single_particle_weight = 1.0\n".format(name))
             f.write("\n".format(name))
diff --git a/Examples/Tests/photon_pusher/inputs_3d b/Examples/Tests/photon_pusher/inputs_3d
index a9ff8ac7b..e2c75981e 100644
--- a/Examples/Tests/photon_pusher/inputs_3d
+++ b/Examples/Tests/photon_pusher/inputs_3d
@@ -28,112 +28,112 @@ p_xp_1.charge = -q_e
 p_xp_1.mass = m_e
 p_xp_1.injection_style = "SingleParticle"
 p_xp_1.single_particle_pos = 0.0 0.0 0.0
-p_xp_1.single_particle_vel = 1.0 0.0 0.0
+p_xp_1.single_particle_u = 1.0 0.0 0.0
 p_xp_1.single_particle_weight = 1.0
 
 p_xn_1.charge = -q_e
 p_xn_1.mass = m_e
 p_xn_1.injection_style = "SingleParticle"
 p_xn_1.single_particle_pos = 0.0 0.0 0.0
-p_xn_1.single_particle_vel = -1.0 0.0 0.0
+p_xn_1.single_particle_u = -1.0 0.0 0.0
 p_xn_1.single_particle_weight = 1.0
 
 p_yp_1.charge = -q_e
 p_yp_1.mass = m_e
 p_yp_1.injection_style = "SingleParticle"
 p_yp_1.single_particle_pos = 0.0 0.0 0.0
-p_yp_1.single_particle_vel = 0.0 1.0 0.0
+p_yp_1.single_particle_u = 0.0 1.0 0.0
 p_yp_1.single_particle_weight = 1.0
 
 p_yn_1.charge = -q_e
 p_yn_1.mass = m_e
 p_yn_1.injection_style = "SingleParticle"
 p_yn_1.single_particle_pos = 0.0 0.0 0.0
-p_yn_1.single_particle_vel = 0.0 -1.0 0.0
+p_yn_1.single_particle_u = 0.0 -1.0 0.0
 p_yn_1.single_particle_weight = 1.0
 
 p_zp_1.charge = -q_e
 p_zp_1.mass = m_e
 p_zp_1.injection_style = "SingleParticle"
 p_zp_1.single_particle_pos = 0.0 0.0 0.0
-p_zp_1.single_particle_vel = 0.0 0.0 1.0
+p_zp_1.single_particle_u = 0.0 0.0 1.0
 p_zp_1.single_particle_weight = 1.0
 
 p_zn_1.charge = -q_e
 p_zn_1.mass = m_e
 p_zn_1.injection_style = "SingleParticle"
 p_zn_1.single_particle_pos = 0.0 0.0 0.0
-p_zn_1.single_particle_vel = 0.0 0.0 -1.0
+p_zn_1.single_particle_u = 0.0 0.0 -1.0
 p_zn_1.single_particle_weight = 1.0
 
 p_dp_1.charge = -q_e
 p_dp_1.mass = m_e
 p_dp_1.injection_style = "SingleParticle"
 p_dp_1.single_particle_pos = 0.0 0.0 0.0
-p_dp_1.single_particle_vel = 1.0 1.0 1.0
+p_dp_1.single_particle_u = 1.0 1.0 1.0
 p_dp_1.single_particle_weight = 1.0
 
 p_dn_1.charge = -q_e
 p_dn_1.mass = m_e
 p_dn_1.injection_style = "SingleParticle"
 p_dn_1.single_particle_pos = 0.0 0.0 0.0
-p_dn_1.single_particle_vel = -1.0 -1.0 -1.0
+p_dn_1.single_particle_u = -1.0 -1.0 -1.0
 p_dn_1.single_particle_weight = 1.0
 
 p_xp_10.charge = -q_e
 p_xp_10.mass = m_e
 p_xp_10.injection_style = "SingleParticle"
 p_xp_10.single_particle_pos = 0.0 0.0 0.0
-p_xp_10.single_particle_vel = 10.0 0.0 0.0
+p_xp_10.single_particle_u = 10.0 0.0 0.0
 p_xp_10.single_particle_weight = 1.0
 
 p_xn_10.charge = -q_e
 p_xn_10.mass = m_e
 p_xn_10.injection_style = "SingleParticle"
 p_xn_10.single_particle_pos = 0.0 0.0 0.0
-p_xn_10.single_particle_vel = -10.0 0.0 0.0
+p_xn_10.single_particle_u = -10.0 0.0 0.0
 p_xn_10.single_particle_weight = 1.0
 
 p_yp_10.charge = -q_e
 p_yp_10.mass = m_e
 p_yp_10.injection_style = "SingleParticle"
 p_yp_10.single_particle_pos = 0.0 0.0 0.0
-p_yp_10.single_particle_vel = 0.0 10.0 0.0
+p_yp_10.single_particle_u = 0.0 10.0 0.0
 p_yp_10.single_particle_weight = 1.0
 
 p_yn_10.charge = -q_e
 p_yn_10.mass = m_e
 p_yn_10.injection_style = "SingleParticle"
 p_yn_10.single_particle_pos = 0.0 0.0 0.0
-p_yn_10.single_particle_vel = 0.0 -10.0 0.0
+p_yn_10.single_particle_u = 0.0 -10.0 0.0
 p_yn_10.single_particle_weight = 1.0
 
 p_zp_10.charge = -q_e
 p_zp_10.mass = m_e
 p_zp_10.injection_style = "SingleParticle"
 p_zp_10.single_particle_pos = 0.0 0.0 0.0
-p_zp_10.single_particle_vel = 0.0 0.0 10.0
+p_zp_10.single_particle_u = 0.0 0.0 10.0
 p_zp_10.single_particle_weight = 1.0
 
 p_zn_10.charge = -q_e
 p_zn_10.mass = m_e
 p_zn_10.injection_style = "SingleParticle"
 p_zn_10.single_particle_pos = 0.0 0.0 0.0
-p_zn_10.single_particle_vel = 0.0 0.0 -10.0
+p_zn_10.single_particle_u = 0.0 0.0 -10.0
 p_zn_10.single_particle_weight = 1.0
 
 p_dp_10.charge = -q_e
 p_dp_10.mass = m_e
 p_dp_10.injection_style = "SingleParticle"
 p_dp_10.single_particle_pos = 0.0 0.0 0.0
-p_dp_10.single_particle_vel = 10.0 10.0 10.0
+p_dp_10.single_particle_u = 10.0 10.0 10.0
 p_dp_10.single_particle_weight = 1.0
 
 p_dn_10.charge = -q_e
 p_dn_10.mass = m_e
 p_dn_10.injection_style = "SingleParticle"
 p_dn_10.single_particle_pos = 0.0 0.0 0.0
-p_dn_10.single_particle_vel = -10.0 -10.0 -10.0
+p_dn_10.single_particle_u = -10.0 -10.0 -10.0
 p_dn_10.single_particle_weight = 1.0
 
 # Diagnostics
diff --git a/Examples/Tests/plasma_lens/analysis.py b/Examples/Tests/plasma_lens/analysis.py
index 80fdc4181..34478c248 100755
--- a/Examples/Tests/plasma_lens/analysis.py
+++ b/Examples/Tests/plasma_lens/analysis.py
@@ -95,9 +95,9 @@ elif 'lattice.elements' in ds.parameters:
 x0 = float(ds.parameters.get('electrons.multiple_particles_pos_x').split()[0])
 y0 = float(ds.parameters.get('electrons.multiple_particles_pos_y').split()[1])
 z0 = float(ds.parameters.get('electrons.multiple_particles_pos_z').split()[0])
-ux0 = float(ds.parameters.get('electrons.multiple_particles_vel_x').split()[0])*c
-uy0 = float(ds.parameters.get('electrons.multiple_particles_vel_y').split()[1])*c
-uz0 = eval(ds.parameters.get('electrons.multiple_particles_vel_z').split()[0])*c
+ux0 = float(ds.parameters.get('electrons.multiple_particles_ux').split()[0])*c
+uy0 = float(ds.parameters.get('electrons.multiple_particles_uy').split()[1])*c
+uz0 = eval(ds.parameters.get('electrons.multiple_particles_uz').split()[0])*c
 
 tt = 0.
 xx = x0
diff --git a/Examples/Tests/plasma_lens/inputs_3d b/Examples/Tests/plasma_lens/inputs_3d
index 48f2263f0..7763013ad 100644
--- a/Examples/Tests/plasma_lens/inputs_3d
+++ b/Examples/Tests/plasma_lens/inputs_3d
@@ -31,9 +31,9 @@ electrons.injection_style = "MultipleParticles"
 electrons.multiple_particles_pos_x = 0.05 0.
 electrons.multiple_particles_pos_y = 0. 0.04
 electrons.multiple_particles_pos_z = 0.05 0.05
-electrons.multiple_particles_vel_x = 0. 0.
-electrons.multiple_particles_vel_y = 0. 0.
-electrons.multiple_particles_vel_z = vel_z/clight vel_z/clight
+electrons.multiple_particles_ux = 0. 0.
+electrons.multiple_particles_uy = 0. 0.
+electrons.multiple_particles_uz = vel_z/clight vel_z/clight
 electrons.multiple_particles_weight = 1. 1.
 
 particles.E_ext_particle_init_style = repeated_plasma_lens
diff --git a/Examples/Tests/plasma_lens/inputs_boosted_3d b/Examples/Tests/plasma_lens/inputs_boosted_3d
index c997136d4..ecdf3b3ba 100644
--- a/Examples/Tests/plasma_lens/inputs_boosted_3d
+++ b/Examples/Tests/plasma_lens/inputs_boosted_3d
@@ -33,9 +33,9 @@ electrons.injection_style = "MultipleParticles"
 electrons.multiple_particles_pos_x = 0.05 0.
 electrons.multiple_particles_pos_y = 0. 0.04
 electrons.multiple_particles_pos_z = 0.05 0.05
-electrons.multiple_particles_vel_x = 0. 0.
-electrons.multiple_particles_vel_y = 0. 0.
-electrons.multiple_particles_vel_z = vel_z/clight vel_z/clight
+electrons.multiple_particles_ux = 0. 0.
+electrons.multiple_particles_uy = 0. 0.
+electrons.multiple_particles_uz = vel_z/clight vel_z/clight
 electrons.multiple_particles_weight = 1. 1.
 
 particles.E_ext_particle_init_style = repeated_plasma_lens
diff --git a/Examples/Tests/plasma_lens/inputs_lattice_3d b/Examples/Tests/plasma_lens/inputs_lattice_3d
index 3473f1420..a41ca8bae 100644
--- a/Examples/Tests/plasma_lens/inputs_lattice_3d
+++ b/Examples/Tests/plasma_lens/inputs_lattice_3d
@@ -31,9 +31,9 @@ electrons.injection_style = "MultipleParticles"
 electrons.multiple_particles_pos_x = 0.05 0.
 electrons.multiple_particles_pos_y = 0. 0.04
 electrons.multiple_particles_pos_z = 0.05 0.05
-electrons.multiple_particles_vel_x = 0. 0.
-electrons.multiple_particles_vel_y = 0. 0.
-electrons.multiple_particles_vel_z = vel_z/clight vel_z/clight
+electrons.multiple_particles_ux = 0. 0.
+electrons.multiple_particles_uy = 0. 0.
+electrons.multiple_particles_uz = vel_z/clight vel_z/clight
 electrons.multiple_particles_weight = 1. 1.
 
 lattice.elements = drift1 plasmalens1 drift2 plasmalens2 drift3 plasmalens3 drift4 plasmalens4
diff --git a/Examples/Tests/plasma_lens/inputs_short_3d b/Examples/Tests/plasma_lens/inputs_short_3d
index 4beb65ef9..8fa111df3 100644
--- a/Examples/Tests/plasma_lens/inputs_short_3d
+++ b/Examples/Tests/plasma_lens/inputs_short_3d
@@ -35,9 +35,9 @@ electrons.injection_style = "MultipleParticles"
 electrons.multiple_particles_pos_x = 0.05 0.
 electrons.multiple_particles_pos_y = 0. 0.04
 electrons.multiple_particles_pos_z = 0.05 0.05
-electrons.multiple_particles_vel_x = 0. 0.
-electrons.multiple_particles_vel_y = 0. 0.
-electrons.multiple_particles_vel_z = vel_z/clight vel_z/clight
+electrons.multiple_particles_ux = 0. 0.
+electrons.multiple_particles_uy = 0. 0.
+electrons.multiple_particles_uz = vel_z/clight vel_z/clight
 electrons.multiple_particles_weight = 1. 1.
 
 particles.E_ext_particle_init_style = repeated_plasma_lens
diff --git a/Examples/Tests/pml/inputs_rz b/Examples/Tests/pml/inputs_rz
index b05db804a..9ca7e4894 100644
--- a/Examples/Tests/pml/inputs_rz
+++ b/Examples/Tests/pml/inputs_rz
@@ -37,7 +37,7 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "singleparticle"
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 10. 0. 0.
+electron.single_particle_u = 10. 0. 0.
 electron.single_particle_weight = 1.
 
 #################################
diff --git a/Examples/Tests/radiation_reaction/test_const_B_analytical/analysis_classicalRR.py b/Examples/Tests/radiation_reaction/test_const_B_analytical/analysis_classicalRR.py
index f242f970b..e123562fb 100755
--- a/Examples/Tests/radiation_reaction/test_const_B_analytical/analysis_classicalRR.py
+++ b/Examples/Tests/radiation_reaction/test_const_B_analytical/analysis_classicalRR.py
@@ -190,7 +190,7 @@ def generate():
             f.write('{}.injection_style = "SingleParticle"\n'.format(cc.name))
             f.write("{}.single_particle_pos = {} {} {}\n".
                 format(cc.name, init_pos[0], init_pos[1], init_pos[2]))
-            f.write("{}.single_particle_vel = {} {} {}\n".
+            f.write("{}.single_particle_u = {} {} {}\n".
                 format(cc.name, cc.init_mom[0], cc.init_mom[1], cc.init_mom[2]))
             f.write("{}.single_particle_weight = {}\n".format(cc.name, very_small_weight))
             f.write("{}.do_classical_radiation_reaction = 1\n".format(cc.name))
diff --git a/Examples/Tests/radiation_reaction/test_const_B_analytical/inputs_3d b/Examples/Tests/radiation_reaction/test_const_B_analytical/inputs_3d
index 3df828fd2..d18145168 100644
--- a/Examples/Tests/radiation_reaction/test_const_B_analytical/inputs_3d
+++ b/Examples/Tests/radiation_reaction/test_const_B_analytical/inputs_3d
@@ -27,35 +27,35 @@ particles.species_names =  ele_para0 ele_perp0 ele_perp1 ele_perp2 pos_perp2
 ele_para0.species_type = "electron"
 ele_para0.injection_style = "SingleParticle"
 ele_para0.single_particle_pos = 0.0 0.0 0.0
-ele_para0.single_particle_vel = 285.7142857142856 428.5714285714285 857.142857142857
+ele_para0.single_particle_u = 285.7142857142856 428.5714285714285 857.142857142857
 ele_para0.single_particle_weight = 1e-08
 ele_para0.do_classical_radiation_reaction = 1
 
 ele_perp0.species_type = "electron"
 ele_perp0.injection_style = "SingleParticle"
 ele_perp0.single_particle_pos = 0.0 0.0 0.0
-ele_perp0.single_particle_vel = -47.91574237499548 6.388765649999412 12.777531299998806
+ele_perp0.single_particle_u = -47.91574237499548 6.388765649999412 12.777531299998806
 ele_perp0.single_particle_weight = 1e-08
 ele_perp0.do_classical_radiation_reaction = 1
 
 ele_perp1.species_type = "electron"
 ele_perp1.injection_style = "SingleParticle"
 ele_perp1.single_particle_pos = 0.0 0.0 0.0
-ele_perp1.single_particle_vel = 1.4274296030894868e-14 178.88543819998318 -89.44271909999159
+ele_perp1.single_particle_u = 1.4274296030894868e-14 178.88543819998318 -89.44271909999159
 ele_perp1.single_particle_weight = 1e-08
 ele_perp1.do_classical_radiation_reaction = 1
 
 ele_perp2.species_type = "electron"
 ele_perp2.injection_style = "SingleParticle"
 ele_perp2.single_particle_pos = 0.0 0.0 0.0
-ele_perp2.single_particle_vel = -958.3148474999095 127.77531299998823 255.55062599997612
+ele_perp2.single_particle_u = -958.3148474999095 127.77531299998823 255.55062599997612
 ele_perp2.single_particle_weight = 1e-08
 ele_perp2.do_classical_radiation_reaction = 1
 
 pos_perp2.species_type = "positron"
 pos_perp2.injection_style = "SingleParticle"
 pos_perp2.single_particle_pos = 0.0 0.0 0.0
-pos_perp2.single_particle_vel = -958.3148474999095 127.77531299998823 255.55062599997612
+pos_perp2.single_particle_u = -958.3148474999095 127.77531299998823 255.55062599997612
 pos_perp2.single_particle_weight = 1e-08
 pos_perp2.do_classical_radiation_reaction = 1
 
diff --git a/Examples/Tests/repelling_particles/inputs_2d b/Examples/Tests/repelling_particles/inputs_2d
index 77c3850bd..a4abd67c9 100644
--- a/Examples/Tests/repelling_particles/inputs_2d
+++ b/Examples/Tests/repelling_particles/inputs_2d
@@ -31,7 +31,7 @@ electron1.charge = q_e
 electron1.mass = m_e
 electron1.injection_style = "singleparticle"
 electron1.single_particle_pos = 1.e-6 0. 0.
-electron1.single_particle_vel = 0. 0. 0.
+electron1.single_particle_u = 0. 0. 0.
 electron1.single_particle_weight = 5.e12
 electron1.initialize_self_fields = 1
 
@@ -39,7 +39,7 @@ electron2.charge = q_e
 electron2.mass = m_e
 electron2.injection_style = "singleparticle"
 electron2.single_particle_pos = -1.e-6 0. 0.
-electron2.single_particle_vel = 0. 0. 0.
+electron2.single_particle_u = 0. 0. 0.
 electron2.single_particle_weight = 5.e12
 electron2.initialize_self_fields = 1
 
diff --git a/Examples/Tests/single_particle/inputs_2d b/Examples/Tests/single_particle/inputs_2d
index bee31c020..71ca2101c 100644
--- a/Examples/Tests/single_particle/inputs_2d
+++ b/Examples/Tests/single_particle/inputs_2d
@@ -23,7 +23,7 @@ electron.charge = -q_e
 electron.mass = m_e
 electron.injection_style = "SingleParticle"
 electron.single_particle_pos = 0.0   0.0 0.0
-electron.single_particle_vel = 1.e20 0.0 0.0  # gamma*beta
+electron.single_particle_u = 1.e20 0.0 0.0  # gamma*beta
 electron.single_particle_weight = 1.0
 
 # Diagnostics
diff --git a/Examples/Tests/vay_deposition/inputs_2d b/Examples/Tests/vay_deposition/inputs_2d
index a47e1f832..c4928831a 100644
--- a/Examples/Tests/vay_deposition/inputs_2d
+++ b/Examples/Tests/vay_deposition/inputs_2d
@@ -26,7 +26,7 @@ electron.charge = -q_e
 electron.injection_style = "SingleParticle"
 electron.mass = 9e-22
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 0.5 0. 0.5
+electron.single_particle_u = 0.5 0. 0.5
 electron.single_particle_weight = 1.
 
 # geometry
@@ -40,7 +40,7 @@ ion.charge = q_e
 ion.injection_style = "SingleParticle"
 ion.mass = 9e-22
 ion.single_particle_pos = 0. 0. 0.
-ion.single_particle_vel = -0.5 0. -0.5
+ion.single_particle_u = -0.5 0. -0.5
 ion.single_particle_weight = 1.
 
 # max_step
diff --git a/Examples/Tests/vay_deposition/inputs_3d b/Examples/Tests/vay_deposition/inputs_3d
index 663da28e4..a672049e1 100644
--- a/Examples/Tests/vay_deposition/inputs_3d
+++ b/Examples/Tests/vay_deposition/inputs_3d
@@ -26,7 +26,7 @@ electron.charge = -q_e
 electron.injection_style = "SingleParticle"
 electron.mass = 9e-22
 electron.single_particle_pos = 0. 0. 0.
-electron.single_particle_vel = 0.5 0.5 0.5
+electron.single_particle_u = 0.5 0.5 0.5
 electron.single_particle_weight = 1.
 
 # geometry
@@ -40,7 +40,7 @@ ion.charge = q_e
 ion.injection_style = "SingleParticle"
 ion.mass = 9e-22
 ion.single_particle_pos = 0. 0. 0.
-ion.single_particle_vel = -0.5 -0.5 -0.5
+ion.single_particle_u = -0.5 -0.5 -0.5
 ion.single_particle_weight = 1.
 
 # max_step