anchor.diff:          gcurr(i)=10.0d0*gtol(i)
           gprev(i)=10.0d0*gtol(i)
         end do
  +cdebug      write(6,*) "reform,ireform:",reform,ireform
   c
   c...loop over iterations
   c
  
positive.code: c
        subroutine iterate(
       & alnz,pcg,zcg,dprev,ja,                                           ! sparse
       & bextern,btraction,bgravity,bconcforce,bprestress,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
       & ielno,iside,ihistry,pres,pdir,                                   ! tractn
       & prop,mhist,infmat,infmatmod,tminmax,                             ! materl
       & gauss,sh,shj,infetype,                                           ! eltype
  
negative.code:      &     x,d,numnp,                                                   ! global
       &     dx,numslp,                                                   ! slip
       &     tfault,numfn,                                                ! fault
       &     state,dstate,state0,dmat,ien,lm,lmx,lmf,infiel,iddmat,       ! elemnt
       &     nstatesz,nstatesz0,ndmatsz,numelt,nconsz,nprestrflag,ipstrs, ! elemnt
       &     ipauto,nstate0,                                              ! elemnt
       &     prop,mhist,infmat,infmatmod,numat,npropsz,tminmax,           ! materl
       &     gauss,sh,shj,infetype,                                       ! eltype
       &     histry,rtimdat,rgiter,ntimdat,nhist,lastep,stress_cmp,       ! timdat
  