anchor.diff:        end
   c
   c version
  -c $Id: read_wink.f,v 1.4 2005/04/14 00:59:45 willic3 Exp $
  +c $Id: read_wink.f,v 1.5 2005/04/16 00:45:25 willic3 Exp $
   c
   c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
   c
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine read_wink(wink,wscal,iwink,id,numnp,nwink,nwinke,
       & kr,wfile,ierr,errstrng)
  c
  c....program for reading and data on winkler restoring forces
  c
  c          wink(ndof,numnp) = values of winkler restoring spring
  c                             constant, force(i,j)=-wink(i,j)*disp(i,j)
  c
  c          iwtmp(ndof,numnp) = application mode:
  c                               iwink = 0, no winkler forces,
  c                               iwink = 1, applied throuthout computation
  c                               iwink = -n, uses load history factor n
  c
  c        After input, the iwtmp array is transformed into the iwink
  c        array, which has the following components:
  c          iwink(1,i)       = nonzero application mode from iwtmp
  c          iwink(2,i)       = global equation number to which force is
  c                             applied
  c
  c     Error codes:
  c         0:  No error
  
negative.code: c.......read winkler force data and output results, if desired
  c
        call pskip(kr)
        do i=1,nwinke
          read(kr,*,err=30,end=30) iwinkid(i),(iwinkdef(j,i),j=1,ndof),
       &   (winkdef(j,i),j=1,ndof)
          do j=1,ndof
            winkdef(j,i)=wscal(j)*winkdef(j,i)
          end do
        end do
        close(kr)
  