anchor.diff:       &         getshapb,bmatrixb,gload_cmp_ss,td_strs_cmp_ss,           ! external
        &         td_strs_mat_cmp_ss,                                      ! external
        &         ierr,errstrng)                                           ! errcode
  +cdebug              write(6,*) "After iterate (2):"
   c
                 if(ierr.ne.izero) return
   c
  
positive.code:      &         ndmatsz,numelt,nconsz,                                   ! elemnt
       &         prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,       ! materl
       &         gauss,sh,shj,infetype,                                   ! eltype
       &         histry,rtimdat,ntimdat,nhist,lastep,td_matinit_cmp_ss,   ! timdat
       &         skew,numrot,                                             ! skew
       &         getshapb,bmatrixb,                                       ! bbar
       &         ierr,errstrng)                                           ! errcode
  c
                if(ierr.ne.izero) return
  c
  
negative.code:             if(fulout.and.idsk.eq.ione) write(kp) ntot
              if(fulout.and.idsk.eq.ione) write(kp) time
              write(kto,5000) time,ntot,lastep*ncycle
              call flush(kto)
  c*            call flush(kw)
  c*            call flush(kp)
  c
  c...  apply boundary conditions
  c
  cdebug            write(6,*) "Before load:"
              call load(id,ibond,bond,d,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
  