24小时热门版块排行榜    

查看: 552  |  回复: 0

房然

铁虫 (小有名气)


[求助] 求助,一个bash脚本中的命令调用一个模块,然后再次输入命令的问题!!

我写了一个脚本,但是遇到了点问题:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#!/bin/bash
#$ -S /bin/bash
#$ -cwd
#echo 'date'
#$ -o /dev/null
#$ -e /dev/null
      cd /nibs02/sbfang/SPA_maker/pre_preparation/
      parmchk -i 4d8c-1.mol2 -f mol2 -o 4d8c-1-lig.frcmod
tleap
      source leaprc.gaff
      source leaprc.ff99SB
      UNK = loadmol2 4d8c-1.mol2
      loadamberparams 4d8c-1-lig.frcmod
      saveamberparm UNK 4d8c-1-lig.promtop 4d8c-1-lig.inpcrd
      DDD = loadpdb rec.pdb
      saveAmberParm DDD 4d8c-1-rec.top 4d8c-1-rec.crd
      4d8c-1-complex = combine {DDD,UNK}
      saveAmberParm 4d8c-1-complex 4d8c-1-complex.top 4d8c-1-complex.crd
quit
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
这里面 tleap 命令打开一个新的终端
-I: Adding /opt/amber12/dat/leap/prep to search path.
-I: Adding /opt/amber12/dat/leap/lib to search path.
-I: Adding /opt/amber12/dat/leap/parm to search path.
-I: Adding /opt/amber12/dat/leap/cmd to search path.

Welcome to LEaP!
(no leaprc in search path)
>

需要在新的终端中输入它下面的命令
     > source leaprc.gaff
     >  source leaprc.ff99SB
     > UNK = loadmol2 4d8c-1.mol2
     > loadamberparams 4d8c-1-lig.frcmod
     >saveamberparm UNK 4d8c-1-lig.promtop 4d8c-1-lig.inpcrd
     > DDD = loadpdb rec.pdb
     > saveAmberParm DDD 4d8c-1-rec.top 4d8c-1-rec.crd
     > 4d8c-1-complex = combine {DDD,UNK}
     > saveAmberParm 4d8c-1-complex 4d8c-1-complex.top 4d8c-1-complex.crd

如果一条一条的输入,太麻烦,可不可以说一下,怎么用以个命令让它一气呵成跑到底呢?
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 房然 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉大化学322求调剂 +3 17501029541 2026-03-23 4/200 2026-03-23 23:47 by Txy@872106
[考研] 279分求调剂 一志愿211 +17 chaojifeixia 2026-03-19 19/950 2026-03-23 23:26 by 呆呆师姐
[考研] 085600材料与化工调剂 +7 A-哆啦Z梦 2026-03-23 12/600 2026-03-23 23:16 by 星空星月
[考研] 328求调剂,英语六级551,有科研经历 +7 生物工程调剂 2026-03-17 12/600 2026-03-23 18:18 by YMU施老师
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
信息提示
请填处理意见