anchor.diff:        double precision g,vmag,ds,det
         double precision sh(4,8),ptmp(3,4),v1(3),v2(3),vc(3),rgs(3)
   c
  -c*      write(6,*) "Hello from presurql_f!"
  +cdebug      write(6,*) "Hello from presurql_f!"
   c
         g=root3i
   c
  
positive.code:         rgs(1)=rf(iiside)+ri(iiside)*g*rg(ipnt)
          rgs(2)=sf(iiside)+si(iiside)*g*sg(ipnt)
          rgs(3)=tf(iiside)+ti(iiside)*g*tg(ipnt)
          call shapql(rgs,xl,det,sh,ien,nen,nsd,infin,iopt,n,idout,kto,kw)
  c
  c...  compute surface element for use in integration
  c
  
negative.code:       end
  c
  c version
  c $Id: presurql.f,v 1.2 2004/08/12 02:22:17 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  