anchor.diff:          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)
  +        call shapql(rgs,xl,det,sh,ien,nen,infin,iopt,n,idout,kto,kw)
   c
   c...  compute surface element for use in integration
   c
  
positive.code:       end
  c
  c version
  c $Id: presurql.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:       double precision g,vmag,ds,det
        double precision sh(4,8),ptmp(3,4),v1(3),v2(3),vc(3),rgs(3)
  c
  cdebug      write(6,*) "Hello from presurql_f!"
  c
        g=root3i
  c
  