24小时热门版块排行榜    

查看: 551  |  回复: 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)

我们都爱小木虫

相关版块跳转 我要订阅楼主 房然 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +5 某某某某位 2026-03-21 5/250 2026-03-23 23:36 by 热情沙漠
[考研] 一志愿华东理工大学081700,初试分数271 +4 kotoko_ik 2026-03-23 5/250 2026-03-23 23:23 by 呆呆师姐
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 352求调剂 +3 大米饭! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考研] 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 人来盛
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
信息提示
请填处理意见