Dear Gromacs users,
We used g_mindist analysis as follows:
g_mindist -f *.xtc -s *.tpr -o atom-pairs.out -od mindist.xvg -or
res_mindist.xvg
and surprisingly in the atom-pairs.out only several atoms in one
residue (from group 1) were at minimum distance from group 2
throughout the entire simulation. Checking this output manually showed
it to be incorrect.
res_mindist.xvg seems to hold the correct data.
When running the same analysis, only omitting the -or flag, we get a
*correct* atom-pairs.out.
g_mindist -f *.xtc -s *.tpr -o atom-pairs.out -od mindist.xvg
Comparing mindist.xvg from both scenarios shows no difference (files
are identical).
Is the -or flag supposed to affect the other output files of the
g_mindist function?
Supplementary details:
1. The same thing happened when g_mindist -s was supplied with either
.gro, .pdb or .tpr files.
2. Index-groups 1 and 2 consisted upon full molecules with no
ommisions (full protein and not only C-alpha).
3. We ruled out PBC by using -nopi
4. This is common with the following Gromacs versions: 3.3.3., 4.0.5, 4.0.7.
Regards,
-Shay
--
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