anchor.diff:  c
                 call iterate(
        &         alnz,pcg,zcg,dprev,ja,                                   ! sparse
  -     &         bextern,btraction,bgravity,bconcforce,bprestress,bintern,! force
  -     &         bresid,bwork,dispvec,nforce,grav,                        ! force
  +     &         bextern,btraction,bgravity,bconcforce,bintern,bresid,    ! force
  +     &         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
  -     &         state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,npar,     ! elemnt
  +     &         state,dstate,state0,dmat,ien,lm,lmx,lmf,infiel,iddmat,   ! elemnt
  +     &         npar,                                                    ! elemnt
        &         ielno,iside,ihistry,pres,pdir,                           ! tractn
        &         prop,mhist,infmat,infmatmod,tminmax,                     ! materl
        &         gauss,sh,shj,infetype,                                   ! eltype
  
positive.code: clater     &           dx,iwinkx,winkx,numslp,numsn,nwinkx,                   ! slip
  clater     &           tfault,numfn,                                          ! fault
  clater     &           s,stemp,                                               ! stiff
  clater     &           state,dstate,dmat,ien,lm,lmx,lmf,infiel,iddmat,        ! elemnt
  clater     &           nstatesz,ndmatsz,numelt,nconsz,                        ! elemnt
  clater     &           prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,     ! materl
  clater     &           gauss,sh,shj,infetype,                                 ! eltype
  clater     &           histry,rtimdat,ntimdat,nhist,lastep,td_matinit_cmp_ld, ! timdat
  
negative.code:               reform=(mod(j,ntdinitp).eq.izero)
              end if
              reform=reform.or.ltim
              if(reform) ireform=ione
              ntimdat(9)=ireform
  cdebug            write(6,*) "ntdinitp,reform,ireform,ltim:",
  cdebug     &       ntdinitp,reform,ireform,ltim
              fulout=.false.
              if(ntot.eq.iprint(indexx)) fulout=.true.
  c
  