|
|
本帖最后由 gaogaoamm 于 2025-5-7 14:05 编辑
（一）问题背景（适用计算场景）
问题一：
笔者最近在做单点能计算，因想临时更换计算精度，需要对Gaussian输入文件更换关键词，需要批量对输入文件进行编辑，因此写了一个脚本辅助处理。
问题二：
在单点能计算过程中，笔者只需要输出文件中的MP2和QCISD(T)数据，所以也写了一个python脚本进行批处理。
（二）代码介绍
2.1 关键词更换脚本（命名为replace.py放在附件中）
只需将此py文件和待处理的输入文件放在同一文件夹下，运行即可生成。
2.2 单点能后处理脚本（命名为LINE.py和MP2.py放在附件中）
只需将此py文件和待处理的输入文件放在同一文件夹下，运行即可生成。
思路说明：
（1）代码实现逻辑为直接读取文件中“MP2=”和“QCISD(T)=”以及后面带的数据，并输出到excel表中。
（2）基于逻辑与笔者代码水平，故先用LINE.py文件读取输出文件的后50行，排除了文件中其他“MP2=”和“QCISD(T)=”关键词的干扰，保证MP2.py读取的为最终计算出的单点能数据。
（3）局限性：因为有些单点能数据恰好呈现下图上图中的样式（即要读取的数据不在同一行）：所以读取会读取出空白如下图。实测下来，这种格式发生率在2%左右，可以自行打开文件输入。
附件：
LINE.py
(2.98 KB, 下载次数 Times of downloads: 3)
MP2.py
(1.28 KB, 下载次数 Times of downloads: 4)
replace.py
(1.28 KB, 下载次数 Times of downloads: 6)
|
评分 Rate
-
查看全部评分 View all ratings
|
发表于 Post on 2025-5-7 14:00:10
收藏 Add to favorites
楼主 Author