anchor.diff:  c
   c...  local variables
   c
  -      integer i,l,nshsize,ngssize
  +      integer nen,ngauss,nec,nee,i,l,nshsize,ngssize
         double precision rr,ss,tt
   c
   c...  definitions
  
positive.code: c
        nen=ieight
        ngauss=ione
        gauss(1,1)=zero
        gauss(2,1)=zero
        gauss(3,1)=zero
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine plinhex(sh,shj,gauss,infetype,intord)
  c
  c... Subroutine to compute shape functions in natural coordinates,
  c    integration points, and weights for a trilinear hexahedron.
  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(8),s(8),t(8)
        data r/-1d0, 1d0, 1d0,-1d0,-1d0, 1d0, 1d0,-1d0/
  