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:         end do
          gauss(3,2)=-root3i
        end if 
        do l=1,ngauss
          do i=1,nen
            rr=one-r(i)+r(i)*gauss(1,l)
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine plinwedge(sh,shj,gauss,infetype,intord)
  c
  c... Subroutine to compute shape functions in natural coordinates,
  c    integration points, and weights for a linear wedge.
  c
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nshape.inc"
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  subroutine arguments
  c
        integer intord
        integer infetype(4)
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision gauss(nsd+1,ngaussmax)
  c
  c...  local constants
  c
        double precision r(6),s(6),t(6),u(6)
        data r/ 1d0, 0d0, 0d0, 1d0, 0d0, 0d0/
  