anchor.diff:       & s,stemp,                                                         ! stiff
        & state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz,         ! elemnt
        & ndmatsz,numelt,nconsz,                                           ! elemnt
  -     & prop,infmat,nprop,matgpt,elas_matinit,td_matinit,matchg,         ! materl
  +     & prop,infmat,nprop,matgpt,elas_matinit,td_matinit,matchg,tminmax, ! materl
        & gauss,sh,shj,infetype,                                           ! eltype
        & rtimdat,ntimdat,                                                 ! timdat
        & skew,numrot,                                                     ! skew
  
positive.code: c
        nmatel=infmat(2)
        imatvar=infmat(4)
  c
  c...  compute d-matrix for first element in the group.  If imatvar is
  c     zero, there is only one d-matrix for the group, and that is all
  
negative.code:       double precision dx(ndof,numnp),tfault(ndof,numfn)
        double precision s(neemax*neemax),stemp(neemax*neemax)
        double precision state(nstr,nstatesz),dstate(nstr,nstatesz)
        double precision dmat(nddmat,ndmatsz),prop(nprop),tminmax
        double precision gauss(nsd+1,ngaussmax,netypes)
        double precision sh(nsd+1,nenmax,ngaussmax,netypes)
        double precision shj(nsd+1,nenmax,ngaussmax,netypes)
  