anchor.diff:            gauss(1,l)=g1
             gauss(2,l)=g1
             gauss(3,l)=g1
  -          if(l.ne.ngauss) gauss(l,l)=g2
  +          if(l.ne.ione) gauss(l-ione,l)=g2
             gauss(4,l)=fourth*tetvol
           end do
         end if
  
positive.code:         rr=gauss(1,l)
          ss=gauss(2,l)
          tt=gauss(3,l)
          uu=g1
          if(l.eq.4) uu=g2
          if(intord.eq.itwo) uu=fourth
  cdebug        write(6,*) "l:",l
  cdebug        write(6,*) "rr,ss,tt,uu:",rr,ss,tt,uu
          sh(4,1,l)=rr*(two*rr-one)
          sh(1,1,l)=four*rr-one
          sh(2,1,l)=zero
          sh(3,1,l)=zero
          sh(4,2,l)=ss*(two*ss-one)
          sh(1,2,l)=zero
          sh(2,2,l)=four*ss-one
          sh(3,2,l)=zero
          sh(4,3,l)=tt*(two*tt-one)
          sh(1,3,l)=zero
          sh(2,3,l)=zero
          sh(3,3,l)=four*tt-one
          sh(4,4,l)=uu*(two*uu-one)
          sh(1,4,l)=-(four*uu-one)
          sh(2,4,l)=-(four*uu-one)
          sh(3,4,l)=-(four*uu-one)
          sh(4,5,l)=four*rr*ss
          sh(1,5,l)=four*ss
          sh(2,5,l)=four*rr
          sh(3,5,l)=zero
          sh(4,6,l)=four*ss*tt
          sh(1,6,l)=zero
          sh(2,6,l)=four*tt
          sh(3,6,l)=four*ss
          sh(4,7,l)=four*rr*tt
          sh(1,7,l)=four*tt
          sh(2,7,l)=zero
          sh(3,7,l)=four*rr
          sh(4,8,l)=four*rr*uu
          sh(1,8,l)=four*(uu-rr)
          sh(2,8,l)=-(four*rr)
          sh(3,8,l)=-(four*rr)
          sh(4,9,l)=four*ss*uu
          sh(1,9,l)=-(four*ss)
          sh(2,9,l)=four*(uu-ss)
          sh(3,9,l)=-(four*ss)
          sh(4,10,l)=four*tt*uu
          sh(1,10,l)=-(four*tt)
          sh(2,10,l)=-(four*tt)
          sh(3,10,l)=four*(uu-tt)
  c*          if(i.le.4) then
  c*            sh(4,i,l)=rr*(two*rr-one)+ss*(two*ss-one)+tt*(two*tt-one)+
  c*     &       uu*(two*uu-one)
  c*            sh(1,i,l)=r(i)*(four*rr-one)-u(i)*(four*uu-one)
  c*            sh(2,i,l)=s(i)*(four*ss-one)-u(i)*(four*uu-one)
  c*            sh(3,i,l)=t(i)*(four*tt-one)-u(i)*(four*uu-one)
  c*          else
  c*            sh(4,i,l)=max(rr,one)*max(ss,one)*max(tt,one)*max(uu,one)
  c*            drr=eight*r(i)*max(ss,one)*max(tt,one)*max(uu,one)
  c*            dss=eight*s(i)*max(rr,one)*max(tt,one)*max(uu,one)
  c*            dtt=eight*t(i)*max(rr,one)*max(ss,one)*max(uu,one)
  c*            duu=(-eight)*u(i)*max(rr,one)*max(ss,one)*max(tt,one)
  c*            sh(1,i,l)=drr-duu
  c*            sh(2,i,l)=dss-duu
  c*            sh(3,i,l)=dtt-duu
  c*          end if
  c*        end do
  cdebug        dbsum=zero
  cdebug        do idb=1,nen
  cdebug          write(6,*) "sh:",(sh(jdb,idb,l),jdb=1,4)
  
negative.code:       end
  c
  c version
  c $Id: pquadtet.f,v 1.4 2004/08/12 02:19:36 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  