anchor.diff:  c
   c...compute the global displacement increment vector using a
   c   preconditioned conjugate gradients iterative solver.  Upon
  -c   return the vector gvec2 contains the displacements.
  +c   return the vector dispvec contains the displacements.
   c
           call pcginv(alnz,bresid,dispvec,bwork,pcg,zcg,dprev,rmin,rmult,
        &   gcurr,gprev,ja,nsiter,neq,nprevdflag,nnz,ndtot,idout,kto,kw,
  
positive.code: c
  c...  subroutine arguments
  c
        integer ierr
        integer ja(*)
        integer id(*),iwink(*)
        integer idx(*),iwinkx(*),idslp(*),ipslp(*)
  
negative.code: c...loop over iterations
  c
        do iter=1,itmaxp
          updats=iter.eq.1.or.lgdefp.ge.1
          skc=iter.gt.1.and.lgdefp.ge.1.and.(numslp.ne.0.and.
       &   (iskopt.eq.2.or.(iskopt.le.0.and.abs(iskopt).eq.nstep)))
          nittot=nittot+1
  