anchor.diff:        include "ndimens.inc"
         include "nshape.inc"
         include "materials.inc"
  +      include "nconsts.inc"
  +      include "rconsts.inc"
   c
   c...  subroutine arguments
   c
  
positive.code: c...  included dimension and type statements
  c
        include "labels_dim.inc"
  c
  c...  defined constants
  c
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  intrinsic functions
  c
  
negative.code: c
  c...  output increments/rates, if desired.
  c
          if(istatout(2,i).ne.izero) then
            tmult=one
            if(istatout(2,i).eq.ione.and.delt.gt.zero) tmult=one/delt
            nout=izero
            do iel=1,numelt
              imat=infiel(2,iel)
              ietype=infiel(3,iel)
  