个人中心
文章发布
退出
作者热门文章
- 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
reqwest v0.9 将 serde v1.0 作为依赖项，因此实现 converting serde_json errors into reqwest error . 在我的代码中，我使用 se
我有这个代码: let file = FileStorage { // ... }; file.write("Test", bytes.as_ref()) .map_err(|e| Mu
我只是尝试用angular-cli创建一个新项目，然后运行服务器，但是它停止并显示一条有趣的消息:Error: No errors。 我以这种方式更新了(希望有帮助):npm uninstall -g
我从我的 javascript 发送交易 Metamask 打开传输对话框 我确定 i get an error message in metamask (inpage.js:1 MetaMask -
这个问题在这里已经有了答案: How do you define custom `Error` types in Rust? (3 个答案) How to get a reference to a
我想知道两者之间有什么大的区别 if let error = error{} vs if error != nil?或者只是人们的不同之处，比如他们如何用代码表达自己？ 例如，如果我使用这段代码: u
当我尝试发送超过 50KB 的图像时，我在 Blazor 服务器应用程序上收到以下错误消息 Error: Connection disconnected with error 'Error: Serv
我有一个error-page指令，它将所有异常重定向到错误显示页面 我的web.xml: [...] java.lang.Exception /vi
我有这样的对象: address: { "phone" : 888, "value" : 12 } 在 WHERE 中我需要通过 address.value 查找对象，但是在 SQL 中有函数
每次我尝试编译我的代码时，我都会遇到大量错误。这不是我的代码的问题，因为它在另一台计算机上工作得很好。我尝试重新安装和修复，但这没有帮助。这是整个错误消息: 1>------ Build starte
在我的代码的类部分，如果我写一个错误，则在不应该的情况下，将有几行报告为错误。我将'| error'放在可以从错误中恢复的良好/安全位置，但是我认为它没有使用它。也许它试图在某个地方恢复中间表情？ 有
我遇到了 csv 输入文件整体读取故障的问题，我可以通过在 read_csv 函数中添加 "error_bad_lines=False" 来删除这些问题来解决这个问题。 但是我需要报告这些造成问题的文
在 Spring 中，验证后我们在 controller 中得到一个 BindingResult 对象。 很简单，如果我收到验证错误，我想重新显示我的表单，并在每个受影响的字段上方显示错误消息。 因此
我不知道出了什么问题，因为我用 Java 编程了大约一年，从来没有遇到过这个错误。在一分钟前在 Eclipse 中编译和运行工作，现在我得到这个错误: #A fatal error has been
SELECT to_char(messages. TIME, 'YYYY/MM/DD') AS FullDate, to_char(messages. TIME, 'MM/DD
我收到这些错误: AnonymousPath\Anonymized.vb : error BC30037: Character is not valid. AnonymousPath\Anonymiz
我刚刚安装了 gridengine 并在执行 qstat 时出现错误: error: commlib error: got select error (Connection refused) erro
嗨，我正在学习 PHP，我从 CRUD 系统开始，我在 Windows 上安装了 WAMP 服务器，当我运行它时，我收到以下错误消息。 SCREAM: Error suppression ignore
我刚刚开始一个新项目，我正在学习核心数据教程，可以找到:https://www.youtube.com/watch?v=zZJpsszfTHM 我似乎无法弄清楚为什么会抛出此错误。我有一个名为“Exp
当我使用 Jenkins 运行新构建时，出现以下错误: "FilePathY\XXX.cpp : fatal error C1853: 'FilePathZ\XXX.pch' precompiled
我是一名优秀的程序员,十分优秀！