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: c...  local variables
  c
        integer nen,ngauss,nec,nee,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
  c...  function definitions
  c
  
negative.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
  cdebug      write(6,*) "Hello from pquadhex_f!"
  c
  c
  c...  definitions
  c
  