[gmx-users] Had anybody modied the source code of gromacs

[聂雪川]

Hello gmx-users,
    I just have one question about  modifing the source code of gromacs4.0.7.For I want to control the motion of one certain atom (e.g. atom 200).I added some codes like "if(n==my_atom)xprime[n][d]   = x[n][d]+...;"  in the functionstatic void do_update_md() of the file gromacs-4.0.7\src\mdlib\update.c.But when I run it Parallel (mpi run),I find that the atom serial number n isn't the real atom serial number n on each CPU.So the code "if(n==my_atom)" is futile.How can I achieve it.
                                                              Thank your.
                                                                  chuan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100419/d5104b69/attachment.html>