24小时热门版块排行榜    

查看: 3359  |  回复: 2

暖雨

铁虫 (小有名气)

[求助] 求助:lammps怎么输出每一步的能量?

我用lammps模拟反应,想要输出每一步的能量,也就是下面TotEng后面那个数值,请问应该怎么做到啊?是在in文件里添加某条命令实现吗?多谢了,各位大神~~~
---------------- Step   160544 ----- CPU =      6.2692 (sec) ----------------
TotEng   =    -29885.2633 KinEng   =      2676.3909 Temp     =      2302.2402
PotEng   =    -32561.6542 E_bond   =         0.0000 E_angle  =         0.0000
E_dihed  =         0.0000 E_impro  =         0.0000 E_vdwl   =    -29138.8029
E_coul   =     -3422.8512 E_long   =         0.0000 Press    =      6157.4715
---------------- Step   160545 ----- CPU =      6.5220 (sec) ----------------
TotEng   =    -29893.5418 KinEng   =      2676.3909 Temp     =      2302.2402
PotEng   =    -32569.9327 E_bond   =         0.0000 E_angle  =         0.0000
E_dihed  =         0.0000 E_impro  =         0.0000 E_vdwl   =    -29146.2070
E_coul   =     -3423.7257 E_long   =         0.0000 Press    =      6043.8727
---------------- Step   160546 ----- CPU =      6.7705 (sec) ----------------
TotEng   =    -29901.9510 KinEng   =      2676.3909 Temp     =      2302.2402
PotEng   =    -32578.3419 E_bond   =         0.0000 E_angle  =         0.0000
E_dihed  =         0.0000 E_impro  =         0.0000 E_vdwl   =    -29153.3749
E_coul   =     -3424.9670 E_long   =         0.0000 Press    =      6439.5617
---------------- Step   160547 ----- CPU =      7.0162 (sec) ----------------
TotEng   =    -29910.3065 KinEng   =      2676.3909 Temp     =      2302.2402
PotEng   =    -32586.6973 E_bond   =         0.0000 E_angle  =         0.0000
E_dihed  =         0.0000 E_impro  =         0.0000 E_vdwl   =    -29159.7367
E_coul   =     -3426.9606 E_long   =         0.0000 Press    =      6080.9623
---------------- Step   160548 ----- CPU =      7.2654 (sec) ----------------
TotEng   =    -29917.4197 KinEng   =      2676.3909 Temp     =      2302.2402
PotEng   =    -32593.8106 E_bond   =         0.0000 E_angle  =         0.0000
E_dihed  =         0.0000 E_impro  =         0.0000 E_vdwl   =    -29164.7130
E_coul   =     -3429.0975 E_long   =         0.0000 Press    =      6763.0574
---------------- Step   160549 ----- CPU =      7.4530 (sec) ----------------
TotEng   =    -29924.6364 KinEng   =      2676.3909 Temp     =      2302.2402
PotEng   =    -32601.0273 E_bond   =         0.0000 E_angle  =         0.0000
E_dihed  =         0.0000 E_impro  =         0.0000 E_vdwl   =    -29168.1781
E_coul   =     -3432.8492 E_long   =         0.0000 Press    =      5433.0208
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

issp-hao

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
暖雨: 金币+5, ★★★很有帮助 2013-10-28 15:03:54
月只蓝: 金币+3, 感谢指导! 2013-10-29 14:43:09
我觉得可以考虑把in文件中原来的那个thermo——style前面加#
后面重写一个只有能量的
比如:
thermo_style custom step etotal
就会只输出总能量了
2楼2013-10-24 08:47:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

botar

木虫 (正式写手)


月只蓝: 金币+1, 感谢指导! 2013-10-29 14:43:16
用matlab好像可以提取出来,并绘图
努力,唯有努力。
3楼2013-10-24 11:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 暖雨 的主题更新
信息提示
请填处理意见