anchor.diff:        fac=e/(pr2*pr3)
         dd=pr1*fac
         od=pr*fac
  -      ss=pr3*fac
  +      ss=half*pr3*fac
         do i=1,3
           dmat(iddmat(i,i))=dd
           dmat(iddmat(i+3,i+3))=ss
  
positive.code:       end
  c
  c
        subroutine update_state_1(state,dstate,nstate)
  c
  c...  routine to update state variables at the end of a time step.
  
negative.code: c       
  
  c version
  c $Id: mat_1.f,v 1.18 2005/03/24 02:08:57 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  