anchor.diff:  c...  open input file
   c
         ierr=izero
  -      numelt=izero
  -      nconsz=izero
  -      call ifill(infmat,izero,ithree*numat)
  +      numelv=izero
  +      ietypev=izero
  +      call ifill(ifamilyv,izero,ithree*nvfamilies)
         open(kr,file=cfile,status="old",err=20)
   c
  -c...  set up part of the infmat array that has already been determined
  +c...  set up part of the ifamilyv array that has already been determined
   c     by python material property input routine.
   c
  -      do i=1,numat
  -        infmat(1,i)=matmod(i)
  -        infmat(3,i)=indprop(i)
  +      do i=1,nvfamilies
  +        ifamilyv(2,i)=matmod(i)
  +        ifamilyv(3,i)=indprop(i)
         end do
   c
   c... scan the file, counting the number of entries for each type of
  
positive.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine scan_connect(neni,infmat,infmatmod,matmod,
       & indprop,numat,numelt,nconsz,kr,cfile,ierr,errstrng)
  c
  c...  subroutine to perform an initial scan of the element connectivity
  c     file to determine the total number of elements and the number of
  c     elements of each material type.
  c
  c     Note that the variable netypesi (and the associated arrays) only
  c     defines the number of 'primitive' element types, and does not
  
negative.code: c    comment.
  c
        call pskip(kr)
        read(kr,*,end=10,err=30) n,ietypei,imat,infin,
       & (ien(j),j=1,neni(ietypei))
        ietypev=ietypei
        backspace(kr)
        if(ietypev.le.izero.or.ietypev.gt.netypesi) then
          ierr=106
          errstrng="scan_connect"
          return
        end if
   40   continue
          read(kr,*,end=10,err=30) n,ietypei,imat,infin,
       &   (ien(j),j=1,neni(ietypei))
  c
          matm=ifamilyv(2,imat)
  c
          if(ietypei.ne.ietypev) then
            ierr=106
            errstrng="scan_connect"
            return
  