c
c     declarations and common block for the coordinates of the
c     Gauss point [natural coordinates of a triangle which are
c     also the shape functions for a three-node triangular
c     element]
c
      real*8 xi1, xi2, xi3
      common /gausspt_ref/ xi1, xi2, xi3 ! xi1 + xi2 + xi3 = 1
