diff options
Diffstat (limited to 'Source/Diagnostics/ParticleIO.cpp')
| -rw-r--r-- | Source/Diagnostics/ParticleIO.cpp | 16 |
1 files changed, 11 insertions, 5 deletions
diff --git a/Source/Diagnostics/ParticleIO.cpp b/Source/Diagnostics/ParticleIO.cpp index a226bb590..a7b389d78 100644 --- a/Source/Diagnostics/ParticleIO.cpp +++ b/Source/Diagnostics/ParticleIO.cpp @@ -41,9 +41,10 @@ using namespace amrex; void RigidInjectedParticleContainer::ReadHeader (std::istream& is) { - is >> charge >> mass; - WarpX::GotoNextLine(is); + // Call parent class + WarpXParticleContainer::ReadHeader( is ); + // Read quantities that are specific to rigid-injected species int nlevs; is >> nlevs; WarpX::GotoNextLine(is); @@ -52,24 +53,29 @@ RigidInjectedParticleContainer::ReadHeader (std::istream& is) for (int i = 0; i < nlevs; ++i) { - int zinject_plane_tmp; + amrex::Real zinject_plane_tmp; is >> zinject_plane_tmp; zinject_plane_levels.push_back(zinject_plane_tmp); WarpX::GotoNextLine(is); } + is >> vzbeam_ave_boosted; + WarpX::GotoNextLine(is); } void RigidInjectedParticleContainer::WriteHeader (std::ostream& os) const { - // no need to write species_id - os << charge << " " << mass << "\n"; + // Call parent class + WarpXParticleContainer::WriteHeader( os ); + + // Write quantities that are specific to the rigid-injected species int nlevs = zinject_plane_levels.size(); os << nlevs << "\n"; for (int i = 0; i < nlevs; ++i) { os << zinject_plane_levels[i] << "\n"; } + os << vzbeam_ave_boosted << "\n"; } void |