Side_discontinuous writing option for Exodus

Refs idaholab#20191

framework/include/outputs/Exodus.h

@@ -192,6 +192,9 @@ class Exodus : public OversampleOutput
   /// Flag to output discontinuous format in Exodus
   bool _discontinuous;
 
+  /// Flag to output added disjoint fictitious sides for side_discontinuous variables
+  bool _side_discontinuous;
+
   /// Flag to output HDF5 format (when available) in Exodus
   bool _write_hdf5;
 };

framework/src/outputs/Exodus.C

@@ -70,6 +70,10 @@ Exodus::validParams()
   params.addParam<bool>(
       "discontinuous", false, "Enables discontinuous output format for Exodus files.");
 
+  // Flag for outputting added side elements (for side-discontinuous data) to Exodus
+  params.addParam<bool>(
+      "side_discontinuous", false, "Enables adding side-discontinuous output in Exodus files.");
+
   // Flag for outputting Exodus data in HDF5 format (when libMesh is
   // configured with HDF5 support).  libMesh wants to do so by default
   // (for backwards compatibility with libMesh HDF5 users), but we
@@ -98,6 +102,7 @@ Exodus::Exodus(const InputParameters & parameters)
     _overwrite(getParam<bool>("overwrite")),
     _output_dimension(getParam<MooseEnum>("output_dimension").getEnum<OutputDimension>()),
     _discontinuous(getParam<bool>("discontinuous")),
+    _side_discontinuous(getParam<bool>("side_discontinuous")),
     _write_hdf5(getParam<bool>("write_hdf5"))
 {
   if (isParamValid("use_problem_dimension"))
@@ -224,6 +229,9 @@ Exodus::outputSetup()
     _exodus_io_ptr->set_hdf5_writing(false);
   }
 
+  if (_side_discontinuous)
+    _exodus_io_ptr->write_added_sides(true);
+
   // Increment file number and set appending status, append if all the following conditions are met:
   //   (1) If the application is recovering (not restarting)
   //   (2) The mesh has NOT changed