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