anchor.diff:  c
         include "implicit.inc"
   c
  +c...  parameter definitions
  +c
  +      include "ndimens.inc"
  +      include "nconsts.inc"
  +      include "rconsts.inc"
  +c
   c...  subroutine arguments
   c
  -      integer ndof,numnp,nen,nee,nskdim
  +      integer numnp,nen,nee
         integer ien(nen)
         double precision s(nee,nee),stemp(nee,nee),skew(nskdim,numnp)
  -c
  -c...  defined constants
  -c
  -      include "nconsts.inc"
  -      include "rconsts.inc"
   c
   c...  local variables
   c
  
positive.code:       do i=1,nen
          k=ien(i)
          if((skew(1,k).ne.zero).and.(skew(nskdim,k).ne.zero)) then
            call formrt(skew(1,k),rot,nskdim)
  	  call dcopy(nee*nee,s,ione,stemp,ione)
            ll=ndof*(i-1)
  c***********  this part still needs to be fixed for BLAS.
  
negative.code:       end
  c
  c version
  c $Id: rstiff.f,v 1.2 2004/07/05 20:12:11 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  