anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine pquadpyr(sh,shj,gauss,ngauss,nen,nsd,nenmax,ngaussmax,
  -     & intord)
  +      subroutine pquadpyr(sh,shj,gauss,infetype,intord)
   c
   c... Subroutine to compute shape functions in natural coordinates,
   c    integration points, and weights for a quadratic pyramid.
   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 nsd,nenmax,ngaussmax,intord
  -      integer ngauss,nen
  +      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...  defined constants
  -c
  -      include "nconsts.inc"
  -      include "rconsts.inc"
  +c...  local constants
   c
         double precision c1,c2,c3,c4,c5
         parameter(c1 =  15.0d0,
  
positive.code: c
  c...  local variables
  c
        integer i,l,nshsize,ngssize
        double precision q(9),h(3),dq(9,2),dh(3)
        double precision rr,ss,tt,g,w
  c
  
negative.code:         h(1)=half*(tt*tt-tt)
          h(2)=half*(tt*tt+tt)
          h(3)=one-tt*tt
          dq(1,1)=(-fourth)*((one-ss)*(-rr-ss-one)+(one-rr)*(one-ss))
          dq(2,1)=( fourth)*((one-ss)*( rr-ss-one)+(one+rr)*(one-ss))
          dq(3,1)=( fourth)*((one+ss)*( rr+ss-one)+(one+rr)*(one+ss))
          dq(4,1)=(-fourth)*((one+ss)*(-rr+ss-one)+(one-rr)*(one+ss))
          dq(5,1)=(-rr)*(one-ss)
          dq(6,1)=( half)*(one-ss*ss)
          dq(7,1)=(-rr)*(one+ss)
          dq(8,1)=(-half)*(one-ss*ss)
          dq(1,2)=(-fourth)*((one-rr)*(-rr-ss-one)+(one-rr)*(one-ss))
          dq(2,2)=(-fourth)*((one+rr)*( rr-ss-one)+(one+rr)*(one-ss))
          dq(3,2)=( fourth)*((one+rr)*( rr+ss-one)+(one+rr)*(one+ss))
          dq(4,2)=( fourth)*((one-rr)*(-rr+ss-one)+(one-rr)*(one+ss))
          dq(5,2)=(-half)*(one-rr*rr)
          dq(6,2)=(-ss)*(one+rr)
          dq(7,2)=( half)*(one-rr*rr)
          dq(8,2)=(-ss)*(one-rr)
          dq(9,1)=zero
          dq(9,2)=zero
          dh(1)=half*(two*tt-one)
          dh(2)=half*(two*tt+one)
          dh(3)=(-two)*tt
          do i=1,nen
            sh(4,i,l)=q(iq(i))*h(ih(i))
            sh(1,i,l)=dq(iq(i),1)*h(ih(i))
  