anchor.diff:              write(kw,2000)(head(j),j,j=1,nen)
               write(kw,2500)
             end if
  -          write(kw,3000) i,imat,ietype,ngauss,(ien(i),i=i1,i2)
  +          write(kw,3000) i,imat,ietype,ngauss,(ien(j),j=i1,i2)
           end do
           write(kw,4000)
           close(kw)
  
positive.code:       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
          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
  
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
  