anchor.diff:       &   dx,iwinkx,winkx,numslp,numsn,nwinkx,                           ! slip
        &   tfault,numfn,                                                  ! fault
        &   s,stemp,                                                       ! stiff
  -     &   state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz,       ! elemnt
  -     &   ndmatsz,numelt,nconsz,                                         ! elemnt
  +     &   state,dstate,state0,dmat,ien,lm,lmx,lmf,infiel,iddmat,nstatesz,! elemnt
  +     &   nstatesz0,ndmatsz,numelt,nconsz,nprestrflag,ipstrs,ipauto,     ! elemnt
  +     &   nstate0,                                                       ! elemnt
        &   prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,             ! materl
        &   gauss,sh,shj,infetype,                                         ! eltype
        &   histry,rtimdat,ntimdat,rgiter,nhist,lastep,                    ! timdat
  
positive.code: c
  clater        call iterate(
  clater     &   alnz,pcg,zcg,dprev,ja,                                         ! sparse
  clater     &   bextern,btraction,bgravity,bconcforce,bprestress,bintern,      ! force
  clater     &   bresid,bwork,dispvec,nforce,grav,                              ! force
  clater     &   x,d,deld,dcur,id,iwink,wink,nsysdat,                           ! global
  clater     &   dx,deldx,dxcur,idx,iwinkx,winkx,idslp,ipslp,                   ! slip
  clater     &   nfault,dfault,tfault,                                          ! fault
  clater     &   s,stemp,                                                       ! stiff
  clater     &   state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,npar,           ! elemnt
  clater     &   ielno,iside,ihistry,pres,pdir,                                 ! tractn
  clater     &   prop,mhist,infmat,infmatmod,tminmax,                           ! materl
  clater     &   gauss,sh,shj,infetype,                                         ! eltype
  
negative.code:       integer istatout(*)
        double precision alnz(*),pcg(*),zcg(*),dprev(*)
        double precision bextern(*),btraction(*),bgravity(*),bconcforce(*)
        double precision bintern(*),bresid(*),bwork(*),dispvec(*),grav(*)
        double precision x(*),d(*),deld(*),dcur(*),wink(*)
        double precision bond(*)
        double precision dx(*),deldx(*),dxcur(*),diforc(*),winkx(*)
        double precision fault(*),dfault(*),tfault(*)
        double precision s(*),stemp(*)
        double precision state(*),dstate(*),state0(*),dmat(*)
        double precision  pres(*),pdir(*)
        double precision prop(*)
        double precision gauss(*),sh(*),shj(*)
  