anchor.diff:        call fill(state,zero,nstr*nstatesz)
   c***  Note:  This should be modified for prestresses.
         call fill(dstate,zero,nstr*nstatesz)
  -      if(numfn.ne.0) call fill(tfault,zero,numfn*ndof)
  +      if(numfn.ne.izero) call fill(tfault,zero,numfn*ndof)
   c
         write(kto,600)
   c*      call flush(kto)
         fulout=.true.
  -      ireform=1
  -      igroup=1
  -      nstep=0
  +      ireform=ione
  +      igroup=ione
  +      nstep=izero
         ntimdat(1)=nstep
  -      naxstp=0
  -      nittot=0
  +      naxstp=izero
  +      nittot=izero
         ntimdat(7)=nittot
  -      nrftot=0
  +      nrftot=izero
         ntimdat(8)=nrftot
  -      ndtot=0
  +      ndtot=izero
         ntimdat(9)=ndtot
         ntimdat(10)=ireform
         call const(maxstp,delt,alfa,maxit,ntdinit,lgdef,ibbar,utol,ftol,
        & etol,itmax,nintg,igroup,naxstp,nfirst,rtimdat,deltp,alfap,
        & ntimdat,nstep,maxitp,ntdinitp,lgdefp,ibbarp,itmaxp,gtol)
  -      if(skc) call skclear(idslp,skew,numsn,numnp)
  -      if(skc) call skcomp(x,d,skew,idslp,ipslp,ipstrs,numsn,numnp,nstep,
  -     & lgdefp,ierr,errstrng)
  -      if(ierr.ne.izero) return
  +      if(skc) then
  +        call skclear(idslp,skew,numsn,numnp)
  +        call skcomp(x,d,skew,idslp,ipslp,ipstrs,numsn,numnp,nstep,
  +     &   lgdefp,ierr,errstrng)
  +        if(ierr.ne.izero) return
  +      end if
   c
   c...transfer boundary conditions into global load vector btot(neq)
   c   and displacement increment vector deld(ndof,numnp)
  
positive.code: c
          if(ierr.ne.izero) return
  c
  clater      else if(lgdefp.eq.1.and.ibbarp.eq.0) then
  clater        call matinit_drv(
  clater     &   alnz,ja,nnz,neq,                                               ! sparse
  clater     &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
negative.code:       end
  c
  c version
  c $Id: elastc.f,v 1.4 2004/07/13 17:40:12 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  