<div dir="ltr"><div><div>Dear Elaine and Eric:<br><br>Thank you for your help.<br></div>Actually the following previous posts were directly helpful to accomplish what I wanted to do, i.e. utilizing MoveDialog.<br></div><div>
<div><a href="http://www.cgl.ucsf.edu/pipermail/chimera-dev/2013/000968.html" target="_blank">http://www.cgl.ucsf.edu/pipermail/chimera-dev/2013/000968.html</a><br><a href="http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-August/006634.html" target="_blank">http://plato.cgl.ucsf.edu/pipermail/chimera-users/2011-August/006634.html</a><br>
<br></div><div>Now hbond information seems properly saved upon updated frame with LoadFrame.<br><br></div><div>Thanks a lot!<br></div><div><br>#chimera --script Hbond.py<br>#<br>from chimera import runCommand<br>from Movie.gui import MovieDialog<br>
from chimera.extension import manager<br><br># pr trajectory movie open by metafile<br>meta = open("metafile", 'w')<br>
meta_input = """ amber<br>mypdb.prmtop<br>pr.crd<br>"""<br>meta.write(meta_input)<br>meta.close()<br>#####<br>def findMDMovie():<br> mdms = [inst for inst in manager.instances if isinstance(inst, MovieDialog)]<br>
if not mdms:<br> raise AssertionError("No MD Movie Instances!")<br> return mdms[-1]<br>#####<br>runCommand("open movie:metafile; sel protein")<br>mdm = findMDMovie()<br>for x in range(1,1000+1,5):<br>
frame = x<br> mdm.LoadFrame(x)<br> runCommand("hbond selRestrict any relax true savefile hbond.%s" % x)<br> runCommand("wait")<br></div>runCommand("stop") # move onto next MD trajectory<br>
</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jan 14, 2014 at 11:06 AM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Inhee Park,<br>
A big limitation of the "coordset" command is that most formats of MD trajectory would have to be played through entirely using the graphical interface before it would work to use "coordset", as mentioned here:<br>
<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/coordset.html" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/coordset.html</a>><br>
<br>
In other words, probably the frame is not updating when you use coordset. Another thing that confuses me in your script is the use of "wait". I don't think it is needed. Assuming "coordset" is working, I'd use something like<br>
<br>
perframe myhb; coordset #0 1,1000,10<br>
~perframe<br>
<br>
Sorry about that limitation. Maybe someone else can suggest how to use python to play the trajectory.<br>
Best,<br>
Elaine<br>
----------<br>
Elaine C. Meng, Ph.D.<br>
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<div><div class="h5"><br>
<br>
On Jan 13, 2014, at 5:18 PM, inhee park <<a href="mailto:ipark.c@gmail.com">ipark.c@gmail.com</a>> wrote:<br>
<br>
> Dear Chimera developers and users:<br>
><br>
><br>
> I am using Chimera (production version 1.8 (build 38824) 2013-06-07 21:09:20 UTC).<br>
><br>
> I wanted to process MD trajectories to extract their detailed H-bond information.<br>
> To do so, I would like to process by script with Hbond.py file.<br>
> $ chimera --script Hbond.py<br>
><br>
> #---[Hbond.py]<br>
> from chimera import runCommand<br>
><br>
> # pr trajectory movie<br>
> with open("metafile", 'w') as m:<br>
> meta_input = """amber<br>
> mypdb.prmtop<br>
> pr_dry_netcdf.crd<br>
> """<br>
> print >> m, meta_input<br>
> #<br>
> runCommand("open movie:metafile")<br>
> runCommand("alias ^myhb hbond relax true savefile hbond.$1")<br>
> runCommand("perframe myhb")<br>
> runCommand("coordset #0 1,1000,10; wait 100")<br>
> #-----------------------<br>
><br>
> As I have many MD trajectories to be processed, I was looking for a commandline<br>
> in lieu of actually clicking buttons in MD movie GUI to play.<br>
> I thought that "coordset #0 1,10000,10; wait 100" may operate like automatically<br>
> starting the MD movie to load frame 1 to 1000 at every 100 interval.<br>
> But resulting saved H-bond files contain same information all about the first frame<br>
> (because graphically displayed was the first frame without any update<br>
> unless clicking the button in MD movie GUI).<br>
><br>
> Is there a way to achieve an automatic play via commandline?<br>
><br>
> Thank you,<br>
> inhee park<br>
</div></div>> _______________________________________________<br>
> Chimera-users mailing list<br>
> <a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br>
> <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>_________<br>inhee park
</div>