anchor.diff:  c
   c...  open input file
   c
  -      ierr=0
  -      nwinke=0
  -      nwink=0
  +      ierr=izero
  +      nwinke=izero
  +      nwink=izero
         open(kr,file=wfile,status="old",err=10)
   c
   c... scan the file, counting the number of entries.
  
positive.code:  40   continue
          read(kr,*,end=10,err=30) n,(iwink(j),j=1,ndof),
       &   (wink(j),j=1,ndof)
          nwinke=nwinke+1
          nnz=0
          do j=1,ndof
            if(iwink(j).ne.0) nnz=nnz+1
          end do
          nwink=nwink+nnz
          go to 40
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine scan_wink(nwinke,nwink,kr,wfile,ierr,errstrng)
  c
  c...  subroutine to perform an initial scan of the winkler BC
  c     file to determine the number of winkler BC.
  