个人中心
文章发布
退出
作者热门文章
- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
24
4
<分区>
我需要一些帮助来运行这段代码。我从 ( http://easybioinfo.free.fr/?q=content/amber-trajectory-gromacs-xtc-conversion) 中获取了这段代码。我正在尝试将 amber trajectory 转换为 gromacs trajectory。
当我执行这段代码时，出现了一些错误。我将错误粘贴到此代码下方:
#!/usr/bin/python
#Workflow based on Trajectory Converter - v1.5 by: Justin Lemkul
#completely reimplemented and improved by Peter Schmidtke & Jesus Seco
import sys,os,re,fnmatch
if len(sys.argv)>4 :
f=sys.argv[1]
if not os.path.exists(f):
sys.exit(" ERROR : Something you provided does not exist. Breaking up.\n\nUSAGE : python trajconv_peter.py amberCrd amberTop trajDir trajPattern outPutPrefix\n\n \
Example : python amber2xtc.py mdcrd.crd mdcrd.top md *.x.gz md_gromacs\n")
else :
sys.exit(" \n USAGE : python amber2xtc.py AMBERCRD AMBERTOP TRAJDIR TRAJPATTERN OUTPUTPREFIX\n\
Example : python amber2xtc.py mdcrd.crd mdcrd.top md *.x.gz md_gromacs\n\
Note that the AmberCrd can also be a PDB file.\n")
crd=sys.argv[1]
top=sys.argv[2]
trajdir=sys.argv[3]
pattern=sys.argv[4]
outputPref=sys.argv[5]
traj_files=fnmatch.filter(os.listdir(trajdir),pattern) #get the fpocket output folders
RE_DIGIT = re.compile(r'(\d+)') #set a pattern to find digits
ALPHANUM_KEY = lambda s: [int(g) if g.isdigit() else g for g in RE_DIGIT.split(s)] #create on the fly function (lambda) to return numbers in filename strings
traj_files.sort(key=ALPHANUM_KEY) #sort by these numbers in filenames
print "Will convert the following files : "
print traj_files
csn=1
for file in traj_files :
ptrajtmp=open("ptraj_tmp.ptr","w")
print "currently converting "+file
ptrajtmp.write("trajin "+trajdir+os.sep+file+"\n")
ptrajtmp.write("reference "+crd+"\n")
ptrajtmp.write("center ~:WAT,CIO mass origin\n")
ptrajtmp.write("image origin center :* byres familiar\n")
ptrajtmp.write("trajout pdb_tmp/mdcrd.pdb pdb")
ptrajtmp.close()
if not os.path.exists("pdb_tmp"):
os.mkdir("pdb_tmp")
os.system("ptraj "+top +" ptraj_tmp.ptr >/dev/null 2>&1")
if not os.path.exists("xtc_tmp"):
os.mkdir("xtc_tmp")
#move to *.pdb
os.system("cd pdb_tmp; ls *.pdb.* | cut -f3 -d\".\" | awk '{print \"mv mdcrd.pdb.\"$0\" mdcrd_\"$0\".pdb\" }' | sh ; cd ../")
pdb_files=fnmatch.filter(os.listdir("pdb_tmp"),"*.pdb")
pdb_files.sort(key=ALPHANUM_KEY) #sort by these numbers in filenames
if csn==1:
os.system("editconf -f pdb_tmp/mdcrd_1.pdb -o "+outputPref+"_t1_top.gro >/dev/null 2>&1")
for pdb in pdb_files:
os.system("echo \"0\" | trjconv -s pdb_tmp/"+pdb+" -f pdb_tmp/"+pdb+" -o xtc_tmp/traj_"+str(csn)+".pdb.xtc -t0 "+str(csn)+" >/dev/null 2>&1")
csn+=1
if os.path.exists(outputPref+"_traj.xtc"):
os.system("trjcat -f "+outputPref+"_traj.xtc xtc_tmp/*.pdb.xtc -o "+outputPref+"_traj.xtc >& trajcat.log")
else :
os.system("trjcat -f xtc_tmp/*.pdb.xtc -o "+outputPref+"_traj.xtc >& trajcat.log")
os.system("rm -rf pdb_tmp/*.pdb")
os.system("rm -rf xtc_tmp/*.xtc")
os.remove("ptraj_tmp.ptr")
os.system("rmdir pdb_tmp")
os.system("rmdir xtc_tmp")
错误如下:
vijay@glycosim:~/Simulation-Folder-Feb2013/chapter5-thermo-paper2-Vj/analysis-malto-/28-difusion-coeff-malto-thermo/convert-gromacs-format$ python2.7 amber2xtc.py malto-THERMO.crd malto-THERMO.top TRAJDIR malto*.traj md_gromacss
Will convert the following files :
['malto-thermo.set11.traj', 'malto-thermo.set12.traj', 'malto-thermo.set13.traj', 'malto-thermo.set14.traj', 'malto-thermo.set15.traj']
currently converting malto-thermo.set11.traj
ls: cannot access *.pdb.*: No such file or directory
sh: Syntax error: Bad fd number
currently converting malto-thermo.set12.traj
ls: cannot access *.pdb.*: No such file or directory
sh: Syntax error: Bad fd number
currently converting malto-thermo.set13.traj
ls: cannot access *.pdb.*: No such file or directory
sh: Syntax error: Bad fd number
currently converting malto-thermo.set14.traj
ls: cannot access *.pdb.*: No such file or directory
sh: Syntax error: Bad fd number
currently converting malto-thermo.set15.traj
ls: cannot access *.pdb.*: No such file or directory
sh: Syntax error: Bad fd number
vijay@glycosim:~/Simulation-Folder-Feb2013/chapter5-thermo-paper2-Vj/analysis-malto-/28-difusion-coeff-malto-thermo/convert-gromacs-format$
有关信息，我使用的是 Ubuntu 11.10(64 位)。
如何纠正这个错误？感谢任何帮助。谢谢。
24
4
0
我发现了以下区别: catalina.sh 运行 ---> 在前台启动 tomcat。 startup.sh ---> 在写入响应时在后台启动 tomcat在日志文件中。 我找不到与此命令的比较: t
这个问题已经有答案了: 已关闭11 年前。 Possible Duplicate: Whats the difference between running a shell script as ./s
我正在尝试使用 bash 和数组。但是执行示例脚本时，我收到了意外的语法错误消息:example.sh: 3: example.sh: Syntax error: "("unexpected。这是脚本
我有文件 config.sh: data_dir=/home/olle/investing_data bin_dir=/home/olle/investing_bin 我想在两个单独的脚本a.sh 和
我正在处理一个 sh 脚本，其中我处于 WHILE 循环中，其中每次迭代都会增加一个变量，并且我想每五次迭代保存一个文件。 我通常会做的(用 C 语言说)是做一个 if ( n % 5 == 0)测试
这个 STR="Hello\nWorld" echo $STR 作为输出产生 Hello\nWorld 而不是 Hello World 如何在字符串中添加换行符？ 注意:这个问题与echo无关。我知道
这个 STR="Hello\nWorld" echo $STR 作为输出产生 Hello\nWorld 代替 Hello World 我应该怎么做才能在字符串中换行？ 注意:这个问题不是关于echo的
我有一个看起来像这样的脚本 #!/bin/bash function something() { echo "hello world!!" } something | tee logfile 我已经
在这里，我有一个名为 load.sh 的 shell 脚本。 它用监督启动我名为“demo”的程序，当我用 sh load.sh start | 运行它时停止，它运行良好。 但是，当我使用 ./loa
我刚刚从官方网站下载了 robomongo 64 位的 ubuntu。然后我打开我的终端并查看文件夹结构并从 bin 文件夹执行命令。 sh robomongo.sh 我遇到这样的错误 robomon
这个问题已经有答案了: What is the difference between `./example.sh` and `sh example.sh` (3 个回答) 已关闭 8 年前。 使用 .
我正在尝试从 c 程序执行在 Linux 中启动新终端的命令 sh，但我收到此错误/bin/sh:/bin/sh: 当我使用 execv () 时无法执行二进制文件结果是 fork，然后是 execv
我想在我的 sh 脚本中使用数组 我的目标是为以下值创建数组 - a b c并按数组打印所有值 我成功打印每个数组，但未能打印数组中的所有值 请看下面的例子 设置arr中的每个值 n=1 eval a
我基本上正在寻找一种方法来做到这一点 list=$(command) while read -r arg do ... done output & while read -r arg
我部署了一个 Surge.sh 项目并成功发布，但是，我想对该项目进行一些更新，更新已在本地应用并且可以看到工作，但它尝试再次将其发布到同一个域但它没有'没有出现更新。我还尝试拆除该项目并重新上传它仍
我在 CentOS 上运行调用远程 sh 文件的 sh 脚本时遇到问题。在 user@host1 我有 start.sh 文件，里面有下一个命令 NODE1_SSH_PATH=user@host2 P
我是surge.sh plus客户，拥有自定义域，比方说 http://bobswidgets.com 但是，Bob 要求我为他的新定制服务设置一个子域。他要custom.bobswidgets.co
$ sh sh-3.2$ if > ps -ef | grep apple ; > then > echo APPLE > fi ; lazer 7584 7571 0 04:36 pts/4
我正在使用 bash shell 导出一些路径，如下所示: 环境.sh: export GOPATH=$GOPATH:$QBOXROOT/pili-zeus export PATH=$PATH:$
我在采购别名方面遇到了问题。这确实是具体案例。我使用 open_new.sh 打开带有两个选项卡的新终端窗口: #!/bin/bash gnome-terminal --tab --t
我是一名优秀的程序员,十分优秀！