I'm dealing with a large system of proteins in water molecules of SPC model on
GMX 3.3.1. When doing position restraints, it always gives such error info as
follows and then stops:
t = 0.002 ps: Water molecule starting at atom 72750 can not be settled.
Check for bad contacts and/or reduce the timestep.
This error comes from src/mdlib/constr.c:
if (!bOK && bDumpOnError==1)
fprintf(stdlog,"\nt = %.3f ps: Water molecule starting at atom %d can
not be "
"settled.\nCheck for bad contacts and/or reduce the timestep.",
ir->init_t+step*ir->delta_t,owptr[error]+1);
I want to check this naughty guy out and test whether gmx can run without this
bad water molecule.The question is, there are ~500,000 atoms in the system and
the figure 72750 shows up 5 times. How can I know which one is the bad water
molecule? Is it the ?72750th (?=0,1,2,3,4) atom in the whole system or just in
water molecules?
//description of proteins
ATOM 1 N ALA 1 65.430 66.090 97.040 1.00 0.00
...
//description of water
ATOM 73347 OW SOL 8016 5.700 12.790 11.830 1.00 0.00
ATOM 73348 HW1 SOL 8016 4.860 12.540 11.310 1.00 0.00
ATOM 73349 HW2 SOL 8016 5.770 13.800 11.860 1.00 0.00
...
Any help would be appreciable.Thank you all in adcance.
Ying Ren
2007-11-16
_______________________________________________
gmx-users mailing list [email protected]
http://www.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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
