24小时热门版块排行榜    

查看: 3422  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 暖雨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 十三加油 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 一志愿070300浙大化学358分,求调剂! +3 酥酥鱼.. 2026-03-21 3/150 2026-03-22 11:31 by 杨杨杨紫
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
信息提示
请填处理意见