[gmx-users] grompp bug (system charge)
David
spoel at xray.bmc.uu.se
Sun Apr 27 15:10:01 CEST 2003
On Sun, 2003-04-27 at 12:26, Christoph Freudenberger wrote:
> hi there,
>
> There seems to be a bug in the way grompp calculates the total
> charge of the system.
> For my box (9 kations/9 anions/some 10000 solvents) grompp prints
> out a warning:
> WARNING 4 [file "topol.top", line 22]:
> System has non-zero total charge: 7.999992e+00
>
> I dublechecked the charges in the resulting tpr using the following
> command-pipe:
> > gmxdump -s test.tpr |grep "atom\["|awk -F, '{print $5}'|awk -F= '{sum+=$2}END{print sum}'
>
> and get:
> -2.55351e-15
> which is neutral to me.
can you include the topology?
--
Groeten, David.
________________________________________________________________________
Dr. David van der Spoel, Dept. of Cell and Molecular Biology
Husargatan 3, Box 596, 75124 Uppsala, Sweden
phone: 46 18 471 4205 fax: 46 18 511 755
spoel at xray.bmc.uu.se spoel at gromacs.org http://xray.bmc.uu.se/~spoel
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
More information about the gromacs.org_gmx-users
mailing list