[gmx-users] gro does not match top
Shima Arasteh
shima_arasteh2001 at yahoo.com
Sun Mar 24 13:46:14 CET 2013
Hi,
I am solving a protein in water. This protein has a charge of +3.
After solvation there is 3272 water molecules.
Next, I use genion command to add NACL 1M and 3 extra CL ions to neutralize the system.
#genion -s ions.tpr -o cyclic_solv_ions.gro -p topol.top -nname CL -nn 3 -conc 1
After this, the top file shows:
Protein_chain_A 1
SOL 3147
NA 61
CL 61
CL 3
But when I check the cyclic_solv_ions.gro file, there are 60 NA and 60 CL ions! Therefore the top doesnt match with gro file and I get the fatal error.
What' s your suggestions? How does it come?
Thanks in advance.
Sincerely,
Shima
More information about the gromacs.org_gmx-users
mailing list