24小时热门版块排行榜    

查看: 455  |  回复: 1

leon583520

新虫 (初入文坛)

[求助] 求助 我的脚本运行不了,求高手指点

我写的脚本如下:
#!/bin/bash -x
# Two Step RESP script
#cp AAadd.dat add.dat
espgen -i $1_resp.out -o $1.esp
antechamber -fi gout -fo ac -i $1_resp.out -o $1.ac -j $2 > antechamber.out
respgen -i $1.ac -o RESP.INPUT1 -a add.dat -f resp1 -n 1
respgen -i $1.ac -o RESP.INPUT2 -a add.dat -f resp2 -n 1

sed '/ioutopt/ a\
  qwt = 0.005,
  ' RESP.INPUT1 > RESP.INPUT1a

resp -O -i RESP.INPUT1a -o RESP.OUTPUT1a -q QIN -e $1.esp
#resp -O -i RESP.INPUT1 -o RESP.OUTPUT1 -q QIN -e $1.esp
resp -O -i RESP.INPUT2 -o RESP.OUTPUT2 -q qout -e $1.esp
antechamber -fi ac -fo ac -i $1.ac -o $1_out.ac -c rc -cf qout -at amber -j $2 >> antechamber.out
cp RESP.OUTPUT2 $1.respout
#tar -czf $1.tgz  RESP* $1*.out $1*.ac $1.respout
grep "Dipole Moment" RESP.OUTPUT2
grep ' Tot=' $1_resp.out

grep ATOM $1_out.ac |cut -c 55-69 > 5k.q
cat 5k.q

然后把高斯的输出文件和add.dat文件一并拷贝到一个文件夹里,运行脚本后,出现到结果是:

+ espgen -i HFO-1234yf00._resp.out -o HFO-1234yf00..esp

Cannot open the input file HFO-1234yf00._resp.out, exit+ antechamber -fi gout -fo ac -i HFO-1234yf00._resp.out -o HFO-1234yf00..ac -j
Cannot open the input file HFO-1234yf00._resp.out in rgout(), exit
+ respgen -i HFO-1234yf00..ac -o RESP.INPUT1 -a add.dat -f resp1 -n 1
Cannot open the ac input file: HFO-1234yf00..ac, exit
+ respgen -i HFO-1234yf00..ac -o RESP.INPUT2 -a add.dat -f resp2 -n 1
Cannot open the ac input file: HFO-1234yf00..ac, exit
+ sed '/ioutopt/ a\
  qwt = 0.005,
  ' RESP.INPUT1
sed: can't read RESP.INPUT1: No such file or directory
+ resp -O -i RESP.INPUT1a -o RESP.OUTPUT1a -q QIN -e HFO-1234yf00..esp
At line 334 of file resp.f (unit = 5, file = 'RESP.INPUT1a')
Fortran runtime error: End of file
+ resp -O -i RESP.INPUT2 -o RESP.OUTPUT2 -q qout -e HFO-1234yf00..esp

  Unit    5 Error on OPEN: RESP.INPUT2                                                                     
+ antechamber -fi ac -fo ac -i HFO-1234yf00..ac -o HFO-1234yf00._out.ac -c rc -cf qout -at amber -j
Cannot open the ac input file: HFO-1234yf00..ac, exit
+ cp RESP.OUTPUT2 HFO-1234yf00..respout
+ grep 'Dipole Moment' RESP.OUTPUT2
+ grep ' Tot=' HFO-1234yf00._resp.out
grep: HFO-1234yf00._resp.out: No such file or directory
+ grep ATOM HFO-1234yf00._out.ac
+ cut -c 55-69
grep: HFO-1234yf00._out.ac: No such file or directory
+ cat 5k.q
zhengxy@zhengxy-ThinkCentre-XXXX:~/src/amber10/HFO-1234yf_resp/pop=mk$ ./one5k.sh HFO-1234yf00
+ espgen -i HFO-1234yf00_resp.out -o HFO-1234yf00.esp

Cannot open the input file HFO-1234yf00_resp.out, exit+ antechamber -fi gout -fo ac -i HFO-1234yf00_resp.out -o HFO-1234yf00.ac -j
Cannot open the input file HFO-1234yf00_resp.out in rgout(), exit
+ respgen -i HFO-1234yf00.ac -o RESP.INPUT1 -a add.dat -f resp1 -n 1
Cannot open the ac input file: HFO-1234yf00.ac, exit
+ respgen -i HFO-1234yf00.ac -o RESP.INPUT2 -a add.dat -f resp2 -n 1
Cannot open the ac input file: HFO-1234yf00.ac, exit
+ sed '/ioutopt/ a\
  qwt = 0.005,
  ' RESP.INPUT1
sed: can't read RESP.INPUT1: No such file or directory
+ resp -O -i RESP.INPUT1a -o RESP.OUTPUT1a -q QIN -e HFO-1234yf00.esp
At line 334 of file resp.f (unit = 5, file = 'RESP.INPUT1a')
Fortran runtime error: End of file
+ resp -O -i RESP.INPUT2 -o RESP.OUTPUT2 -q qout -e HFO-1234yf00.esp

  Unit    5 Error on OPEN: RESP.INPUT2                                                                     
+ antechamber -fi ac -fo ac -i HFO-1234yf00.ac -o HFO-1234yf00_out.ac -c rc -cf qout -at amber -j
Cannot open the ac input file: HFO-1234yf00.ac, exit
+ cp RESP.OUTPUT2 HFO-1234yf00.respout
+ grep 'Dipole Moment' RESP.OUTPUT2
+ grep ' Tot=' HFO-1234yf00_resp.out
grep: HFO-1234yf00_resp.out: No such file or directory
+ grep ATOM HFO-1234yf00_out.ac
+ cut -c 55-69
grep: HFO-1234yf00_out.ac: No such file or directory
+ cat 5k.q

请教下这是什么原因呢,我设置了PATH环境变量,把AMBER的bin路径也加了进来,还是算不出来,小弟菜鸟,求高手指点,不胜感激

另:我到高斯输出文件名是HFO-1234yf00.out

[ Last edited by ghcacj on 2011-8-16 at 06:50 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

刚开始做分子模拟,正在努力,多多交流
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leon583520

新虫 (初入文坛)

为什么木人理偶,是不是悬赏少了,我就有那一金币呢......
刚开始做分子模拟,正在努力,多多交流
2楼2011-08-15 22:06:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leon583520 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-11 4/200 2026-11-12 02:56 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Z7qfCsdIiD1T 2026-11-11 4/200 2026-11-12 02:53 by 9fOTqj074tIg
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 02:23 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[基金申请] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 02:14 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 02:01 by 9fOTqj074tIg
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 3/150 2026-11-12 01:50 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 5/250 2026-11-12 01:17 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 00:55 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 00:55 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 3/150 2026-11-11 18:26 by TgvkbZZl9tqB
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-11 17:25 by TgvkbZZl9tqB
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:09 by TgvkbZZl9tqB
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 0TEMgxdT1eIw 2026-11-09 7/350 2026-11-11 16:58 by TgvkbZZl9tqB
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-11 16:41 by TgvkbZZl9tqB
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:41 by gRND4NVAZ8VG
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
信息提示
请填处理意见