anchor.diff:            gauss(4,l)=fourth*tetvol
           end do
         end if
  +c
  +      infetype(1)=ngauss
  +      infetype(2)=nen
  +      infetype(3)=nec
  +      infetype(4)=nee
  +c
         do l=1,ngauss
           do i=1,nen
             rr=r(i)*gauss(1,l)
  
positive.code: c...  Quadratic tet definition
  c
        nen=10
        if(intord.eq.2) then
          ngauss=ione
          gauss(1,1)=fourth
  
negative.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
  