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: c...  definitions
  c
        tetvol=sixth
        nshsize=(nsd+1)*nenmax*ngaussmax
        ngssize=(nsd+1)*ngaussmax
  c
  
negative.code:         gauss(2,1)=fourth
          gauss(3,1)=fourth
          gauss(4,1)=tetvol
  cdebug        gauss(1,1)=fourth-0.01d0
  cdebug        gauss(2,1)=fourth-0.01d0
  cdebug        gauss(3,1)=fourth-0.01d0
        else
          ngauss=ifour
          g1=(five-sqrt(five))/20.0d0
  