anchor.diff:       & prop,mhist,infmat,infmatmod,ismatmod,                            ! materl
        & gauss,sh,shj,infetype,                                           ! eltype
        & histry,rtimdat,ntimdat,nvisdat,maxstp,delt,alfa,maxit,ntdinit,   ! timdat
  -     & lgdef,ibbar,utol,ftol,etol,itmax,                                ! timdat
  +     & lgdef,utol,ftol,etol,itmax,                                      ! timdat
        & rgiter,gcurr,gi,gprev,gtol,rmin,rmult,nsiter,                    ! iterate
        & skew,                                                            ! skew
        & ncodat,nunits,nprint,istatout,                                   ! ioinfo
  
positive.code: c
          if(ierr.ne.izero) return
  c
  clater      else if(lgdefp.eq.ione.and.ibbarp.eq.izero) then
  clater        call matinit_drv(
  clater     &   alnz,ja,nnz,neq,                                               ! sparse
  clater     &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
negative.code: c******  Another thing to consider is that I will need to pass in some extra
  c******  info that isn't currently needed for the small strain case.
  c
        if(lgdefp.eq.izero.and.intord.ne.ithree) then
          call matinit_drv(
       &   alnz,ja,nnz,neq,                                               ! sparse
       &   x,d,iwink,wink,numnp,nwink,                                    ! global
  