anchor.diff:  clater     &         s,stemp,                                                 ! stiff
   clater     &         state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,npar,     ! elemnt
   clater     &         ielno,iside,ihistry,pres,pdir,                           ! tractn
  -clater     &         prop,mhist,infmat,infmatmod,                             ! materl
  +clater     &         prop,mhist,infmat,infmatmod,tminmax,                     ! materl
   clater     &         gauss,sh,shj,infetype,                                   ! eltype
   clater     &         histry,rtimdat,ntimdat,nvisdat,                          ! timdat
   clater     &         rgiter,gcurr,gi,gprev,gtol,rmin,rmult,nsiter,            ! iterate
  
positive.code: c     requested, compute forces due to applied displacements and split
  c     nodes, and perform iterative solution.
  c
              if(lgdefp.eq.izero.and.ibbarp.eq.izero) then
                if(reform) call matinit_drv(
       &         alnz,ja,nnz,neq,                                         ! sparse
       &         x,d,iwink,wink,numnp,nwink,                              ! global
  
negative.code: c
                if(ierr.ne.izero) return
  c
  clater            else if(lgdefp.eq.ione.and.intord.ne.ithree) then
  clater              if(reform) call matinit_drv(
  clater     &         alnz,ja,nnz,neq,                                         ! sparse
  clater     &         x,d,iwink,wink,numnp,nwink,                              ! global
  