Hello,

I am having trouble running scripts which vary selections on different states of a PyMOL object. I've converted an MD trajectory into a pdb with 500 models that PyMOL loads as 500 states of 1 object. For example, I would then like to run the anglesbetweenhelices.py script (from PyMOL wiki) on each state separately. A simple script example:

   import anglebetweenhelices as abh
   cmd.load('output.pdb', object="cdp")
   for i in range(1,501):
        cmd.set('state',i,'cdp')
        cmd.select('hel1','/cdp///1-22')
        cmd.select('hel2','/cdp///25-37')
        x=abh.angle_between_helices('hel1','hel2', visualize=0, quiet=1)
        print x

In theory the select command has a state option but this doesn't seem to work (as is confirmed by http://www.mail-archive.com/pymol-users@lists.sourceforge.net/msg08224.html). I then tried using the set state command, but this also does not seem to affect the actual coordinates of the selection. The only thing that does work is if in the GUI I go to the desired state by pressing the "<" and ">" animation buttons. But is there a way to select the state from within a script?

Thanks!
Pawel
------------------------------------------------------------------------------
Get your SQL database under version control now!
Version control is standard for application code, but databases havent 
caught up. So what steps can you take to put your SQL databases under 
version control? Why should you start doing it? Read more to find out.
http://pubads.g.doubleclick.net/gampad/clk?id=49501711&iu=/4140/ostg.clktrk
_______________________________________________
PyMOL-users mailing list (PyMOL-users@lists.sourceforge.net)
Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users
Archives: http://www.mail-archive.com/pymol-users@lists.sourceforge.net

Reply via email to