24小时热门版块排行榜    

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

我们都爱小木虫

相关版块跳转 我要订阅楼主 房然 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 290调剂生物0860 +28 哇哈哈,。 2026-04-11 31/1550 2026-04-13 01:16 by 幸免 ..
[考研] 322求调剂 +6 123安康 2026-04-12 13/650 2026-04-12 15:51 by 123安康
[考研] 305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[考研] 11408。358求调剂 +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by 氮气气气
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 293求调剂 +6 勇远库爱314 2026-04-08 6/300 2026-04-11 10:08 by zhq0425
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 考研调剂 +13 冰冰,,, 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 生物学学硕,初试351分,求调剂 +4 …~、王…~ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
信息提示
请填处理意见