2
c $Id: list_bonds.f 19707 2010-10-29 17:59:36Z d3y133 $
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'