24小时热门版块排行榜    

查看: 393  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见