anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine pquadhex(sh,shj,gauss,ngauss,nen,intord)
  +      subroutine pquadhex(sh,shj,gauss,infetype,intord)
   c
   c... Subroutine to compute shape functions in natural coordinates,
   c    integration points, and weights for a quadratic (20-node)
  
positive.code: c
        nen=20
        ngauss=ieight
        do l=1,ngauss
          gauss(1,l)=r(l)*root3i
          gauss(2,l)=s(l)*root3i
  
negative.code: 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(20),s(20),t(20)
        data r/-1d0, 1d0, 1d0,-1d0,-1d0, 1d0, 1d0,-1d0,
  