24小时热门版块排行榜    

查看: 3421  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 暖雨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 求调剂 +6 十三加油 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
信息提示
请填处理意见