anchor.diff:        read(kr,*) cscale,(idir(i),i=1,nsd)
         if(cscale.eq.0.0d0) cscale=1.0d0
   c
  -c...  box dimensions and bc
  +c...  bc info
   c
         call pskip(kr)
  -      read(kr,*) ((xlim(j,i),j=1,2),i=1,nsd)
  -      do i=1,nsides
  +      read(kr,*) nbc
  +      do i=1,nbc
           call pskip(kr)
  -        read(kr,*) (ibc(j,i),bc(j,i),j=1,nsd)
  +        read(kr,*) ibcode(i),(ibc(j,i),bc(j,i),j=1,nsd)
         end do
   c...  connectivity order option
   c
  
positive.code:       program readnetgen
  c
  c...  quick and dirty code to translate netgen neutral format to a
  c     format suitable for lithomop.  The domain is assumed to be a
  c     rectangular solid.  Geometry assumptions:
  c             xmin = left
  c             xmax = right
  c             ymin = front
  c             ymax = back
  c             zmin = bottom
  c             zmax = top
  c
  c     At present, code is only set up for linear tetrahedra, although
  c     this would be easy to change.
  
negative.code: c
  c...  parameters and variables read from netgen file
  c
        integer numnp,numel,numflt,nbcfac,ibcflt,ibcfac
        integer ien(nen,maxelmts),mat(maxelmts),ibcvert(3)
        double precision x(nsd,maxnodes)
  c
  c...  external routines
  