anchor.diff:        i1=nnblnk(fileroot)
         i2=nchar(fileroot)
         call ifill(ibcnode,izero,maxnodes*maxbnds)
  -      call ifill(nbcnodes,izero,maxbnds)
  +      call ifill(iwbcnode,izero,maxnodes*maxwsets)
  +      call fill(wink,zero,maxnodes*maxwsets)
   c
   c...  read parameter file
   c
  
positive.code: c
        integer icflip(4)
        data icflip/1,3,2,4/
        integer kwb,kwc,kwfb(maxflts),kwfc(maxflts)
        data kwb/12/
        data kwc/13/
        data kwfb/12,13,14,15,16,17/
        data kwfc/18,19,20,21,22,23/
        integer izero
        data izero/0/
  c
  c...  parameters read from parameter file
  c
        integer nbc,iconopt
        integer ibcode(maxbnds),ibc(nsd,maxbnds),isn(nsd)
        integer iftype(maxflts),nmatf(2,maxflts),ifmat(2,maxflts,maxmatf)
        double precision bc(nsd,maxbnds),fsplit(nsd,2,maxflts)
        double precision cscale
  c
  c...  filenames
  c
        character fileroot*200,pfile*200,ifile*200,nfile*200,cfile*200
        character bcfile*200,bccfile*200
        character fbcfile*200,fbccfile*200
  c
  c...  parameters and variables read from netgen file
  
negative.code:             end if
            end do
          end if
          do j=1,nwsets
            if(ibcfac.eq.iwcode(j)) then
              area=acomp(x,ibcvert,nsd,numnp,iwdir(j))
              do k=1,3
                iwbcnode(ibcvert(k),j)=ibcfac
                wink(ibcvert(k),j)=wink(ibcvert(k),j)+
       &         area*rhow(j)*gw(j)/3.0d0
              end do
            end if
          end do
        end do
        close(kr)
        nfltnodes=0
        do i=1,numnp
          if(ifltnode(i).eq.1) then
  