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:           rr=gquad(gauss(1,l),r(i))
            ss=gquad(gauss(2,l),s(i))
            tt=gquad(gauss(3,l),t(i))
            drr=dgquad(gauss(1,l),r(i))
            dss=dgquad(gauss(2,l),s(i))
            dtt=dgquad(gauss(3,l),t(i))
            sh(4,i,l)=rr*ss*tt
            sh(1,i,l)=drr*ss*tt
            sh(2,i,l)=dss*rr*tt
            sh(3,i,l)=dtt*rr*ss
          end do
          do i=1,4
            sh(i,1,l)=sh(i,1,l)-half*(sh(i,9,l)+sh(i,12,l)+ sh(i,17,l))
            sh(i,2,l)=sh(i,2,l)-half*(sh(i,9,l)+sh(i,10,l)+sh(i,18,l))
            sh(i,3,l)=sh(i,3,l)-half*(sh(i,10,l)+sh(i,11,l)+sh(i,19,l))
            sh(i,4,l)=sh(i,4,l)-half*(sh(i,11,l)+sh(i,12,l)+sh(i,20,l))
            sh(i,5,l)=sh(i,5,l)-half*(sh(i,13,l)+sh(i,16,l)+sh(i,17,l))
            sh(i,6,l)=sh(i,6,l)-half*(sh(i,13,l)+sh(i,14,l)+sh(i,18,l))
            sh(i,7,l)=sh(i,7,l)-half*(sh(i,14,l)+sh(i,15,l)+sh(i,19,l))
            sh(i,8,l)=sh(i,8,l)-half*(sh(i,15,l)+sh(i,16,l)+sh(i,20,l))
          end do
        end do
        call dcopy(nshsize,sh,ione,shj,ione)
  
negative.code:       end
  c
  c version
  c $Id: pquadhex.f,v 1.2 2004/07/06 19:14:36 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  