lammps 2个温度修正命令，使模拟结果更精确

大家好，我是小马老师。
根据分子动力学理论，体系的温度与原子的速度有关。
温度由公式 KE = dim/2NKT计算得到
其中：KE为原子组的总动能（1/2 mv^2)，dim为模拟维度，N为原子总数，K为玻尔兹曼常数，T为温度。
在一些模拟中，体系的温度并不能完全按照以上公式进行计算。
例如，在lammps流体模拟中，流体原子按照一定的速度流动，如果把流体流动速度加入到上述公式进行计算，得到的温度可能与实际温度不符。
一般来说，流体的温度计算公式中应该去掉流体的流动速度。
在lammps中，有两个命令可实现这个功能：
（1）compute temp/partial
该命令语法为：

compute ID group-ID temp/partial xflag yflag zflag

其中，group-ID为需要计算温度的原子组。
xflag yflag zflag确定是否计算该方向的温度。取值为1表示计算温度是包含该方向的速度，取值为0表示不包含该方向的速度。