Fixed NBONDS in psf writers for bead models

beadModel.py

@@ -902,7 +902,7 @@ class beadModel():
 
             ## Write out bonds
             bonds = self._getBonds()
-            fh.write("{:>8d} !NBOND\n".format(len(bonds)-1))
+            fh.write("{:>8d} !NBOND\n".format(len(bonds)))
             counter = 0
             for n,seps in bonds:
                 fh.write( "{:d} {:d} ".format(n[0].idx+1,n[1].idx+1) )

beadModelTwist.py

@@ -1000,7 +1000,7 @@ class beadModelTwist():
 
             ## Write out bonds
             bonds = self.bonds
-            fh.write("{:>8d} !NBOND\n".format(len(bonds)-1))
+            fh.write("{:>8d} !NBOND\n".format(len(bonds)))
             counter = 0
             for n1,n2,pot in bonds:
                 fh.write( "{:d} {:d} ".format(n1.idx+1,n2.idx+1) )