diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index bfeb0597fcee8a370ffc0921f5bb0fead2ca9896..15cc2092be5affcb5cf4d4e4054a4b7954fc2f7b 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -364,6 +364,10 @@ int main (int argc, char *argv[]) try
         FieldVector<double,dim> resultantForce, resultantTorque;
         resultantForce  = rodAssembler.getResultantForce(couplingBoundary, rodX, resultantTorque);
 
+        // Flip orientation
+        resultantForce  *= -1;
+        resultantTorque *= -1;
+        
         std::cout << "resultant force: " << resultantForce << std::endl;
         std::cout << "resultant torque: " << resultantTorque << std::endl;