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
  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
        double precision uu,dur,dus,dut,v,vi,rp,ri,sp,si,tp,ti
  c
  
negative.code:           end do
          end do
        end if
  c
        infetype(1)=ngauss
        infetype(2)=nen
        infetype(3)=nec
        infetype(4)=nee
  c
        do l=1,ngauss
          do i=1,nen
            rr=gquad(gauss(1,l),r(i))
  