菜鸡论文模仿03（滑动摩擦）：帮我检查一下错误呗！！！！！！！

1Gpa下1m/s的滑动摩擦

加载过程中，我的F控温和F不控温状态下，温度都能达到一个平衡值，但是一到滑动就会出现EF，CD层温度下降的趋势，而且也不能最终平衡到某一温度附近，不知道是我的代码原因，还是速度太小？
代码需要实现的功能：1Gpa下的1m/s的滑动摩擦
各位帮我看看呗，哪儿有问题，私信或邮箱我，谢谢！！！！！！！！

#-----------------------EF速度恒定
velocity		E_F set NULL ${input_velocity} NULL sum yes

#-----------------------B层恒温300K
fix			1 B nvt temp 300.0 300.0 50.0
reset_timestep		0
timestep		0.5
run			1000

#-----------------------C_D层nve
fix			2 C_D nve
fix			3 C_D wall/reflect zhi ${reflect2} units box	#尽量保持密度一定

#-----------------------E_F层nve
fix			4 E_F nve/limit 0.05
fix			5 E_F addforce 0.0 0.0 -${f_to_addforce}					
fix			6 F temp/rescale 10 300.0 300.0 10.0 1.0
reset_timestep		0
timestep		0.3
run			140000