anchor.diff:        end
   c
   c version
  -c $Id: plinwedge.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  +c $Id: plinwedge.f,v 1.2 2004/07/07 19:34:02 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c
        nen=isix
        ngauss=ione
        gauss(1,1)=zero
        gauss(2,1)=zero
        gauss(3,1)=zero
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine plinwedge(sh,shj,gauss,infetype,intord)
  c
  c... Subroutine to compute shape functions in natural coordinates,
  c    integration points, and weights for a linear wedge.
  c
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nshape.inc"
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  subroutine arguments
  c
        integer intord
        integer infetype(4)
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision gauss(nsd+1,ngaussmax)
  c
  c...  local constants
  c
        double precision r(6),s(6),t(6),u(6)
        data r/ 1d0, 0d0, 0d0, 1d0, 0d0, 0d0/
  