anchor.diff:  c
         nen=isix
         ngauss=ione
  +      nec=nsd*nen
  +      nee=ndof*nen
         gauss(1,1)=zero
         gauss(2,1)=zero
         gauss(3,1)=zero
  
positive.code: c
  c...  local variables
  c
        integer i,l,nshsize,ngssize
        double precision rr,ss,tt,uu,drr,dss,dtt
        double precision trivol
  c
  
negative.code:         end do
          gauss(3,2)=-root3i
        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=one-r(i)+r(i)*gauss(1,l)
  