2
c $Id: list_bonds.f,v 1.1 2008-04-18 17:48:12 marat Exp $
5
SUBROUTINE list_bonds(nbond,ntbond,bondatm)
15
integer nbond,ntbond,bondatm
17
dimension bondatm(mxbond,2)
25
if(atmtype(i).eq.bondatm(j,1))then
28
if(ntbond.gt.mxbond2)then
29
write(output,*)'Increase mxbond2 to ',ntbond
33
k=bondatm(j,2)-bondatm(j,1)
34
if(typmol(atmtype(i)).ne.typmol(atmtype(i+k)))then
35
write(output,*)'Bonded atoms must be from the same molecule'