anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine getshapn(xl,sh,shj,shd,shbar,det,gauss,vol,iel,nen,nsd,
  -     & ngauss,ierr)
  +      subroutine getshapn(xl,sh,shj,shd,shbar,det,gauss,vol,iel,nen,
  +     & ngauss,ierr,errstrng)
   c
   c...  Subroutine to compute shape functions and derivatives at
   c     Gauss points.
  
positive.code: c
        include "implicit.inc"
  c
  c...  dimension parameters
  c
        include "nshape.inc"
  c
  c...  subroutine arguments
  c
        integer iel,nen,nsd,ngauss,ierr
        double precision xl(nsd,nen),sh(nsd+1,nenmax,ngaussmax)
        double precision shj(nsd+1,nenmax,ngaussmax)
        double precision shd(nsd+1,nenmax,ngaussmax),shbar(nsd+1,nenmax)
        double precision det(ngauss),gauss(nsd+1,ngaussmax),vol
  c
  c...  defined constants
  c
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  local variables
  c
  cdebug      integer idb,jdb
        integer l
        double precision xs(3,3)
  c
  
negative.code:       end
  c
  c version
  c $Id: getshapn.f,v 1.4 2004/07/05 20:04:32 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  