| 查看: 598 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
大家帮忙看看像这样的文件如何提取想要的内容。 已有2人参与
|
|||
|
- That atom : 272 E-Eref m_perp ftot fpar fperp eigen delr npart evalf a1 ( eV ) ( eV/Ang ) ( eV/Ang**2 ) 0 K= 0 0.0120 1 1 0.6542 -0.5858 0.2911 0.0000 0.054 0 2 0.00 1 K= 1 0.2074 2 2 2.2216 -2.1115 0.6907 0.0000 0.290 1 5 0.00 ... SADDLE 1001 CONVERGED |ret 20134 |delta energy= 0.0255 |force_(tot,par,perp)= 0.0854 0.0217 0.0000 |eigenval= -3.5023 |npart= 1 |delr= 0.343 |evalf= 327 | - Configuration stored in file : sad1001 - Total energy Saddle (eV) : -4.1048037236E+03 RELAXATION ... 139 M= 5 0.0116 0 0 0.2256 0.1374 0.0000 0.0000 0.377 2 333 0.00 144 M= 10 0.0093 0 0 0.2595 0.1610 0.0000 0.0000 0.380 2 338 0.00 ... - Configuration stored in file : min1001 - Total energy Minimum (eV) : -4.1052381485E+03 MINIMUM 1001 ACCEPTED |E(fin-ini)= -0.4090 |E(fin-sad)= -0.4344 |npart= 18 |delr= 1.186 |evalf= 395 | - That atom : 744 E-Eref m_perp ftot fpar fperp eigen delr npart evalf a1 ( eV ) ( eV/Ang ) ( eV/Ang**2 ) 0 K= 0 0.0110 1 1 0.5893 -0.5266 0.2644 0.0000 0.054 0 2 0.00 1 K= 1 0.1973 2 2 2.0867 -1.9790 0.6615 0.0000 0.287 0 5 0.00 2 K= 2 0.5864 3 4 3.1487 -2.5694 1.8201 0.0000 0.579 12 10 0.00 SADDLE 1002 FAILED |ret 60132 |delta energy= 0.1470 |force_(tot,par,perp)= 0.0861 0.0179 0.0000 |eigenval= 0.1918 |npart= 20 |delr= 1.050 |evalf= 288 | 我想实现，从文件中提取这些信息 1、如果saddle xxxx 显示的是CONVERGED 那么记录下当前原子标号 ，例如：that atom ： 272 并记录delta energy 例如：0.0255 2、这个文件特点是文件中saddle 显示的是CONVERGED,下面就是 RELAXATION过程 要记录E(fin-ini) ，E(fin-sad) 最后想要得到的是信息：272 0.0255 -0.4090 -0.4344 3、如果 saddle 显示的是 FAILED 那么就跳过不记录任何信息 不知道如何实现，是不是在原子标号处做上标记，在看saddle的情况，如果是CONVERGED ，就记录信息。如果不是再往下走，如何写成代码，还是没想好，希望高人指点！谢谢，我只有五个金币，全部奉上。 |
回复此楼» 猜你喜欢
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
-
环境工程 085701，267求调剂
已经有14人回复
-
308求调剂
已经有12人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有21人回复
-
311求调剂
已经有8人回复
-
294分080500材料科学与工程求调剂
已经有13人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有19人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有23人回复
-
一志愿085600中科院宁波所276分求调剂
已经有17人回复
4楼2015-07-23 08:24:37
karon1988
兑换贵宾
- 应助: 37 (小学生)
- 金币: 694.8
- 散金: 80
- 红花: 8
- 帖子: 187
- 在线: 179.4小时
- 虫号: 1690933
- 注册: 2012-03-14
- 专业: 人工智能与知识工程
2楼2015-07-22 09:56:44
jerkwin
无虫
-
专家经验: +14 - 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
3楼2015-07-22 21:07:12
5楼2015-07-23 08:42:09
