anchor.diff:        if(idout.gt.ione) close(kw)
         close(kp)
   c
  + 650  format(//,"Reforming the stiffness matrix:",/)
    700  format('STEP ',i7)
    800  format(/," Total number of equilibrium iterations        = ",i7,/,
        &         " Total number of stiffness matrix reformations = ",i7,/,
  
positive.code:               reform=(mod(j,ntdinitp).eq.izero)
              end if
              reform=reform.or.ltim
              if(reform) ireform=ione
              ntimdat(9)=ireform
              fulout=.false.
              if(ntot.eq.iprint(indexx)) fulout=.true.
  c
  
negative.code:               if(ierr.ne.izero) return
  c
  clater            else if(lgdefp.eq.ione.and.intord.ne.ithree) then
  clater              if(reform) then
  clater                write(kto,650)
  clater                call matinit_drv(
  clater     &           alnz,ja,nnz,neq,                                       ! sparse
  clater     &           x,d,iwink,wink,numnp,nwink,                            ! global
  clater     &           dx,iwinkx,winkx,numslp,numsn,nwinkx,                   ! slip
  clater     &           tfault,numfn,                                          ! fault
  clater     &           s,stemp,                                               ! stiff
  clater     &           state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,        ! elemnt
  clater     &           nstatesz,ndmatsz,numelt,nconsz,                        ! elemnt
  clater     &           prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,     ! materl
  clater     &           gauss,sh,shj,infetype,                                 ! eltype
  clater     &           histry,rtimdat,ntimdat,nhist,lastep,td_matinit_cmp_ld, ! timdat
  clater     &           skew,numrot,                                           ! skew
  clater     &           getshapn,bmatrixn,                                     ! bbar
  clater     &           ierr,errstrng)                                         ! errcode
  clater              end if
  c
  clater              if(ierr.ne.izero) return
  c
  