anchor.diff:   40   continue
           read(kr,*,end=10,err=30) n,(iwinkx(j),j=1,ndof),
        &   (winkx(j),j=1,ndof)
  -        nwinkxe=nwinkxe+1
  -        nnz=0
  +        nwinkxe=nwinkxe+ione
  +        nnz=izero
           do j=1,ndof
  -          if(iwinkx(j).ne.0) nnz=nnz+1
  +          if(iwinkx(j).ne.0) nnz=nnz+ione
           end do
           nwinkx=nwinkx+nnz
           go to 40
  
positive.code: c
   30   continue
          ierr=3
          close(kr)
          return
  c
        end
  c
  c version
  c $Id: scan_winkx.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine scan_winkx(numslp,nwinkxe,nwinkx,kr,wxfile,ierr,
       & errstrng)
  c
  c...  subroutine to perform an initial scan of the differential winkler
  c     BC file to determine the number of differential winkler BC.
  