anchor.diff:        open(kr,file=ifile,status="old",err=20)
         call pskip(kr)
         do i=1,numelt
  -        if(i.ne.ione) infiel(1,i)=infiel(1,i-1)+neni(ietypei)-1
  +        if(i.ne.ione) infiel(1,i)=infiel(1,i-1)+neni(ietypei)
           read(kr,*,end=30,err=30) n,ietypei,imat,inf,
        &   (ien(j),j=infiel(1,i),infiel(1,i)+neni(ietypei)-1)
           infiel(2,i)=imat
  
positive.code: c
  1000  format(30i7)
  2000  format(1x,///,
       x' e l e m e n t  d a t a ',//,5x,
       x' element    material  element  number',20(a4,i2,4x))
  2500  format(5x,'  number      type    gauss',/,
       x       5x,'                     points',/)
  3000  format(6x,i7,25(5x,i7))
  4000  format(//)
  c
        end
  c
  c version
  c $Id: read_connect.f,v 1.2 2004/07/09 20:22:41 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c     total size of the state variable and material matrices.
  c
          infiel(5,i)=nstatesz+ione
          ngauss=infetype(1,infiel(3,i))
          nstate=infmatmod(2,matmod)
          imatvar=infmatmod(4,matmod)
          nstatesz=nstatesz+nstate*ngauss
          if(imatvar.eq.izero) then
            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
  