anchor.diff:        end
   c
   c version
  -c $Id: rstiff.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  +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
  
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: c
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  subroutine arguments
  c
        integer numnp,nen,nee
        integer ien(nen)
        double precision s(nee,nee),stemp(nee,nee),skew(nskdim,numnp)
  c
  c...  local variables
  c
  