anchor.diff:              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
  -     &         dx,iwinkx,winkx,numslp,numsn,nwinkx,                     ! slip
  -     &         tfault,numfn,                                            ! fault
  -     &         s,stemp,                                                 ! stiff
  -     &         state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz, ! elemnt
  -     &         ndmatsz,numelt,nconsz,                                   ! elemnt
  -     &         prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,       ! materl
  -     &         gauss,sh,shj,infetype,                                   ! eltype
  -     &         histry,rtimdat,ntimdat,rgiter,nhist,lastep,              ! timdat
  -     &         td_matinit_cmp_ss,                                       ! timdat
  -     &         skew,numrot,                                             ! skew
  -     &         getshapb,bmatrixb,                                       ! bbar
  -     &         ierr,errstrng)                                           ! errcode
  +              if(reform) then
  +                write(kto,650)
  +                call matinit_drv(
  +     &           alnz,ja,nnz,neq,                                       ! sparse
  +     &           x,d,iwink,wink,numnp,nwink,                            ! global
  +     &           dx,iwinkx,winkx,numslp,numsn,nwinkx,                   ! slip
  +     &           tfault,numfn,                                          ! fault
  +     &           s,stemp,                                               ! stiff
  +     &           state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,        ! elemnt
  +     &           nstatesz,ndmatsz,numelt,nconsz,                        ! elemnt
  +     &           prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,     ! materl
  +     &           gauss,sh,shj,infetype,                                 ! eltype
  +     &           histry,rtimdat,ntimdat,rgiter,nhist,lastep,            ! timdat
  +     &           td_matinit_cmp_ss,                                     ! timdat
  +     &           skew,numrot,                                           ! skew
  +     &           getshapb,bmatrixb,                                     ! bbar
  +     &           ierr,errstrng)                                         ! errcode
  +              end if
   cdebug              write(6,*) "After matinit_drv (2):"
   c
                 if(ierr.ne.izero) return
  
positive.code:             if(lgdefp.eq.izero.and.intord.ne.ithree) then
  cdebug              write(6,*) "Before matinit_drv (1):"
  cdebug              write(6,*) "reform:",reform
                if(reform) call matinit_drv(
       &         alnz,ja,nnz,neq,                                         ! sparse
       &         x,d,iwink,wink,numnp,nwink,                              ! global
       &         dx,iwinkx,winkx,numslp,numsn,nwinkx,                     ! slip
       &         tfault,numfn,                                            ! fault
       &         s,stemp,                                                 ! stiff
       &         state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz, ! elemnt
       &         ndmatsz,numelt,nconsz,                                   ! elemnt
       &         prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,       ! materl
       &         gauss,sh,shj,infetype,                                   ! eltype
       &         histry,rtimdat,ntimdat,rgiter,nhist,lastep,              ! timdat
       &         td_matinit_cmp_ss,                                       ! timdat
       &         skew,numrot,                                             ! skew
       &         getshapn,bmatrixn,                                       ! bbar
       &         ierr,errstrng)                                           ! errcode
  cdebug              write(6,*) "After matinit_drv (1):"
  c
                if(ierr.ne.izero) return
  
negative.code:       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,/,
  