anchor.diff:              do l=1,ngauss
                 call dcopy(nstr,dstate(1,indstateg),ione,
        &         state(1,indstateg),ione)
  -              call daxpy(nstr,sub,state(1,indstateg+nstr),ione,
  -     &         dstate(1,indstateg+nstr),ione)
  -              call daxpy(nstr,one,dstate(1,indstateg+nstr),ione,
  -     &         state(1,indstateg+nstr),ione)
  +              call daxpy(nstr,sub,state(1,indstateg+ione),ione,
  +     &         dstate(1,indstateg+ione),ione)
  +              call daxpy(nstr,one,dstate(1,indstateg+ione),ione,
  +     &         state(1,indstateg+ione),ione)
  +              indstateg=indstateg+nstate
               end do
  -            indstateg=indstateg+incstate
             end do
           else if(nstate.eq.ithree) then
             do ielg=matgpt,matgpt+nmatel-1
  
positive.code:       end
  c
  c version
  c $Id: update_state.f,v 1.4 2004/07/12 21:07:20 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:         matmodel=infmat(1,imat)
          nmatel=infmat(2,imat)
          nstate=infmatmod(2,matmodel)
  c
  c...  loop over elements in group, updating appropriate state variables
  c
  