anchor.diff:            sh(3,i,l)=dtt*rr*ss*uu+dut*rr*ss*tt
           end do
         end do
  +cdebug      do idb=1,ngauss
  +cdebug        do jdb=1,nen
  +cdebug          write(6,*) "sh:",(sh(kdb,jdb,idb),kdb=1,4)
  +cdebug        end do
  +cdebug      end do
         call dcopy(nshsize,sh,ione,shj,ione)
   c
         return
         end
   c
   c version
  -c $Id: pquadhex.f,v 1.3 2004/07/07 19:40:22 willic3 Exp $
  +c $Id: pquadhex.f,v 1.4 2004/08/02 21:20:44 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code:       dhquadr(ri,si,ti)=abs(ri*si*ti)*ri
        dhquads(ri,si,ti)=abs(ri*si*ti)*si
        dhquadt(ri,si,ti)=abs(ri*si*ti)*ti
  c
  c...  definitions
  c
  
negative.code:       if(intord.ne.2) then
          ngauss=27
          g1=sqrt(three/five)
          w(1)=five/nine
          w(2)=eight/nine
          w(3)=five/nine
          l=0
          do l3=1,3
            tt=dble(l3-itwo)
            ttw=w(l3)
            do l2=1,3
              ss=dble(l2-itwo)
              ssw=w(l2)
              do l1=1,3
                rr=dble(l1-itwo)
                rrw=w(l1)
                l=l+1
                gauss(1,l)=rr*g1
                gauss(2,l)=ss*g1
  