anchor.diff:  c
    20   continue
           ierr=1
  +        errstrng="read_wink"
           close(kr)
           return
   c
  
positive.code: c
  c...  open output file
  c
        if(idout.gt.0) open(kw,file=ofile,status="old",access="append")
  c
  c...  open input file
  c
        ierr=0
        if(nwink.eq.0) return
        open(kr,file=wfile,status="old",err=20)
  c
  c.......read winkler force data and output results, if desired
  c
        call pskip(kr)
        nwtot=0
        nlines=0
        npage=50
        do i=1,nwinke
          read(kr,*,err=30,end=30) n,(iwtmp(j),j=1,ndof),
       &   (wtmp(j),j=1,ndof)
          nnz=0
          do j=1,ndof
            wtmp(j)=wscal(j)*wtmp(j)
            if(iwtmp(j).ne.0) then
              nnz=nnz+1
              nwtot=nwtot+1
              iwink(1,nwtot)=iwtmp(j)
  
negative.code:             wink(nwtot)=wtmp(j)
            end if
          end do
          if(idout.gt.izero.and.nnz.ne.izero) then
            if(mod(nlines,npage).eq.izero) then
              write(kw,1000,err=50)
              write(kw,3000,err=50) (labeld(iline),iline=1,ndof)
              write(kw,2005,err=50)
            end if
            nlines=nlines+1
            write(kw,2020,err=50) n,(wtmp(iline),iwtmp(iline),
       &     iline=1,ndof)
          end if 
        end do
        close(kr)
        if(idout.gt.izero) close(kw)
  c
  c...  normal return
  c
  