anchor.diff:       &         getshapn,bmatrixn,gload_cmp_ss,td_strs_cmp_ss,           ! external
        &         td_strs_mat_cmp_ss,                                      ! external
        &         ierr,errstrng)                                           ! errcode
  +cdebug              write(6,*) "After iterate (1):"
   c
                 if(ierr.ne.izero) return
   c
               else if(lgdefp.eq.izero.and.intord.eq.ithree) then
  +cdebug              write(6,*) "Before matinit_drv (2):"
  +cdebug              write(6,*) "reform:",reform
                 if(reform) call matinit_drv(
        &         alnz,ja,nnz,neq,                                         ! sparse
        &         x,d,iwink,wink,numnp,nwink,                              ! global
  
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:      &         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
  