24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 351  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +9 汲亦昊 2026-03-29 9/450 2026-04-02 11:44 by 运气yunqi
[考研] 08开头看过来!!! +4 wwwwffffff 2026-03-31 6/300 2026-04-02 11:42 by 均值回归
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 265求调剂11408 +4 刘小鹿lu 2026-03-27 4/200 2026-04-02 11:27 by guyan1000
[考研] 311求调剂 +13 蓝月亮亮 2026-03-30 13/650 2026-04-02 11:26 by wanghuan83214
[考研] 337求调剂 +11 《树》 2026-03-29 11/550 2026-04-02 10:20 by 不吃魚的貓
[考研] 求调剂,一志愿南京师范大学计算机专硕,初试373,六级通过, +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 【求调剂】085601材料工程专硕 | 总分272 | +10 脚滑的守法公民 2026-03-27 10/500 2026-04-01 17:23 by pies112
[考研] 330分求调剂 +11 qzenlc 2026-03-29 11/550 2026-04-01 14:32 by chenqifeng666
[考研] 085600,321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 材料与化工调剂一志愿大连海事085600,349 +9 吃的不少 2026-03-30 9/450 2026-04-01 11:24 by wangjy2002
[考研] 085410人工智能 初试316分 求调剂 +3 残星拂曙 2026-03-31 3/150 2026-04-01 11:09 by 小熊raider
[考研] 322求调剂 +8 三水sss 2026-04-01 8/400 2026-04-01 10:19 by 唐沐儿
[考研] 352分-085602-一志愿985 +6 海纳百川Ly 2026-03-29 6/300 2026-03-31 21:06 by yuq
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 304材料求调剂 +5 钟llll 2026-03-26 5/250 2026-03-30 19:37 by 源_2020
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 16:05 by shuang5186
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 一志愿上海理工能源动力(085800)310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
信息提示
请填处理意见