anchor.diff:  c
   c...  local variables
   c
  -      integer i,l,nshsize,ngssize
  +      integer nen,ngauss,nec,nee,i,l,nshsize,ngssize
         double precision t(6),h(3),dt(6,2),dh(3)
         double precision rr,ss,tt,uu,g1,w1
         double precision trivol
  
positive.code: c
        nen=15
        ngauss=itwo
        do l=1,ngauss
          gauss(1,l)=third
          gauss(2,l)=third
  
negative.code:         dt(2,1)=zero
          dt(3,1)=one-four*uu
          dt(4,1)=four*ss
          dt(5,1)=(-four)*ss
          dt(6,1)=four*(uu-rr)
          dt(1,2)=zero
          dt(2,2)=four*ss-one
          dt(3,2)=one-four*uu
          dt(4,2)=four*rr
          dt(5,2)=four*(uu-ss)
          dt(6,2)=(-four)*rr
          dh(1)=half*(two*tt+one)
          dh(2)=half*(two*tt-one)
          dh(3)=one-two*tt
  