Commit 67cd9e18 authored by Jason Swails's avatar Jason Swails
Browse files

Add TIP4PEw name.

parent e2c80f2e
......@@ -720,7 +720,7 @@ def readAmberSystem(prmtop_filename=None, prmtop_loader=None, shake=None, gbmode
system.addParticle(mass)
# Add constraints.
isWater = [prmtop.getResidueLabel(i) == ('WAT', 'TP4', 'TP5') for i in range(prmtop.getNumAtoms())]
isWater = [prmtop.getResidueLabel(i) == ('WAT', 'TP4', 'TP5', 'T4E') for i in range(prmtop.getNumAtoms())]
if shake in ('h-bonds', 'all-bonds', 'h-angles'):
for (iAtom, jAtom, k, rMin) in prmtop.getBondsWithH():
system.addConstraint(iAtom, jAtom, rMin)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment