New Diagnostics: FieldMaximum (#988)

* Add MaxField reduced diagnostic

* Add test of the MaxField reduced diag

* Minor modifications

* Fix test input file

* Use read_raw_data.py module in test

* Merge FieldEnergy and MaxField reduced diags

* Add specific error message if old syntax is used

* Fix column header numbers

* Update copyright

* Change name of diag in input file from Field to FieldReduced

* Add comment when allocating MultiFAB

* Add BackwardCompatibility () function

* Update Source/Diagnostics/ReducedDiags/MultiReducedDiags.cpp

Co-authored-by: MaxThevenet <mthevenet@lbl.gov>

* Add FieldReduced.cpp in CMakeLists.txt

* Apply suggestions from code review

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Remove magic numbers and update MFIter construct

* Apply suggestions from code review

Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

* Fix typo

* Do a reduction on the fly instead of allocating a MultiFAB

* Interpolate fields to cell center before computing maximum value of |E| and |B|

* Reseparate FieldReduced into FieldEnergy and MaxField

* Fix few typos

* Replace MaxField with FieldMaximum

Co-authored-by: MaxThevenet <mthevenet@lbl.gov>
Co-authored-by: Axel Huebl <axel.huebl@plasma.ninja>

1 files changed, 6 insertions, 1 deletions

diff --git a/Source/Diagnostics/ReducedDiags/MultiReducedDiags.cpp b/Source/Diagnostics/ReducedDiags/MultiReducedDiags.cpp
index bc8a8fb9a..69ffad6ee 100644
--- a/Source/Diagnostics/ReducedDiags/MultiReducedDiags.cpp
+++ b/Source/Diagnostics/ReducedDiags/MultiReducedDiags.cpp
@@ -10,6 +10,7 @@
 #include "BeamRelevant.H"
 #include "ParticleEnergy.H"
 #include "FieldEnergy.H"
+#include "FieldMaximum.H"
 #include "MultiReducedDiags.H"
 
 #include <AMReX_ParmParse.H>
@@ -54,6 +55,11 @@ MultiReducedDiags::MultiReducedDiags ()
             m_multi_rd[i_rd].reset
                 ( new FieldEnergy(m_rd_names[i_rd]));
         }
+        else if (rd_type.compare("FieldMaximum") == 0)
+        {
+            m_multi_rd[i_rd].reset
+                ( new FieldMaximum(m_rd_names[i_rd]));
+        }
         else if (rd_type.compare("BeamRelevant") == 0)
         {
             m_multi_rd[i_rd].reset
@@ -101,7 +107,6 @@ void MultiReducedDiags::WriteToFile (int step)
     // loop over all reduced diags
     for (int i_rd = 0; i_rd < static_cast<int>(m_rd_names.size()); ++i_rd)
     {
-
         // Judge if the diags should be done
         if ( (step+1) % m_multi_rd[i_rd]->m_freq != 0 ) { continue; }