anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine pquadhex(sh,shj,gauss,ngauss,nen,nsd,nenmax,ngaussmax,
  -     & intord)
  +      subroutine pquadhex(sh,shj,gauss,ngauss,nen,intord)
   c
   c... Subroutine to compute shape functions in natural coordinates,
   c    integration points, and weights for a quadratic (20-node)
  
positive.code: c
  c...  user-defined functions
  c
        double precision gquad,dgquad
  c
  c...  local variables
  c
        integer i,l,l1,l2,l3,nshsize,ngssize
        double precision g1,w1,w2,rr,ss,tt,rrw,ssw,ttw,drr,dss,dtt,beta
        double precision betai
  c
  c...  function definitions
  c
        gquad(beta,betai)=(half*(one-abs(betai))+half)*(one+beta*betai+
       & (abs(betai)-one)*beta*beta)
        dgquad(beta,betai)=half*betai-two*beta*(one-abs(betai))
  c
  c...  definitions
  c
  
negative.code: 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 ngauss,nen,intord
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision gauss(nsd+1,ngaussmax)
  c
  c...  local data definitions
  c
        double precision r(20),s(20),t(20)
        data r/-1d0, 1d0, 1d0,-1d0,-1d0, 1d0, 1d0,-1d0,
  