24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3457  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 暖雨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程281还有调剂机会吗 +28 xaw. 2026-04-11 29/1450 2026-04-12 17:55 by 郁郁菲菲
[考研] 药学求调剂 +3 RussHu 2026-04-12 4/200 2026-04-12 17:49 by 陈皮皮
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 环境专硕调剂 +16 会说话的肘子 2026-04-06 16/800 2026-04-10 10:30 by asy1wn
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
信息提示
请填处理意见