Plotfile capability for BTD (#1576)

* add geometry object and plotfile variables for BTD

* temporarily disable species output

* fix indices for output multifab

* add diag geom for output

* rename as m_output_species

* delete unnecessary deletion of m_species

* fixing buffer and snapshot definitions and emptying/refilling of buffers

* Adding Plotfile capability for BTD

* added BTD plotfile impl in cmake

* fix eol

* add comments and remove commented line in OpnPMD

* remove unnecesary print statement

* eol

* remove unusued var

* Read Fab on Disk string from buffer header

* virtual function

* initializing snapshot geom

* read header without bcast within io scope and modify fabname in first fabheader

* remove tab

* add empty line back

1 files changed, 5 insertions, 4 deletions

diff --git a/Source/Diagnostics/FlushFormats/FlushFormatPlotfile.cpp b/Source/Diagnostics/FlushFormats/FlushFormatPlotfile.cpp
index cf7aa32a2..1462cfe29 100644
--- a/Source/Diagnostics/FlushFormats/FlushFormatPlotfile.cpp
+++ b/Source/Diagnostics/FlushFormats/FlushFormatPlotfile.cpp
@@ -49,7 +49,7 @@ FlushFormatPlotfile::WriteToFile (
 
     WriteJobInfo(filename);
 
-    WriteWarpXHeader(filename, particle_diags);
+    WriteWarpXHeader(filename, particle_diags, geom);
 
     VisMF::SetHeaderVersion(current_version);
 }
@@ -179,7 +179,8 @@ FlushFormatPlotfile::WriteJobInfo(const std::string& dir) const
 void
 FlushFormatPlotfile::WriteWarpXHeader(
     const std::string& name,
-    const amrex::Vector<ParticleDiag>& particle_diags) const
+    const amrex::Vector<ParticleDiag>& particle_diags,
+    amrex::Vector<amrex::Geometry>& geom) const
 {
     auto & warpx = WarpX::GetInstance();
     if (ParallelDescriptor::IOProcessor())
@@ -232,11 +233,11 @@ FlushFormatPlotfile::WriteWarpXHeader(
 
         // Geometry
         for (int i = 0; i < AMREX_SPACEDIM; ++i) {
-            HeaderFile << warpx.Geom(0).ProbLo(i) << ' ';
+            HeaderFile << geom[0].ProbLo(i) << ' ';
         }
         HeaderFile << '\n';
         for (int i = 0; i < AMREX_SPACEDIM; ++i) {
-            HeaderFile << warpx.Geom(0).ProbHi(i) << ' ';
+            HeaderFile << geom[0].ProbHi(i) << ' ';
         }
         HeaderFile << '\n';