VMD-L Mailing List
From: Anuradha Mittal (anuradha.mittal_at_gmail.com)
Date: Tue Apr 11 2006 - 13:47:06 CDT
- Next message: Pietro Amodeo: "Re: GLSL rendermode now OK on Linux with NVidia 8756 driver version"
- Previous message: John Stone: "Re: GLSL rendermode now OK on Linux with NVidia 8756 driver version"
- Next in thread: Peter Freddolino: "Re: RMSD"
- Reply: Peter Freddolino: "Re: RMSD"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]
Hi,
I am trying to find the RMSD of protein between each frame of a dcd file and
a reference pdb file. The PDB and DCD files have different number of water
molecules.
While running the script pasted below, it is giving an error of "incorrect
number of atoms".
proc myrmsd { frame } {
global ref sel all
$all move [measure fit $sel $ref]
puts "$frame: [measure rmsd $sel $ref]"
}
mol load psf prot.psf # for prot.dcd
set all [atomselect top all]
set ref [atomselect top "name CA" frame 0]
set sel [atomselect top "name CA"]
animate read pdb protein.pdb
bigdcd myrmsd prot.dcd
Is there any way of calulating rmsd from files with different number of
atoms?
Thanks
Anu
- Next message: Pietro Amodeo: "Re: GLSL rendermode now OK on Linux with NVidia 8756 driver version"
- Previous message: John Stone: "Re: GLSL rendermode now OK on Linux with NVidia 8756 driver version"
- Next in thread: Peter Freddolino: "Re: RMSD"
- Reply: Peter Freddolino: "Re: RMSD"
- Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]