[gmx-users] Software inconsistency error: update_coords called for velocity without VV integrator

[zhongjin]

Hi ,
 When I use g_membed in GMX4.5.3 Beta 3 to insert a protein into a DMPC bilayer,  I came across such an error:
Program g_membed, VERSION 4.5-beta3
Source code file: update.c, line: 1595
Software inconsistency error:
update_coords called for velocity without VV integrator
For more information and tips for troubleshooting, please check the GROMACS

 
g_membed -f gmem.tpr -p 1b8l.top -xyinit 0.5 -xyend 1 -nxy 1000 -ndiff 11 -c OK.pdb
 
Reading file gmem.tpr, VERSION 4.5-beta3 (single precision)
Reading file gmem.tpr, VERSION 4.5-beta3 (single precision)
Select a group to embed in the membrane:
Group     0 (         System) has 16956 elements
Group     1 (        Protein) has  3708 elements
Group     2 (      Protein-H) has  2904 elements
Group     3 (        C-alpha) has   388 elements
Group     4 (       Backbone) has  1164 elements
Group     5 (      MainChain) has  1556 elements
Group     6 (   MainChain+Cb) has  1904 elements
Group     7 (    MainChain+H) has  1940 elements
Group     8 (      SideChain) has  1768 elements
Group     9 (    SideChain-H) has  1348 elements
Group    10 (    Prot-Masses) has  3708 elements
Group    11 (    non-Protein) has 13248 elements
Group    12 (          Other) has 13248 elements
Group    13 (           DMPC) has 13248 elements
Select a group: 1
Selected 1: 'Protein'
Select a group to embed Protein into (e.g. the membrane):
Group     0 (         System) has 16956 elements
Group     1 (        Protein) has  3708 elements
Group     2 (      Protein-H) has  2904 elements
Group     3 (        C-alpha) has   388 elements
Group     4 (       Backbone) has  1164 elements
Group     5 (      MainChain) has  1556 elements
Group     6 (   MainChain+Cb) has  1904 elements
Group     7 (    MainChain+H) has  1940 elements
Group     8 (      SideChain) has  1768 elements
Group     9 (    SideChain-H) has  1348 elements
Group    10 (    Prot-Masses) has  3708 elements
Group    11 (    non-Protein) has 13248 elements
Group    12 (          Other) has 13248 elements
Group    13 (           DMPC) has 13248 elements
Select a group: 13
Selected 13: 'DMPC'
The estimated area of the protein in the membrane is 15.930 nm^2
There are 130 lipids in the membrane part that overlaps the protein.
The area per lipid is 0.6809 nm^2.
Maximum number of lipids that will be removed is 46.
Will resize the protein by a factor of 0.500 in the xy plane and 1.000 in the z direction.
This resizing will be done with respect to the geometrical center of all protein atoms
that span the membrane region, i.e. z between 1.299 and 6.362
Embedding piece 0 with center of geometry: 4.933612 4.933665 3.830500
Will remove 0 Protein_chain_A molecules
Will remove 0 Protein_chain_B molecules
Will remove 0 Protein_chain_C molecules
Will remove 0 Protein_chain_D molecules
Will remove 37 DMPC molecules
Back Off! I just backed up 1b8l.top to ./#1b8l.top.2#
Back Off! I just backed up traj.trr to ./#traj.trr.1#
Back Off! I just backed up traj.xtc to ./#traj.xtc.1#
Back Off! I just backed up ener.edr to ./#ener.edr.1#
starting mdrun 'God Rules Over Mankind, Animals, Cosmos and Such'
1000 steps,      2.0 ps.
-------------------------------------------------------
Program g_membed, VERSION 4.5-beta3
Source code file: update.c, line: 1595
Software inconsistency error:
update_coords called for velocity without VV integrator
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
-------------------------------------------------------
"Look at these, my work-strong arms" (P.J. Harvey)
 
How to solve such an problem? Thanks!
 
Zhongjin He


      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20100824/511ac82e/attachment.html>