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,tminmax, ! materl
  +     & prop,nprop,matgpt,imatvar,nmatel,elas_matinit,td_matinit,matchg, ! materl
  +     & tminmax,                                                         ! materl
        & gauss,sh,shj,infetype,                                           ! eltype
        & rtimdat,ntimdat,                                                 ! timdat
        & skew,numrot,                                                     ! skew
  
positive.code: c...  subroutine arguments
  c
        integer nnz,neq,numnp,numslp,numsn,numfn,nstatesz,ndmatsz,numelt
        integer nconsz,nprop,matgpt,numrot,ierr
        logical matchg
        integer ja(nnz),ien(nconsz),lm(ndof,nconsz),lmx(ndof,nconsz)
        integer lmf(nconsz),infiel(6,numelt),iddmat(nstr,nstr),infmat(6)
        integer infetype(4,netypes)
        character errstrng*(*)
        double precision alnz(nnz),x(nsd,numnp),d(ndof,numnp)
  
negative.code: c
  cdebug      write(6,*) "Hello from elas_matinit_cmp_ss_f!"
  c
        tminmax=big
  c
  c...  compute d-matrix for first element in the group.  If imatvar is
  