Update the reference file for test

Author: ttadano

example/ZnO/qha_relax/reference/ZnO_qha.qha_thermo

@@ -1,13 +1,13 @@
 # The renormalized static potential Phi_0 is also shown.
-# Temperature [K], Cv [in kB unit], F_{vib} (QHA term) [Ry], Phi0 [Ry], F_{total} [Ry]
-        0.000000      0.000000e+00      1.693972e-02      7.585253e-05      1.701558e-02
-      100.000000      4.453944e+00      1.625895e-02      7.109735e-05      1.633005e-02
-      200.000000      7.903822e+00      1.295553e-02      7.522584e-05      1.303076e-02
-      300.000000      9.708692e+00      7.141166e-03      9.678675e-05      7.237953e-03
-      400.000000      1.058941e+01     -7.498083e-04      1.331007e-04     -6.167076e-04
-      500.000000      1.105796e+01     -1.034688e-02      1.835930e-04     -1.016329e-02
-      600.000000      1.133091e+01     -2.137413e-02      2.488419e-04     -2.112529e-02
-      700.000000      1.150209e+01     -3.362757e-02      3.299940e-04     -3.329757e-02
-      800.000000      1.161588e+01     -4.695310e-02      4.284902e-04     -4.652461e-02
-      900.000000      1.169509e+01     -6.123165e-02      5.461399e-04     -6.068551e-02
-     1000.000000      1.175232e+01     -7.636875e-02      6.846873e-04     -7.568406e-02
+# Temperature [K], Cv [in kB unit], F_{vib} (QHA term) [Ry], Phi0 [Ry], F_{total} [Ry], S_{vib} [in kB unit]
+        0.000000      0.000000e+00      1.693972e-02      7.585253e-05      1.701558e-02     -0.000000e+00
+      100.000000      4.453944e+00      1.625895e-02      7.109735e-05      1.633005e-02      3.045469e+00
+      200.000000      7.903822e+00      1.295553e-02      7.522584e-05      1.303076e-02      7.288478e+00
+      300.000000      9.708692e+00      7.141166e-03      9.678675e-05      7.237953e-03      1.088367e+01
+      400.000000      1.058941e+01     -7.498083e-04      1.331007e-04     -6.167076e-04      1.382049e+01
+      500.000000      1.105796e+01     -1.034688e-02      1.835930e-04     -1.016329e-02      1.625110e+01
+      600.000000      1.133091e+01     -2.137413e-02      2.488419e-04     -2.112529e-02      1.830810e+01
+      700.000000      1.150209e+01     -3.362757e-02      3.299940e-04     -3.329757e-02      2.008522e+01
+      800.000000      1.161588e+01     -4.695310e-02      4.284902e-04     -4.652461e-02      2.164735e+01
+      900.000000      1.169509e+01     -6.123165e-02      5.461399e-04     -6.068551e-02      2.304026e+01
+     1000.000000      1.175232e+01     -7.636875e-02      6.846873e-04     -7.568406e-02      2.429699e+01