diff --git a/src/MDI/fix_mdi_qm.cpp b/src/MDI/fix_mdi_qm.cpp
index 25e6870bd8..535359f069 100644
--- a/src/MDI/fix_mdi_qm.cpp
+++ b/src/MDI/fix_mdi_qm.cpp
@@ -348,6 +348,8 @@ void FixMDIQM::post_force(int vflag)
   if (ierr) error->all(FLERR, "MDI: <FORCES data");
   MPI_Bcast(buf3, 3 * n, MPI_DOUBLE, 0, world);
 
+  printf("MDI FORCE: %g %g %g: %g %g %g\n",buf3[0],buf3[1],buf3[2],buf3[3],buf3[4],buf3[5]);
+
   // fqm = fix output for peratom QM forces
   // use atomID of local atoms to index into ordered buf3
 
@@ -358,6 +360,8 @@ void FixMDIQM::post_force(int vflag)
     fqm[i][2] = buf3[3 * index + 2] * mdi2lmp_force;
   }
 
+  printf("MDI FORCE: %g %g %g: %g %g %g\n",buf3[0],buf3[1],buf3[2],buf3[3],buf3[4],buf3[5]);
+
   // optionally add forces to owned atoms
   // use atomID of local atoms to index into ordered buf3
 
@@ -371,6 +375,9 @@ void FixMDIQM::post_force(int vflag)
     }
   }
 
+  double **f = atom->f;
+  printf("LMP FORCE: %g %g %g: %g %g %g\n",f[0][0],f[0][1],f[0][2],f[1][0],f[1][1],f[1][2]);
+
   // optionally request stress tensor from MDI engine, convert to virial
   // qm_virial = fix output for global QM virial