Jones de Andrade wrote:

 >From .gro or .pdb file in the standard way.

     > number of molecules types,
    No
     > number of molecules of each type,
    No
     > number of atoms of each molecule type,
    No


I think those three can become yes in a really easy way: some kind of "state machine" would easily do the trick and "count" them properly.

Yes, but now you have to infer molecule changes from residue name changes, which is doable, but messy.

     > coordinates,
    Yes


 From .xtc or .pdb file.

    charges
    Yes, if you use the right editconf option
    and masses on a .pdb file, correct?
    No.


Only modifying editconf to print those in a adapted .pdb file, throw a file conversion and much disk space "wasted". Or I can learn exactly from editconf how to read it from .tpr or .top (I guess it takes from these files the atomic charges when needed) and adapt into my codes.

The point of editconf is to feed your system *one* pdb structure with atom names, charges and/or whatever - it would be easy to modify editconf to write masses in one column and charges in another. You rely on the user to then supply a .xtc file that corresponds to that .pdb, in the same way that most of the GROMACS analysis tools accept an .xtc and various subsets of the possible structure files. You don't convert your .xtc to a .pdb trajectory, because that would be immensely wasteful, and the static information of the names, masses and charges doesn't vary with time under the MD approximation.

    If you use the editconf suggestion to get charges, you don't need to
    read a .tpr file.


But then or I get locked into the problem of generality of having too much disk usage. But that idea will be also usefull to learn how to read the masses and charges from .tpr files.

Easier is to hack editconf as described above.

Mark
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to