anchor.diff:  c
   c...  local variables
   c
  -      integer i,nshsize,ngssize
  +      integer nen,ngauss,nec,nee,i,nshsize,ngssize
         double precision rr,ss,tt,uu
         double precision tetvol
   c
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine plintet(sh,shj,gauss,ngauss,nen,nsd,nenmax,ngaussmax,
       & intord)
  c
  c... Subroutine to compute shape functions in natural coordinates,
  c    integration points, and weights for a linear tetrahedron.
  c
        include "implicit.inc"
  c
  c...  subroutine arguments
  c
        integer nsd,nenmax,ngaussmax,intord
        integer ngauss,nen
        double precision sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision gauss(nsd+1,ngaussmax)
  c
  c...  defined constants
  c
        include "nconsts.inc"
        include "rconsts.inc"
  c
        double precision r(4),s(4),t(4),u(4)
        data r/ 1d0, 0d0, 0d0, 0d0/
  
negative.code:       end
  c
  c version
  c $Id: plintet.f,v 1.2 2004/07/07 19:28:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  