anchor.diff:  c     cartesian natural coordinate along the length of the wedge.
   c
         call pquadwedge(sh(1,1,1,60),shj(1,1,1,60),gauss(1,1,60),
  -     & ngauss(60),nen(60),nenmax,ngaussmax,intord)
  +     & infetype(1,60),intord)
   c
   c...  Type 61:  quadratic hex with 9 points collapsed to a point
   c     (pyramid).  r and s are cartesian natural coordinates on the base,
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine preshape(sh,shj,gauss,ngauss,nen,nec,nee,ndof,nsd,
       & nenmax,ngaussmax,netypes,intord)
  c
  c...program to compute the shape functions in natural coordinates.
  c   The program also assigns integration points and weights.
  
negative.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),
       & infetype(1,29),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
  