Dear Shay, dear Gromacians >For purposes of debugging this bug you can also use >the flag of g_mindist that gives out atom pairs and >see which atom pairs is it that gromacs thinks are > less the 1 angstrom. I used the "-o atm-pair.out" option to get the atom pairs with the minimal distance. For the all atom trajectory the correct atom pairs are selected. But for the CA trajectory the wrong atom pairs are selected. Instead of the pairs with the shortest distance a much further pair is selected. Although the output pair from "atm-pair.out" has a distance of more then 100 Angstrom the value in the file "mindist.xvg" is smaller then one Angstrom.
I think the "g_mindist" command has a bug. In the trajectory each residue is only represented by a single CA atom (to reduce storage footprint and increase speed), can it be that g_mindist requires that the atoms have to "connected" (no broken bonds)? But for calculating the minimal distance of two groups this should really not be required. Thanks Peter -- Sarah Kreuz, die DSDS-Siegerin der Herzen, mit ihrem eindrucksvollen Debütalbum "One Moment in Time". http://portal.gmx.net/de/go/musik -- gmx-users mailing list gmx-users@gromacs.org http://lists.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-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
