24小时热门版块排行榜    

Znn3bq.jpeg
查看: 365  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +6 popular289 2026-05-18 11/550 2026-05-18 22:40 by iaeyuan
[考博] 博士申请 +4 星…… 2026-05-18 5/250 2026-05-18 17:34 by 炎甲00
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见