diff --git a/source/module_md/md_func.cpp b/source/module_md/md_func.cpp
index 3d40ca3f7e..07bc1f03f2 100644
--- a/source/module_md/md_func.cpp
+++ b/source/module_md/md_func.cpp
@@ -453,6 +453,7 @@ double current_temp(double& kinetic,
 {
     if (3 * natom == frozen_freedom)
     {
+        kinetic = 0.0;
         return 0.0;
     }
     else