anchor.diff:        integer ind,iel,indien,ietype,indstate,inddmat
         integer ngauss,nen,nee,l,indstateg,inddmatg,ngtest
         integer ngaussdim,incstate
  +      double precision tmax
         double precision dl(60),xl(60),scur(162),ee(162),p(60),det(27)
   c
   c...  included variable definitions
  
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:       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)
  