diff --git a/include/Realm.h b/include/Realm.h
index 343e5d5757..047253238f 100644
--- a/include/Realm.h
+++ b/include/Realm.h
@@ -183,7 +183,7 @@ class Realm
   void compute_vrtm(const std::string& = "velocity");
   void compute_l2_scaling();
   void output_converged_results();
-  void provide_output();
+  void provide_output(bool forcedOutput = false);
   void provide_restart_output();
 
   void register_interior_algorithm(stk::mesh::Part* part);
@@ -481,6 +481,7 @@ class Realm
 
   // check if there are negative Jacobians
   bool checkJacobians_;
+  bool outputFailedJacobians_;
 
   // types of physics
   bool isothermalFlow_;
diff --git a/src/Realm.C b/src/Realm.C
index 69743afcb3..0cd57013a3 100644
--- a/src/Realm.C
+++ b/src/Realm.C
@@ -240,6 +240,7 @@ Realm::Realm(Realms& realms, const YAML::Node& node)
     edgesPart_(0),
     checkForMissingBcs_(false),
     checkJacobians_(false),
+    outputFailedJacobians_(false),
     isothermalFlow_(true),
     uniformFlow_(true),
     provideEntityCount_(false),
@@ -670,6 +671,9 @@ Realm::load(const YAML::Node& node)
 
   // check for bad Jacobians in the mesh
   get_if_present(node, "check_jacobians", checkJacobians_, checkJacobians_);
+  get_if_present(
+    node, "output_on_failed_jacobian_check", outputFailedJacobians_,
+    outputFailedJacobians_);
 
   // entity count
   get_if_present(
@@ -3152,7 +3156,7 @@ Realm::overset_field_update(
 //-------- provide_output --------------------------------------------------
 //--------------------------------------------------------------------------
 void
-Realm::provide_output()
+Realm::provide_output(bool forcedOutput)
 {
   stk::diag::TimeBlock mesh_output_timeblock(Simulation::outputTimer());
   const double start_time = NaluEnv::self().nalu_time();
@@ -3169,7 +3173,6 @@ Realm::provide_output()
     const int modStep = timeStepCount - outputInfo_->outputStart_;
 
     // check for elapsed WALL time threshold
-    bool forcedOutput = false;
     if (outputInfo_->userWallTimeResults_.first) {
       const double elapsedWallTime = stk::wall_time() - wallTimeStart_;
       // find the max over all core
@@ -3179,8 +3182,9 @@ Realm::provide_output()
         1);
       // convert to hours
       g_elapsedWallTime /= 3600.0;
-      // only force output the first time the timer is exceeded
-      if (g_elapsedWallTime > outputInfo_->userWallTimeResults_.second) {
+      if (
+        g_elapsedWallTime > outputInfo_->userWallTimeResults_.second ||
+        forcedOutput) {
         forcedOutput = true;
         outputInfo_->userWallTimeResults_.first = false;
         NaluEnv::self().naluOutputP0() << "Realm::provide_output()::Forced "
diff --git a/src/ngp_algorithms/GeometryInteriorAlg.C b/src/ngp_algorithms/GeometryInteriorAlg.C
index 33b210648b..4e84cfaf8b 100644
--- a/src/ngp_algorithms/GeometryInteriorAlg.C
+++ b/src/ngp_algorithms/GeometryInteriorAlg.C
@@ -20,6 +20,7 @@
 #include "SolutionOptions.h"
 #include "utils/StkHelpers.h"
 #include "stk_mesh/base/NgpMesh.hpp"
+#include <stk_util/parallel/ParallelReduce.hpp>
 
 namespace sierra {
 namespace nalu {
@@ -57,8 +58,9 @@ template <typename AlgTraits>
 void
 GeometryInteriorAlg<AlgTraits>::execute()
 {
-  if (realm_.checkJacobians_)
+  if (realm_.checkJacobians_) {
     impl_negative_jacobian_check();
+  }
 
   impl_compute_dual_nodal_volume();
 
@@ -145,7 +147,12 @@ GeometryInteriorAlg<AlgTraits>::impl_negative_jacobian_check()
     },
     reducer);
 
-  if (numNegVol > 0) {
+  size_t globalNegVol = 0;
+  stk::all_reduce_sum(
+    NaluEnv::self().parallel_comm(), &numNegVol, &globalNegVol, 1);
+
+  if (globalNegVol > 0) {
+    realm_.provide_output(realm_.outputFailedJacobians_);
     const stk::topology topology(AlgTraits::topo_);
     throw std::runtime_error(
       "GeometryInteriorAlg encountered " + std::to_string(numNegVol) +