anchor.diff:       &   nfault,dfault,tfault,numfn,                                    ! split
        &   getshapn,bmatrixn,                                             ! bbar
        &   ierr,errstrng)                                                 ! errcode
  +cdebug        write(6,*) "After formf_ss (1):"
   c
           if(ierr.ne.izero) return
   c
  
positive.code: 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
  
negative.code: c
        call load(id,ibond,bond,dcur,deld,btot,histry,deltp,numnp,neq,
       & nhist,nstep,lastep,ierr,errstrng)
  cdebug      write(6,*) "After load:"
        if(ierr.ne.izero) return
  c
  c...compute current split node displacements
  