1      SUBROUTINE list_const(ncons,ntcons,consatm)
2
3      implicit none
4
5      include 'p_array.inc'
6      include 'p_input.inc'
7      include 'cm_atom.inc'
8      include 'cm_cons.inc'
9
10      integer i,j,k
11      integer ncons,ntcons,consatm
12
13      dimension consatm(mxcons,2)
14
15      ntcons=0
16
17      do i=1,natms
18
19       do j=1,ncons
20
21        if(atmtype(i).eq.consatm(j,1))then
22
23         ntcons=ntcons+1
24         if(ntcons.gt.mxcons2)then
25          write(output,*)'Increase mxcons2 to ',ntcons
26          stop
27         endif
28
29         k=consatm(j,2)-consatm(j,1)
30         conlist(ntcons,1)=i
31         conlist(ntcons,2)=i+k
32         conlist(ntcons,3)=j
33
34        endif
35
36       enddo
37
38      enddo
39
40      write(99,*) ntcons
41      do i=1,ntcons
42       write(99,*) conlist(i,1),conlist(i,2)
43
44      end do
45      return
46
47      END
48c $Id$
49