diff --git a/src/ASM/IntegrandBase.C b/src/ASM/IntegrandBase.C
index 2a62cd12..41ff57c8 100644
--- a/src/ASM/IntegrandBase.C
+++ b/src/ASM/IntegrandBase.C
@@ -32,8 +32,8 @@
 LocalIntegral* IntegrandBase::getLocalIntegral (size_t nen, size_t,
                                                 bool neumann) const
 {
-  ElmMats* result = new ElmMats(!neumann && m_mode != SIM::RECOVERY);
-  result->rhsOnly = m_mode == SIM::RHS_ONLY || m_mode == SIM::RECOVERY;
+  ElmMats* result = new ElmMats(!neumann && m_mode < SIM::RECOVERY);
+  result->rhsOnly = m_mode >= SIM::RHS_ONLY;
   result->resize(neumann ? 0 : 1, 1);
   result->redim(npv*nen);
 
diff --git a/src/SIM/SIMenums.h b/src/SIM/SIMenums.h
index 31aa2def..9088cd28 100644
--- a/src/SIM/SIMenums.h
+++ b/src/SIM/SIMenums.h
@@ -37,6 +37,7 @@ namespace SIM //! Simulation scope
     STIFF_ONLY,
     MASS_ONLY,
     RHS_ONLY,
+    INT_FORCES,
     RECOVERY
   };