diff --git a/src/SIM/SIMbase.C b/src/SIM/SIMbase.C
index 67281e7f..f960c2e2 100644
--- a/src/SIM/SIMbase.C
+++ b/src/SIM/SIMbase.C
@@ -1149,7 +1149,7 @@ bool SIMbase::extractLoadVec (Vector& loadVec, size_t idx, const char* hd) const
     return false;
 
 #if SP_DEBUG > 1
-  std::cout <<"\nLoad vector:"<< loadVec;
+  mySam->printVector(std::cout,loadVec,"\nLoad vector");
 #endif
   if (hd)
   {
diff --git a/src/SIM/SIMgeneric.C b/src/SIM/SIMgeneric.C
index cb66a8af..661fa5e2 100644
--- a/src/SIM/SIMgeneric.C
+++ b/src/SIM/SIMgeneric.C
@@ -69,11 +69,21 @@ Vector SIMgeneric::getInterfaceForces (const Vector& sf,
 
   IntVec glbNodes;
   this->getBoundaryNodes(code,glbNodes);
+#if SP_DEBUG > 1
+  std::cout <<"\nInternal nodal forces at interface "<< code
+            <<" with "<< glbNodes.size() <<" nodes"<< std::endl;
+#endif
 
   for (int inod : glbNodes)
   {
     double w = inod <= (int)weights.size() ? weights[inod-1] : 1.0;
     std::pair<int,int> dof = mySam->getNodeDOFs(inod);
+#if SP_DEBUG > 1
+    std::cout <<"Node "<< inod <<":";
+    for (int j = dof.first; j <= dof.second; j++)
+      std::cout <<" "<< sf(j);
+    std::cout << std::endl;
+#endif
     for (unsigned char i = 0; i < nsd; i++, dof.first++)
       if (dof.first <= dof.second && dof.first < (int)sf.size())
         force[i] += w*sf(dof.first);