anchor.diff:              matmodel=infmat(1,imat)
               nstate=infmatmod(2,matmodel)
               indstate=infiel(5,iel)+matmodpt(i,matmodel)
  -            incgauss=nstate*nstr
               do l=1,ngauss
  -              indstateg=indstate+(l-1)*incgauss
  +              indstateg=indstate+(l-1)*nstate
                 if(ismatmod(i,matmodel).eq.izero) then
                   call fill(stmp,zero,nstr)
                 else
  
positive.code:               if(idout.gt.1) then
                  nout=nout+1
                  if(nout.eq.1.or.mod(nout,m50).eq.izero) then
                    write(kw,2000) (labels(j),j=1,nstr)
                    write(kw,*) ' '
                  end if
                  write(kw,3000) iel,l,(stmp(j),j=1,nstr)
  
negative.code: c
        integer matmodpt(nstatesmax,nmatmodmax)
        integer m50,i,iel,imat,ietype,ngauss,matmodel,nstate,indstate,l
        integer indstateg,nout,j,iout
        double precision stmp(6),tmult
        character statedescr*21
  cdebug      integer idb,jdb
  c
  c...  included variable definitions
  c
  