anchor.diff:  c
   c...  parameters and variables read from netgen file
   c
  -      integer numnp,numel,numflt,iconopt,nbcfac,ibcflt,ibcfac
  -      integer ien(nen,maxelmts),mat(maxelmts)
  +      integer numnp,numel,numflt,nbcfac,ibcflt,ibcfac
  +      integer ien(nen,maxelmts),mat(maxelmts),ibcvert(3)
         double precision x(nsd,maxnodes)
   c
   c...  external routines
  
positive.code:       read(kti,"(a200)") fileroot
        i1=nnblnk(fileroot)
        i2=nchar(fileroot)
  c
  c...  read parameter file
  c
  
negative.code: c
  c...  local variables
  c
        double precision det,sgn,xl(nsd,nen),xtmp(nsd)
        integer iadjf(maxflts*maxfnodes,maxnbrs),inodef(maxflts*maxfnodes)
        integer ientmp(nen),nelsf(maxflts*maxfnodes),ifltnode(maxnodes)
        integer ifhist(maxflts),idir(nsd),ibcnode(maxnodes,maxbnds)
        integer nbcnodes(maxbnds)
        integer i,j,k,l,jj,i1,i2,j1,j2,nenl,nsdl,nfltnodes,kk
        integer iel,nflip
        character cstring*15,cbound(maxbnds)*3,cfnum*2
        data cstring/"coord_units = m"/
        data cbound/
       & "b01","b02","b03","b04","b05","b06","b07","b08","b09","b10",
       & "b11","b12","b13","b14","b15","b16","b17","b18","b19","b20",
       & "b21","b22","b23","b24","b25","b26","b27","b28","b29","b30"/
  c
        nenl=nen
        nsdl=nsd
  