anchor.diff:  c
           if(ierr.ne.izero) return
   c
  -      else if(lgdefp.eq.0.and.ibbarp.eq.1) then
  +      else if(lgdefp.eq.izero.and.ibbarp.eq.ione) then
           call matinit_drv(
        &   alnz,ja,nnz,neq,                                               ! sparse
        &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
positive.code: c
  c...add differential forces across internal free interfaces
  c
        if(numdif.ne.0) call loadx(btot,diforc,histry,idx,idhist,neq,
       & numnp,nhist,nstep,lastep,ierr,errstrng)
        if(ierr.ne.izero) return
        ii=1
  c
  c...  initialize elastic material matrices and stiffness matrix, 
  c     compute forces due to applied displacements and split nodes,
  
negative.code: c
  c...  local variables
  c
        integer igroup,naxstp,nfirst
        double precision time,tminmax
        logical*4 fulout,skc
  c
  