anchor.diff:  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=infmat(1,imat)
  +        matm=ifamilyv(2,imat)
   c
  -        if(ietypei.le.izero.or.ietypei.gt.netypesi) then
  +        if(ietypei.ne.ietypev) then
             ierr=106
             errstrng="scan_connect"
             return
  
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:           return
          end if
  c
          numelv=numelv+ione
          ifamilyv(1,imat)=ifamilyv(1,imat)+ione
  c
          go to 40
  c
  