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