anchor.diff:  c
         include "implicit.inc"
   c
  +c...  parameter definitions
  +c
  +      include "ndimens.inc"
  +      include "rconsts.inc"
  +c
   c...  subroutine arguments
   c
  -      integer nen,nsd
  +      integer nen
         double precision det,sh(nsd+1,nen),shd(nsd+1,nen),xs(nsd,nsd)
  -c
  -c...  defined constants
  -c
  -      include "rconsts.inc"
   c
   c...  local variables
   c
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine getder(det,sh,shd,xs,nen,nsd)
  c
  c...  subroutine to compute shape function derivatives in
  c     global coordinates
  
negative.code: c
  cdebug      write(6,*) "Hello from getder_f!"
  c
  c
  c...transform natural derivatives to (x,y,z) derivatives using an
  c   explicit 3x3 matrix inversion routine
  