[gmx-users] Numbering in xtc and in conf.gro. Does it coincide?

David van der Spoel spoel at xray.bmc.uu.se
Tue Jan 20 15:54:23 CET 2009


Vitaly Chaban wrote:
>> Source code file: matio.c, line: 533
>>
>> Fatal error:
>> Lo: 0.000000, Mid: 1.000000, Hi: 1.000000
>>
>>
>> Although nclust.xvg, maxclust.xvg, avclust.xvg are created and seem
>> adequate.
>>
>> However what could the above error indicate?
> 
>> that all atoms belong to the same cluster.
> 
> Yes, it's true in my case.
> 
> 
> And one more question about index file. If xtc-qrps is used and only
> the groups of interest are output in xtc file, should the atoms in the
> index file be renumbered or not?
> 
> For example, if in conf.gro we have atoms of interest with numbers 4,5,6
> and didn't write atoms 1,2,3 to the xtc file, should we then refer to the
> wanted atoms as 4,5,6 (as it was initially) or 1,2,3 (due to the
> shift...)?
> 
> In other words, is information about the initial (in conf.gro)
> numbering of molecules present in the xtc trajectory file? Looking at
> the xtc using gmxdump I see the numbering always start with zero in
> spite of the xtc-qrps arguments.
> 
All files number from 0 to natom-1. You have to do the remapping.
> 
> Thank you.
> Vitaly
> 
> 
> _______________________________________________
> gmx-users mailing list    gmx-users at 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 gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php


-- 
David van der Spoel, Ph.D., Professor of Biology
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:	+46184714205. Fax: +4618511755.
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://folding.bmc.uu.se



More information about the gromacs.org_gmx-users mailing list