anchor.diff:        if(intord.ne.2) then
           ngauss=27
           g1=sqrt(three/five)
  -        w1=five/nine
  -        w2=eight/nine
  +        w(1)=five/nine
  +        w(2)=eight/nine
  +        w(3)=five/nine
           l=0
           do l3=1,3
  -          tt=dble(l3-2)
  -          ttw=tt*w1+(abs(tt)-one)*w2
  +          tt=dble(l3-itwo)
  +          ttw=w(l3)
             do l2=1,3
  -            ss=dble(l2-2)
  -            ssw=ss*w1+(abs(ss)-one)*w2
  +            ss=dble(l2-itwo)
  +            ssw=w(l2)
               do l1=1,3
  -              rr=dble(l1-2)
  -              rrw=rr*w1+(abs(rr)-one)*w2
  +              rr=dble(l1-itwo)
  +              rrw=w(l1)
                 l=l+1
                 gauss(1,l)=rr*g1
                 gauss(2,l)=ss*g1
  
positive.code:           end do
          end do
        end if
  c
        infetype(1)=ngauss
        infetype(2)=nen
  
negative.code: c...  local variables
  c
        integer nen,ngauss,nec,nee,i,l,l1,l2,l3,nshsize,ngssize
        double precision g1,rr,ss,tt,rrw,ssw,ttw,drr,dss,dtt
        double precision uu,dur,dus,dut,v,vi,rp,ri,sp,si,tp,ti
        double precision w(3)
  cdebug      integer idb,jdb,kdb
  c
  c...  function definitions
  c
  