24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 350  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 324分 085600材料与化工 +20 呆鹅oor 2026-03-27 20/1000 2026-04-02 10:13 by oooqiao
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 一志愿南昌大学324求调剂 +12 hanamiko 2026-03-27 12/600 2026-04-01 22:34 by 运气yunqi
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考研] 0856,材料与化工321分求调剂 +13 大馋小子 2026-03-27 14/700 2026-04-01 15:58 by yanflower7133
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 283求调剂 +9 A child 2026-03-28 9/450 2026-04-01 14:20 by Jaylen.
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk¥
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 322求调剂:一志愿湖南大学 材料与化工(085600),已过六级。 +10 XX小邓 2026-03-29 10/500 2026-03-31 16:46 by 不吃魚的貓
[考研] 学硕274求调剂 +17 Li李鱼 2026-03-26 17/850 2026-03-31 15:19 by 客尔美德
[考研] 求收留 +8 1943443204 2026-03-28 8/400 2026-03-31 15:00 by -迷了路啊路
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
[考研] 085600,专业课化工原理,320分求调剂 +6 大馋小子 2026-03-29 6/300 2026-03-31 10:03 by 氯化亚硝酰
[考研] 一志愿华中师范化学332分求调剂 +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
信息提示
请填处理意见