24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 555  |  回复: 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)

我们都爱小木虫

相关版块跳转 我要订阅楼主 房然 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +8 Xxjc1107. 2026-04-13 8/400 2026-04-13 14:47 by ajpv风雷
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 材料工程085601,270求调剂 +38 @ASDF1234 2026-04-08 41/2050 2026-04-13 12:53 by 张zhihao
[考研] 2026硕士调剂_能动_河南农业大学 +3 河南农业大学-能 2026-04-12 3/150 2026-04-13 12:44 by 木鱼.
[考研] 一志愿0703化学招61最终排名62化学求调剂 +29 招61排名62 2026-04-07 33/1650 2026-04-13 10:16 by sxdj2
[考研] 268分085602化学工程调剂 +31 月照花林。 2026-04-09 31/1550 2026-04-12 23:37 by 学zh
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 347求调剂 +4 mhyqyy 2026-04-06 4/200 2026-04-12 02:27 by 秋豆菜芽
[考研] 计算机22408 281分,求调剂 +7 17715607211 2026-04-06 7/350 2026-04-12 00:43 by xuxiang
[考研] 求调剂 +6 小聂爱学习 2026-04-11 9/450 2026-04-11 21:20 by 蓝云思雨
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 085500求调剂材料 +10 易11122 2026-04-09 10/500 2026-04-11 10:39 by maddjdld
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
信息提示
请填处理意见