anchor.diff:  c...  Quadratic pyramid definition
   c
         nen=13
  +      nec=nsd*nen
  +      nee=ndof*nen
         if(intord.eq.2) then
           ngauss=ifive
           g=eight*sqrt(two/c1)/five
  
positive.code:         gauss(3,13)=g5
          gauss(4,13)=w4
        end if
        do l=1,ngauss
          rr=gauss(1,l)
          ss=gauss(2,l)
  
negative.code:       end
  c
  c version
  c $Id: pquadpyr.f,v 1.2 2004/07/07 19:48:45 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  