计算化学公社

 找回密码 Forget password
 注册 Register
Views: 10486|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 如何实现文件名随运行次数同步更新

[复制链接 Copy URL]

20

帖子

0

威望

589

eV
积分
609

Level 4 (黑子)

各位老师好,我想做一个多重循环的MD模拟,初始的结构文件为0.gro, 每步MD的运行需要调用上一步MD新产生的结构文件,想请教一下如何编写shell脚本可以实现产生的结构文件[step].gro能够随着运行的step数同步更改,从而可以被下个MD过程调用?请各位老师指导,谢谢!

313

帖子

2

威望

3900

eV
积分
4253

Level 6 (一方通行)

2#
发表于 Post on 2021-4-4 23:43:48 | 只看该作者 Only view this author

2425

帖子

1

威望

6196

eV
积分
8641

Level 6 (一方通行)

3#
发表于 Post on 2021-4-5 09:11:27 | 只看该作者 Only view this author
app.exe < INPUT > OUTPUT
保持不变。
while 循环,
counter递增

采用 output-$counter.gro之类的文件名。
然后用cp 或者ln创建每次文件名。

我自己会这么做。

会有更好的方案。
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

13

帖子

0

威望

656

eV
积分
669

Level 4 (黑子)

4#
发表于 Post on 2021-4-9 10:57:24 | 只看该作者 Only view this author
#!/bin/sh
for((i=1; i<43; i++));
do
j=$[i+1];
gro1=“..." 此处结构文件路径
if [ -f "$gro1" ]; 判断是否存在结构文件
then
gmx1="gmx grompp -f def.mdp -p si.top -c def"$i".gro -t def"$i".cpt -o def"$j".tpr" 此处输入gmx grompp命令,字符串组合
echo | $gmx1
tpr="....." 此处生成的tpr文件路径
if [ -f "$tpr" ]; 判断是否存在生成tpr文件
then
md1=".....” 输入gmx mdrun命令
echo | $md1
fi
done
我是这么写的,希望能帮到你

20

帖子

0

威望

589

eV
积分
609

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-4-9 14:35:03 | 只看该作者 Only view this author
wwj 发表于 2021-4-9 10:57
#!/bin/sh
for((i=1; i

学习到,非常感谢您的指导

20

帖子

0

威望

589

eV
积分
609

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-4-9 14:35:36 | 只看该作者 Only view this author
abin 发表于 2021-4-5 09:11
app.exe < INPUT > OUTPUT
保持不变。
while 循环,

感谢老师指导!

20

帖子

0

威望

589

eV
积分
609

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2021-4-9 14:36:37 | 只看该作者 Only view this author
lyj714 发表于 2021-4-4 23:43
http://bbs.keinsci.com/thread-22116-1-1.html

明白了,谢谢您的回复

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 20:53 , Processed in 0.277839 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list