anchor.diff:  c
           call fill(dcur,zero,ndof*numnp)
           call fill(dxcur,zero,ndof*numnp)
  -        call disp(gvec2,dcur,id,numnp,neq)
  -        if(numslp.ne.0) call disp(gvec2,dxcur,idx,numnp,neq)
  +        call disp(dispvec,dcur,id,numnp,neq)
  +        if(numslp.ne.0) call disp(dispvec,dxcur,idx,numnp,neq)
   c
   c...rotate skewed coordinates to global system
   c
  
positive.code: c
  c...  included variable definitions
  c
        include "nsysdat_def.inc"
        include "npar_def.inc"
        include "rtimdat_def.inc"
  
negative.code:           write(kto,650)
            call stress_mat_drv(
       &     alnz,ja,nnz,                                                 ! sparse
       &     bintern,neq,                                                 ! force
       &     x,d,iwink,wink,numnp,nwink,                                  ! global
       &     dx,iwinkx,winkx,numslp,numsn,nwinkx,                         ! slip
       &     tfault,numfn,                                                ! fault
  