anchor.diff:        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)
  +      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)
  
positive.code:         do l=1,ngauss
            call td_strs_mat(state(1,indstateg),dstate(1,indstateg),
       &     ee(nstr*(l-1)),dmat(1,inddmatg),prop,rtimdat,rgiter,ntimdat,
       &     nstate,matchg)
            call dcopy(nstr,dstate(1,indstateg),ione,scur(nstr*(l-1)),
       &     ione)
            indstateg=indstateg+incstate
  
negative.code:       include "rgiter_def.inc"
        include "ntimdat_def.inc"
  c
  cdebug      write(6,*) "Hello from td_strs_mat_cmp_ss_f!"
  c
        ngtest=izero
        if(imatvar.eq.izero) ngtest=ngaussmax
  