anchor.diff:       &         nfault,dfault,tfault,numfn,                              ! split
        &         getshapb,bmatrixb,                                       ! bbar
        &         ierr,errstrng)                                           ! errcode
  +cdebug              write(6,*) "After formf_ss (2):"
   c
                 if(ierr.ne.izero) return
   c
  
positive.code: c     nodes, and perform iterative solution.
  c
              if(lgdefp.eq.izero.and.intord.ne.ithree) then
                if(reform) call matinit_drv(
       &         alnz,ja,nnz,neq,                                         ! sparse
       &         x,d,iwink,wink,numnp,nwink,                              ! global
  
negative.code: c
  c...  define constants to control stepping in current group
  c
  cdebug          write(6,*) "Before const:"
            call const(maxstp,delt,alfa,maxit,ntdinit,lgdef,utol,
       &     ftol,etol,itmax,nintg,i,naxstp,nfirst,rtimdat,deltp,alfap,
       &     ntimdat,nstep,maxitp,ntdinitp,lgdefp,itmaxp,gtol)
  cdebug          write(6,*) "After const:"
  cdebug          write(6,*) nintg,i,naxstp,nfirst,deltp,alfap,maxitp,maxitcp,
  cdebug     &     lgdefp,itmaxp,(gtol(idb),idb=1,3)
            ltim=.true.
  