anchor.diff:  c
   c...  local variables
   c
  -      integer matgpt,imat,matmodel,nmatel,nstate,incstate,ielg,iel
  +      integer matgpt,imat,matmodel,nmatel,nstate,ielg,iel
         integer ietype,indstate,indstateg,ngauss,l
         double precision sub
   c
  
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:             do l=1,ngauss
                call dcopy(nstr,dstate(1,indstateg),ione,
       &         state(1,indstateg),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
            end do
          else if(nstate.eq.ithree) then
            do ielg=matgpt,matgpt+nmatel-1
  