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: c
        indmat(1)=ione
        do imat=2,numat
          indmat(imat)=indmat(imat-1)+infmat(2,imat)
        end do
  c
  c...  read connectivity, material number, and infinite element info
  c
  
negative.code:           errstrng="read_connect"
            return
          end if
  cdebug        write(6,*) "i,imat,indmat(imat):",i,imat,indmat(imat)
          infiel(4,indmat(imat))=i
          indmat(imat)=indmat(imat)+ione
  c
  