anchor.diff:  c                                 is identical to sh.
   c       gauss(nsd+1,ngaussmax,netypes)     = Gauss point coordinates and
   c                                            weights
  -c       ngauss(netypes)                    = number of Gauss points
  -c       nen(netypes)                       = number of element nodes
  -c       nec(netypes)                       = number of element
  -c                                            coordinate directions
  -c       nee(netypes)                       = number of element equations
  -c
  +c       infetype(4,netypes)                = element type info array
  +c         infetype(1,i) = number of gauss points for element type
  +c         infetype(2,i) = number of nodes for element type
  +c         infetype(3,i) = number of element coordinate directions for
  +c                         element type
  +c         infetype(4,i) = number of element degrees of freedom for
  +c                         element type
   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 nsd,nenmax,ngaussmax,netypes,intord
  -      integer ngauss(netypes),nen(netypes)
  +      integer intord
  +      integer infetype(4,netypes)
         double precision sh(nsd+1,nenmax,ngaussmax,netypes)
         double precision shj(nsd+1,nenmax,ngaussmax,netypes)
         double precision gauss(nsd+1,ngaussmax,netypes)
   c
  -c...  defined constants
  -c
  -      include "nconsts.inc"
  -      include "rconsts.inc"
  +c...  local constants
   c
         double precision rlh(8),slh(8),tlh(8)
         data rlh/-1d0, 1d0, 1d0,-1d0,-1d0, 1d0, 1d0,-1d0/
  
positive.code: c     cartesian natural coordinate along the length of the wedge.
  c
        call plinwedge(sh(1,1,1,29),shj(1,1,1,29),gauss(1,1,29),
       & ngauss(29),nen(29),nenmax,ngaussmax,intord)
  c
  c...  Type 30:  linear hex with 4 points collapsed to a point (pyramid)
  c     r and s are cartesian natural coordinates on the base, and t is
  
negative.code: c
  c...  Type 59:  quadratic hex with one set of collapsed nodes (wrick)
  c
        infetype(1,59)=infetype(1,32)
        infetype(2,59)=18
        infetype(3,59)=nsd*infetype(2,59)
        infetype(4,59)=ndof*infetype(2,59)
        call dcopy(ngssize,gauss(1,1,32),ione,gauss(1,1,59),ione)
        call dcopy(nshsize,sh(1,1,1,32),ione,sh(1,1,1,59),ione)
        do l=1,infetype(1,59)
          do i=1,4
            sh(i,7,l,59)=sh(i,7,l,59)+sh(i,8,l,59)+sh(i,15,l,59)
            sh(i,8,l,59)=zero
  