anchor.diff:  c
         integer imat,npage,matmod,i,n,j,inf,ietypei,ngauss,nstate,imatvar
         integer ietype,nen,i1,i2
  +cdebug      integer idb,jdb
   c
   cdebug      write(6,*) "Hello from read_connect_f!"
   cdebug      write(6,*) nconsz,numelt,numnp,numat,nstatesz,ndmatsz,kr,kw,kp
  -cdebug      write(6,*) (neni(i),i=1,netypesi)
  -cdebug      do i=1,netypes
  -cdebug        write(6,*) (infetype(j,i),j=1,4)
  -cdebug      end do
  -cdebug      do i=1,numat
  -cdebug        write(6,*) (infmat(j,i),j=1,3)
  -cdebug      end do
  -cdebug      do i=1,nmatmodmax
  -cdebug        write(6,*) (infmatmod(j,i),j=1,5)
  +cdebug      write(6,*) (neni(idb),idb=1,netypesi)
  +cdebug      write(6,*) "infetype:"
  +cdebug      do idb=1,netypes
  +cdebug        write(6,*) (infetype(jdb,idb),jdb=1,4)
  +cdebug      end do
  +cdebug      write(6,*) "infmat:"
  +cdebug      do idb=1,numat
  +cdebug        write(6,*) (infmat(jdb,idb),jdb=1,3)
  +cdebug      end do
  +cdebug      write(6,*) "infmatmod:"
  +cdebug      do idb=1,nmatmodmax
  +cdebug        write(6,*) (infmatmod(jdb,idb),jdb=1,5)
   cdebug      end do
         call ifill(ien,izero,nconsz)
         call ifill(indmat,izero,numat)
  
positive.code:           if(imgrp(imat).eq.izero) then
              infiel(6,i)=ndmatsz+ione
              ndmatsz=ndmatsz+ngaussmax
            else
              infiel(6,i)=infiel(6,infiel(4,i-1))
            end if
          else
            infiel(6,i)=ndmatsz+ione
            ndmatsz=ndmatsz+ngauss
          end if
          imgrp(imat)=imgrp(imat)+ione
  cdebug        write(6,*) i,(infiel(j,i),j=1,6)
        end do
        close(kr)
  c
  c...  output plot info, if desired
  c
  
negative.code:       end
  c
  c version
  c $Id: read_connect.f,v 1.4 2004/08/12 02:24:21 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  