anchor.diff:  c
         indmat(1)=ione
         do imat=2,numat
  -        indmat(imat)=indmat(imat-1)+infmat(2,imat)
  +        indmat(imat)=indmat(imat-1)+infmat(2,imat-1)
         end do
  +cdebug      write(6,*) "indmat:",(indmat(idb),idb=1,numat)
   c
   c...  read connectivity, material number, and infinite element info
   c
  
positive.code:       end
  c
  c version
  c $Id: read_connect.f,v 1.3 2004/08/02 21:25:14 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:           if(imgrp(imat).eq.izero) then
              infiel(6,i)=ndmatsz+ione
              ndmatsz=ndmatsz+ngaussmax
              imgrp(imat)=infiel(6,i)
            else
              infiel(6,i)=imgrp(imat)
  cdebug            infiel(6,i)=infiel(6,infiel(4,i-1))
            end if
          else
            infiel(6,i)=ndmatsz+ione
            ndmatsz=ndmatsz+ngauss
          end if
  cdebug        write(6,*) i,(infiel(jdb,i),jdb=1,6)
        end do
        close(kr)
  cdebug      write(6,*) "infiel:"
  cdebug      do idb=1,numelt
  cdebug        write(6,*) idb,(infiel(jdb,idb),jdb=1,6)
  cdebug      end do
  c
  c...  output plot info, if desired
  c
  