anchor.diff:        close(kw)
   c
   c...  read face BC codes to determine which nodes are associated
  -c     with each boundary and fault.
  +c     with each boundary, fault, and Winkler BC.
   c
         read(kr,*) nbcfac
         do i=1,nbcfac
  
positive.code: c...  parameters
  c
        integer nsd,maxnodes,maxelmts,nen,maxnbrs,maxflts,maxfnodes
        integer maxbnds,maxmatf,ietyp,inf
        parameter(nsd=3,maxnodes=500000,maxelmts=500000,nen=4,
       & maxnbrs=100,maxflts=6,maxfnodes=100000,
       & maxbnds=30,maxmatf=5,ietyp=5,inf=0)
        integer kti,kto,kr,kw
        parameter(kti=5,kto=6,kr=10,kw=11)
        double precision eps
  
negative.code:         read(kr,*) (fsplit(j,1,i),j=1,nsd),(ifmat(1,i,j),j=1,nmatf(1,i))
          call pskip(kr)
          read(kr,*) (fsplit(j,2,i),j=1,nsd),(ifmat(2,i,j),j=1,nmatf(2,i))
        end do
  c
  c...  Winkler BC definitions
  c
        call pskip(kr)
        read(kr,*) nwsets
        do i=1,nwsets
          call pskip(kr)
          read(kr,*) iwcode(i),iwdir(i),modew(i),rhow(i),gw(i)
        end do
        close(kr)
  c
  