Hi again,
my topology is non-consecutively numbered. I was looking for a tool to
renumber it, but didn't find it yet. A hint will be really appreciated.
I was not sure either whether or not is it possible to use the xtc
generated with the original topology with the renumbered one. My main
doubt is how to properly generate the index file for g_order. If each
group contain an atom type in the chain, let's say all C21, another all
C22, etc., I get many warnings of the type:
WARNING: distance between atoms 6767 and 6751 > 0.3 nm (1.406577). Index
file might be corrupt.
If all atoms belonging to one chain are put together in a group, then
the xvg files are empty.
There should be some posting in the archive clarifying the issue, but
was not able to find it yet.
Thanks,
Felipe
On 02/20/2013 04:17 PM, Justin Lemkul wrote:
On 2/20/13 9:41 AM, Felipe Pineda, PhD wrote:
Hi,
is it possible to apply g_order to a trajectory of a lipid bilayer when the
carbon atoms of the acyl chains of the lipid molecules for which the SCD should
be calculated are not numbered consecutively, i.e. the atom numbers
corresponding to C21-C22-C23-... (sn-2 chain) are eg 23 22 19 ... and these
appear in the corresponding group of the index file?
Theoretically the bonded structure should be taken from the .tpr file, but you
should probably do a quick test to be sure (i.e. consecutively numbered topology
vs. nonconsecutive).
-Justin
--
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists
