anchor.diff:  c...  parameters
   c
         integer nsd,maxnodes,maxelmts,nen,maxnbrs,maxflts,maxfnodes
  -      integer nsides,maxmatf,ietyp,inf
  +      integer maxbnds,maxmatf,ietyp,inf
         parameter(nsd=3,maxnodes=500000,maxelmts=500000,nen=4,
  -     & maxnbrs=100,maxflts=6,maxfnodes=10000,nsides=6,maxmatf=5,
  -     & ietyp=5,inf=0)
  +     & 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
  
positive.code: c
        integer icflip(4)
        data icflip/1,3,2,4/
        integer kwb(nsides),kwc(nsides),kwfb(maxflts),kwfc(maxflts)
        data kwb/12,13,14,15,16,17/
        data kwc/18,19,20,21,22,23/
        data kwfb/12,13,14,15,16,17/
        data kwfc/18,19,20,21,22,23/
  c
  c...  parameters read from parameter file
  c
        integer ibc(nsd,nsides),isn(nsd)
        integer iftype(maxflts),nmatf(2,maxflts),ifmat(2,maxflts,maxmatf)
        double precision xlim(2,nsd),bc(nsd,nsides),fsplit(nsd,2,maxflts)
        double precision cscale
  c
  c...  filenames
  
negative.code:       read(kr,*) cscale,(idir(i),i=1,nsd)
        if(cscale.eq.0.0d0) cscale=1.0d0
  c
  c...  bc info
  c
        call pskip(kr)
        read(kr,*) nbc
        do i=1,nbc
          call pskip(kr)
          read(kr,*) ibcode(i),(ibc(j,i),bc(j,i),j=1,nsd)
        end do
  c...  connectivity order option
  c
  