anchor.diff:  c
         integer matmodpt(nstatesmax,nmatmodmax)
         integer m50,i,iel,imat,ietype,ngauss,matmodel,nstate,indstate,l
  -      integer indstateg,nout,j,incgauss
  +      integer indstateg,nout,j,incgauss,iout
         double precision stmp(6),tmult
   c
   c...  included variable definitions
  
positive.code: c
  c...  subroutine arguments
  c
        integer nstatesz,numelt,numat,idout,idsk,kw,kp
        integer infiel(6,numelt),infmat(3,numat),infmatmod(5,nmatmodmax)
        integer ismatmod(nstatesmax,nmatmodmax),infetype(4,netypes)
        integer istatout(2,nstatesmax)
        double precision delt
        double precision state(nstr,nstatesz),dstate(nstr,nstatesz)
  c
  
negative.code: c
  c...  output increments/rates, if desired.
  c
          if(istatout(2,i,iout).ne.izero) then
            tmult=one
            if(istatout(2,i,iout).eq.ione.and.delt.gt.zero) tmult=one/delt
            nout=izero
            do iel=1,numelt
              imat=infiel(2,iel)
  