anchor.diff:  c
   c...  local variables
   c
  -      integer i,l,nshsize,ngssize
  +      integer nen,ngauss,nec,nee,i,l,nshsize,ngssize
         double precision g1,g2,rr,ss,tt,uu,drr,dss,dtt,duu
         double precision tetvol
   c
  
positive.code:             drr=eight*r(i)*max(ss,one)*max(tt,one)*max(uu,one)
              dss=eight*s(i)*max(rr,one)*max(tt,one)*max(uu,one)
              dtt=eight*t(i)*max(rr,one)*max(ss,one)*max(uu,one)
              duu=-eight*u(i)*max(rr,one)*max(ss,one)*max(tt,one)
              sh(1,i,l)=drr-duu
              sh(2,i,l)=dss-duu
              sh(3,i,l)=dtt-duu
  
negative.code:       end
  c
  c version
  c $Id: pquadtet.f,v 1.2 2004/07/07 19:53:01 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  