分子動力学機能のBerendesen thermostatにおいて、自由度の計算がNを原子数としたとき

3*N-3でなければならないところが3*(N-3)になっていたバグが発見されました。

このバグは設定した温度Tに対して実際の温度Tcalが

(3*N-3)/2*k_B*Tcal = 3*(N-3)/2*k_B*T

とされる影響を及ぼします。Nが大きいとその分影響は小さくなります。

140129版のsrc/moldyn.F90を添付のファイルで差し替えてください。