anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine read_wink(wink,wscal,iwink,id,numnp,nwink,nwinke,
  -     & kr,wfile,ierr,errstrng)
  +      subroutine read_wink(winkdef,wscal,iwinkdef,iwinkid,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          winkdef(ndof,nwinke) = values of winkler restoring spring
   c                             constant, force(i,j)=-wink(i,j)*disp(i,j)
   c
  -c          iwtmp(ndof,numnp) = application mode:
  +c          iwinkdef(ndof,nwinke) = 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
  
positive.code: c
        include "ndimens.inc"
        include "nconsts.inc"
  c
  c...  subroutine arguments
  c
        integer numnp,nwink,nwinke,kr,ierr
        integer iwink(2,nwink),id(ndof,numnp)
        double precision wink(nwink),wscal(3)
        character wfile*(*),errstrng*(*)
  c
  c...  included dimension and type statements
  
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)
  