anchor.diff:                if(ierr.ne.izero) return
   c
                 call iterate(
  -     &         alnz,pcg,zcg,ja,                                         ! sparse
  -     &         b,btot,bres,pvec,gvec1,gvec2,grav,                       ! force
  -     &         x,d,deld,dprev,dcur,id,iwink,wink,nsysdat,               ! global
  +     &         alnz,pcg,zcg,dprev,ja,                                   ! sparse
  +     &         bextern,btraction,bgravity,bconcforce,bprestress,bintern,! force
  +     &         bresid,bwork,dispvec,nforce,grav,                        ! force
  +     &         x,d,deld,dcur,id,iwink,wink,nsysdat,                     ! global
        &         dx,deldx,dxcur,idx,iwinkx,winkx,idslp,ipslp,             ! slip
        &         nfault,dfault,tfault,                                    ! fault
        &         s,stemp,                                                 ! stiff
  
positive.code: c
  c...  included dimension and type statements
  c
        include "nsysdat_dim.inc"
        include "npar_dim.inc"
        include "ntimdat_dim.inc"
  
negative.code: c
  c...  subroutine arguments
  c
        integer ierr
        integer ja(*)
        integer id(*),iwink(*)
        integer ibond(*)
        integer idx(*),iwinkx(*),idslp(*),ipslp(*),idhist(*)
        integer nfault(*)
        integer ien(*),lm(*),lmx(*),lmf(*),infiel(*),iddmat(*)
        integer ielno(*),iside(*),ihistry(*)
        integer mhist(*),infmat(*),infmatmod(*),ismatmod(*)
        integer infetype(*)
        integer maxstp(*),maxit(*),ntdinit(*),lgdef(*),itmax(*)
        integer iprint(*)
        integer istatout(*)
        double precision alnz(*),pcg(*),zcg(*),dprev(*)
        double precision bextern(*),btraction(*),bgravity(*),bconcforce(*)
        double precision bprestress(*),bintern(*),bresid(*),bwork(*)
        double precision dispvec(*),grav(*)
        double precision x(*),d(*),deld(*),dcur(*),wink(*)
        double precision bond(*)
        double precision dx(*),deldx(*),dxcur(*),diforc(*),winkx(*)
        double precision fault(*),dfault(*),tfault(*)
  